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A Microflash Unit for Ballistic Photography* 


W. W. McCormick, L. MApANsky,** AND A. F. FAIRBANKS*** 
Randall Laboratory of Physics, University of Michigan, Ann Arbor, Michigan 
(Received August 19, 1947) 


A flash unit has been developed for photographing large projectiles in flight. The light source 
combines high intensity with short duration. Details of the geometry and processing of the 
lamp are given, together with a discussion of the methods of measuring intensity and duration. 
Several of the units described have been in use at the Aberdeen Proving Ground and some 


photographs taken there are reproduced. 





HE condensed spark has long been used as 

a source in ballistic photography to obtain 
shadow photographs of small projectiles in flight. 
For projectiles larger than, say, three inches in 
diameter, however, the area of film required 
becomes excessive since no lenses are used, and 
the field of illumination is insufficient. These 
limitations may be overcome by substituting 
for the spark, which is essentially a point source 
of light, an extended source such as a flash lamp. 
The projectile is then photographed by light 
reflected from the surface into one or more 
cameras advantageously placed about the tra- 
jectory. In addition to furnishing information 
regarding the orientation in space of the pro- 
jectile, this technique permits examination of the 
surface for flaws which sometimes develop during 
firing. It does not, however, yield any informa- 
tion regarding shock wave patterns as does the 
condensed spark source. 


* This paper is based on work completed in 1944 under 
Contract NDCrc-185 between the University of Michigan 
and the Office of Scientific Research and Development and 
submitted for clearance in 1946. 

** National Research Fellow. 

*** Now at the North American Aviation Company. 


At the time that the work reported in this 
paper was begun, some reflection photographs 
using a commercial flash unit had been taken at 
the Aberdeen Proving Ground. As a result of 
these tests it was decided that further investiga- 
tion of the flash lamp source, with particular 
emphasis on reduction of the duration of the 
flash, was desirable. This was undertaken at the 
Randall Laboratory of Physics, University of 
Michigan, in cooperation with the Ballistic 
Research Laboratory at Aberdeen. 


THE LIGHT SOURCE 


Basically, a flash lamp consists of two elec- 
trodes sealed into a glass envelope in a rare gas 
atmosphere at a pressure of from 10 to 60 cm 
of mercury. The flash occurs when a condenser, 
charged to a potential difference of several 
thousand volts, is discharged through the lamp, 
the flash being initiated by a suitable triggering 
arrangement. 

For satisfactory ballistic photography the 
duration of the flash, that is, the time during 
which the illumination of the projectile surface 
is sufficient to produce noticeable blackening on 
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Fic. 1. The light source, approximately to scale. 


the photographic film, must be very short in 
order to prevent blurring of the image. The 
permissible time is of the order of one or two 
microseconds for projectile velocities of 3000 
ft./sec. and faster. 

The light pulse produced by a flash lamp rises 
rapidly to a peak of very high intensity, often in 
a fraction of a microsecond, and then decays 
more slowly. The duration of the flash therefore 
depends principally on the decay time. Some of 
the lamp designs tested, using a _ photo-cell 
receiver coupled to an oscilloscope, showed a 
plateau in the decay curve which prolonged the 
duration beyond the limits of tolerance. Others 
exhibited oscillations, the output consisting of 
four or more distinct flashes. By arranging the 
geometry of the lamp so that the discharge is 








confined to a relatively narrow space between 
two glass surfaces, it is found that these ob- 
jectionable phenomena may be satisfactorily 
reduced. The duration is then fixed by the 
capacitance of the condenser used to fire the 
lamp! and by the time constant of the discharge 
circuit. By using a small enough capacitance the 
quenching becomes sufficiently rapid to yield a 
duration time of the desired magnitude. The 
addition of some hydrogen to the gas content 
has been found to aid in quenching the discharge. 
The energy input to the lamp, upon which the 
peak intensity depends, is also decreased as the 
capacitance is reduced, but, if the voltage to 
which the condenser is charged is increased, the 
intensity may be maintained at the desired 
level. This voltage, of course, must not be high 
enough to cause the lamp to fire without being 
triggered. The allowable voltage depends on the 
kind of gas used in the lamp and upon the 
pressure. The addition of some hydrogen, men- 
tioned previously as an aid in quenching the 
discharge, also raises the breakdown voltage. 
The lamp design finally adopted} may be seen 
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Fic. 2. The electric circuit 
(firing circuit in upper part 
of diagram; triggering cir- 
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'P. M. Murphy and H. E. Edgerton, J. App. Phys. 12, 848 (1941). 


he lamp is now 


turer’s numbers R-4210 and R-4211). 
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being manufactured by Sylvania Electric Products, Inc., Salem, Massachusetts (manufac- 
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in Fig. 1. The Pyrex rod is supported by the 
electrodes, each of which is 3 inch in diameter 
and 3 inch long. Made of steel, they are drilled 
at one end to permit easy insertion of the Pyrex 
rod and at the other end to receive a 100-mil 
tungsten support rod. The Pyrex envelope is 
blown out at each end to form a bulb. It should 
be noted in the figure that the ‘‘flare out’’ for 
each bulb begins before the electrode is reached 
so that the 1-mm spacing is preserved between 
the edge of the electrode and adjacent glass wall. 
This construction permits the gas, heated by the 
discharge, to expand into the bulbs. The tungsten 
support rods are beaded with Nonex and sealed 
into the envelope with uranium glass. The 
dimensions of the various parts of the lamp are 
shown in the figure. Pyrex tubing with a 2-mm 
wall has been used for the envelope, but a 
somewhat thicker wall is permissible and would 
probably extend the life of the lamp. 

The lamps are evacuated, baked at 360° C for 
one-half hour in a furnace, and filled to a pressure 
of 45 to 60 cm of mercury with a mixture of 
hydrogen and commercial krypton (90 percent 
Kr and 10 percent Xe). The amount of hydrogen 
(approximately 25 percent of the total pressure) 
and the total pressure are adjusted by trial so 
that, unless triggered, the lamp will not break 
down when a potential difference of about 10,000 
volts is applied. The life of the lamps when 
operated at such high energy inputs varies con- 
siderably from one specimen to another and is 
generally terminated by mechanical failure. 

An external electrode for triggering (not 
shown) is finally added. It consists of one turn 
of annealed brass wire around the middle of the 
lamp. 


THE ELECTRIC CIRCUIT 


The electric circuit, shown in Fig. 2, was 
developed from that of Edgerton, Germeshausen, 
and Grier.? The triggering circuit was designed 
for a microphone pick-up, to be connected at B, 
with one stage of amplification. The amplified 
signal from the pick-up device causes an FG-98A 
thyratron to become conducting. The condenser 
in the plate circuit then discharges through the 
thyratron and the primary of an automobile-type 


2H. E. Edgerton, K. J. Germeshausen, and H. E. Grier, 
J. App. Phys. 8, 2 (1937). 
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Fic. 3. Photograph of the complete unit. 


spark coil and the high voltage pulse from the 
coil secondary triggers the lamp. The thyratron 
bias voltage may be adjusted as indicated for 
the desired triggering sensitivity. Provision is 
made for by-passing the amplifier (terminal C) 
when a high level signal from some other type of 
pick-up is available. 

The lamp is fired by a capacitance of 0.5 mf 
charged to somewhat more than 10,000 volts. 
The power supply utilizes two 1-mf condensers 
in series in a voltage doubler circuit, these con- 
densers also being used to fire the lamp. The 
leads from the condensers to the lamp must be 
very short for satisfactory performance. The 
weight of the power supply has been reduced by 
the use of a small 4000-volt neon sign trans- 
former. A potential divider across one of the 
firing condensers furnishes a pulse for operating a 
chronograph (terminal D) when the !amp is fired. 

The addition of a small Variac (not shown) to 
the lamp power supply circuit is advisable in 
order to permit operation of the lamp at some- 
what reduced voltage when maximum intensity 
is not required. 

The circuit components are mounted on a 
metal chassis. Brass rods, fastened to the chassis, 
serve to hold a 14-inch metal reflector in which 
the lamp is mounted with shock resistant sup- 
ports. The electrical connections to the lamp are 
passed through holes in the reflector in poly- 
styrene tubes. A sheet metal case with a plate 
glass window completes the assembly (see Fig. 3). 
The total weight of the unit is about 40 lbs. 


MEASUREMENT OF DURATION 


The duration of the flash must, of course, be 
defined in terms of the type of light sensitive 
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Fic. 4. Photograph of a 90-mm shell in flight. 
Velocity 3000 ft./sec. 


receiver employed. Since in ballistic photography 
an image of the projectile is formed on a photo- 
graphic film, the following method of measuring 
duration time was devised. 

A narrow slit, illuminated by radiation from 
the flash lamp, was imaged on a distant screen 
by a lens. A rotating mirror, driven at 3600 
r.p.m., was inserted into the light path near the 
lens in order to sweep the image of the slit across 
the screen. The lamp was triggered at the proper 
instant to cause the image to pass over a strip of 





Fic. 5. Photograph of an 8-inch, armor piercing shell 
in flight. Velocity 2800 ft./sec. 
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Verichrome film tacked to the screen. The length 
of the blackened area on the film in the direction 
of motion of the light beam was then measured, 
and the duration of the flash determined. Values 
obtained in this way from several lamps ranged 
from one to two microseconds. 

Another type of receiver, consisting of a 
caesium-surface, vacuum photo-cell coupled to 
the vertical deflection plates of a cathode-ray 
tube through a wide-band amplifier, was also 
available. Traces taken with this equipment 
indicate a duration of about four microseconds 
measured from the beginning of the rise to a 
point where the intensity has fallen to about 
one-tenth of the peak value. 

The difference between these values of the 
duration time may be due to failure of the film 
to respond to the lower levels of intensity at the 
beginning and at the end of the flash. 


MEASUREMENT OF PEAK INTENSITY 


The peak intensity of the flash has been 
measured with the photo-cell receiver mentioned 
in the preceding section. A 55-candlepower 
tungsten lamp operated at a color temperature 
of 2848 K was used for calibrating the receiver. 
The peak intensities of various lamps fired from 
a 4-mf condenser at 10,000 volts were found to 
lie between 4X 10° and 6X 10° equivalent candle- 
power with respect to a caesium detector. For 
these measurements the lamps were not mounted 
in a reflector. If a photo-tube possessing the 
approximate spectral response for some common 
photographic emulsion together with other essen- 
tial characteristics had been available it might 
have yielded an intensity rating more appro- 
priate for photographic purposes. 


PERFORMANCE 


As an indication of the performance to be 
expected, two photographs, taken at the Aber- 
deen Proving Ground, are reproduced. Figure 4 
shows a 90-mm shell at a velocity of 3000 ft./sec. 
Electrostatic triggering was employed. An 8-inch 
shell at 2800 ft./sec. is shown in Fig. 5. The 
shell was painted white in order to improve the 
reflectivity, but only the “windshield” retained 
the paint during firing. Magnetic triggering was 
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employed, part of the solenoid in which the 
pulse is induced appearing in the photograph. 

Both photographs were taken with a 4X5 
Eastman Graphic View Camera using a 127-mm, 
f:4.7, lens and Ansco Triple S Ortho film. 

The use of several of the flash units suitably 
placed along the trajectory has yielded quanti- 
tative information on spin velocity and decelera- 
tion. 
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Erratum: Remarks on Some Recently Developed Devices for Summing 
Fourier Series for Crystal Structure Analysis 
[J. App. Phys. 18, 1133 (1947)] 


C. L. Curist 
Stamford Research Laboratories, American Cyanamid Company, Stamford, Connecticut 


DOUBLE summation sign was inadvertently omitted in Eq. (1); this 


should read 


p(x, 9) =L L | Fuo| +E | Fino| cosl2a(hx-+ky) — anno}. 
h_ k h ik 
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On the Investigation of Specimen Contamination in the Electron Microscope* 


James HILLIER 
RCA Laboratories, Princeton, New Jersey 


(Received March 11, 1947) 


A number of experiments directed toward determining the nature and cause of deposits 
which occur on electron microscope specimens under electron bombardment are described. 
From these the conclusion is drawn that the material is a result of chemical reaction, occurring 
at any surface simultaneously bombarded by electrons and organic molecules. The organic 
vapors are shown to arise in two ways: (1) by diffusion out of the metal walls, gaskets, and 
greases exposed to the vacuum, (2) by the outgassing of metal parts of the instrument exposed 
to electron bombardment. To eliminate the effect both sources must be removed. 


I, INTRODUCTION 


LECTRON. microscopists have been aware 
that specimens examined in the electron 
microscope undergo changes in character upon 
bombardment by electrons. A few of these 
changes in character are due to effects which are 
characteristic of the specimen or some part of it. 
One phenomenon, however, appears for all speci- 
mens and seems to be more or less independent 
of them. It consists of a deposition upon the 
specimen under bombardment of a light amor- 
phous material. This particular effect has be- 
come known among electron microscopists as 
specimen ‘“‘contamination.” The observed effects 
of this contamination are somewhat different, 
depending on the nature of the specimen and the 
method by which it has been mounted. In the 
case of filmless mounting or of particles mounted 
on a membrane, the effect appears as a continual 
increase in particle size. In the case of very small 
particles, either on or imbedded in a plastic mem- 
brane, the effect appears simply as an increase 
of thickness of the whole membrane accompanied 
by a loss of contrast. The effect has become in- 
creasingly important as the techniques of appli- 
cation and the resolving power of the electron 
microscope have been increased until, at present, 
in this laboratory, it has become one of the major 
technical difficulties. In the presence of this effect 
it is extremely difficult to obtain accurate meas- 
urements of particle size of low density materials 


*The main part of this paper was presented at the 
Pittsburgh 1946 meeting of the Electron Microscope So- 
ciety of America as discussion to a paper by J. H. L. 
Watson, “An effect of electron bombardment on carbon 
black.” (J. H. L. Watson, J. App. Phys. 18, 153 (1947).) 
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in the size range below 100A. In the case of film 
specimens which contain only small density 
differences, it quickly reduces the contrast to the 
vanishing point. This effect has also precluded 
any possibility of preserving interesting fields of 
view for subsequent examination. 

In a recent paper, Watson! has described some 
observations made on the changes resulting from 
this effect in the particle size of various carbon 
black preparations. He offers the hypothesis that 
the deposit is the result of a polymerization and 
condensation of organic vapors and gases under 
electron bombardment. 

In this laboratory, contamination phenomena 
have also been observed, and intensive studies 
have been made over a long period for a great 
variety of specimens and operating conditions. 
While the results obtained cannot be considered 
to be completely conclusive, they do throw con- 
siderably more light on the mechanism of this 
phenomena and give an indication as to a method 
by which it can be minimized. 


Il. EXPERIMENTAL RESULTS 


The experimental results and conclusions pre- 
sented below have been collected over several 
years time and are the results of the efforts of a 
number of workers and from work on a number 
of different instruments in addition to the elec- 
tron microscope. 

(1) The first appearing deposits observed by 
ordinary light are brownish and translucent. 
Heavy deposits obtained after long electron 
bombardment are black and opaque. 

(2) The deposits are insoluble in any of the 
common solvents and are not affected by strong 
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Fic. 1. High magnification electron micrograph of 
carbon black particles on which a layer of contamination 
has been deposited. The image was taken on the long focal 
length side of focus to enhance contrast. The original 
carbon black particles can be just distinguished as slightly 
denser regions inside the mass of contamination. 


acids or alkalies. Deposits on metal surfaces dis- 
appear if heated to red heat in air. Qualitatively, 
their properties appear to be independent of the 
substrate upon which they are deposited. 

(3) They have an electron scattering power 





(a) 


slightly less than that of carbon (Fig. 1). They 
are structureless to the electron microscope and 
amorphous to electron diffraction. 

(4) The deposits formed in the beam of the 
electron microanalyzer' on a silica substrate? 
have been analyzed by that instrument. Carbon 
only was found, but, because of the inability of 
the instrument to detect hydrogen, the observa- 
tion does not exclude the possibility that the 
material contains hydrogen. Similarly, there may 
be other elements present in low concentration. 

(5) The material is not stable in thick layers 
in air. When a membrane and screen is bom- 
barded for several hours so that an excessively 
thick deposit is produced (Fig. 2a) it is often 
found that on standing overnight in a partial 
vacuum of the instrument the deposit loses its 
homogeneity and also shows signs of having 
migrated over considerable distances (Fig. 2b). 

(6) The rate of formation is roughly propor- 
tional to the intensity of bombardment. 

(7) It has occurred on all materials examined, 
and, while it has not been proven that the rate 
of formation is the same for all substrates, quali- 





(b) 


Fic. 2. (a) Electron micrograph of a hole in a collodion 
membrane. The membrane has been greatly thickened by 
contamination. (b) A similar hole in the same membrane 
after it had stood overnight in a partial vacuum of the 
electron microscope. The migrating contamination has 
bridged the hole with a thin film and has collected into a 
granulated layer in the central region. 


! J. Hillier and R. F. Baker, “‘Microanalysis by means of electrons,” J. App. Phys. 15, 663 (1944). 
2 J. Hillier and R. F. Baker, ‘Further developments on the microanalyzer,” J. App. Phys. 16, 266 (1945). 
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Fic. 3. Cross section of the arrangement of a liquid air 
chamber enclosing the specimen in the microanalyzer. In 
this arrangement the electron beam travels upwards on 
the axis of the system. 


tative observations indicate this to be the case if 
all other experimental conditions are carefully 
controlled. Actually, it would be expected that 
the initiation of the deposit might vary for differ- 
ent substrates, but, obviously, once the deposit 
is initiated the substrate is the material formed. 
Thus, once the deposition has been started, it 
would be expected that the addition of further 
material would proceed at a rate determined only 
by the experimental conditions in the instrument. 
(8) Surrounding the specimen (microanalyzer) 
and the beam near the specimen in a liquid air 
trap (Fig. 3) produced the following results: 


a. Introduction of the copper chamber (it had been 
freshly cleaned) around the specimen at room temperature 
quickly reduced the rate of the deposit to approximately 
one-fifth of its original value. 

b. Cooling the chamber to liquid air temperature made 
no further change. 


These observations are consistent with those 
made on instruments including liquid air traps 
in the pumping line, in which case no observable 
change in the rate of contamination of the speci- 
men has been observed. This was even true in a 
particular case where the system was never al- 
lowed to be opened to the oil diffusion pump 
without the liquid air trap in operation. It is 
believed in this laboratory that oil vapors from 
the diffusion pumps are not responsible for con- 
tamination. At first sight, it appears that the 
experiment eliminated the possibility of the phe- 
nomena being due to the polymerization of 
hydrocarbons, since all of these would be con- 


228 ° 





densed out at liquid air temperatures. However, 
the arrangement of the experiment was such that 
the rod bearing the specimen was not cooled 
more than a few degrees below room tempera- 
ture, and hence the rate of diffusion of gases or 
vapors from the metal parts would not be affected 
greatly. 

9. When three membranes separated by 50yz 
were introduced in the microanalyzer, it was 
found that the rate of deposition on the central 
membrane when the probe was focused on it 
would decrease rapidly with time but would re- 
turn to its original value when the beam was 
turned off for a few minutes. The results of this 
experiment, as well as those described in 8a, can 
be explained on the assumption that the material 
from which the contamination comes is in the 
vapor phase and that its partial pressure is very 
low. The action of the electron beam appears to 
reduce this partial pressure rapidly within the 
restricted volume near the specimen. In both ex- 
periments the arrangement is such that it takes 
considerable time for equilibrium to be re- 
established after the electron beam is turned off. 

10. Changing the potential of the specimen 
relative to its surroundings by as much as a 
thousand volts positively or negatively made no 
change in the rate of deposition. This observation 
eliminates the possibility that the deposition is 
the result of electrostatic attraction of gas or 
vapor molecules which have been ionized by the 
electron beam. It is known that in the electron 
microscope, where the bombardment usually in- 
cludes part of a wire screen connected to ground, 
the area of the thin membrane under bombard- 
ment is not chargedThis is not true if that area 
includes extremely large particles (greater than 
10 microns) or.if the electron beam is so concen- 
trated that none of it strikes the supporting wire 
screen. The latter is invariably true in the case 
of the probe diffraction camera’ or microanalyzer. 
Assuming the most favorable case, i.e., when 
equal numbers of negative and positive ions are 
involved, at least a 50 percent reduction in the 
rate of deposition should be observed upon ap- 
plication of the potential. Since no such change 
was observed, it can only be assumed that the 


’J. Hiller and R. F. Baker, “On the improvement of 
resolution in electron diffraction cameras,” J. App. Phys. 
17, 12 (1946). 
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electrostatic attraction of the ions to the bom- 
barded area is not responsible for the observed 
deposits. It should be noted that 1000 volts is 
sufficient to repel ions which have been given 
maximum possible kinetic energy by the electron 
beam. 

11. The introduction of magnetic or electro- 
static fields for the deflection of beams (the 
microanalyzer probe was used), immediately be- 
fore it strikes the specimen, have shown that the 
material is deposited only in the area where the 
deflected beam irradiates the specimen. This ex- 
periment was carried out in order to prove that 
the deposit was not due to an ion beam arising 
in the gun of the instrument. It also strengthens 
the suggestion that the surface bombarded by 
the electron beam takes part in the chemical 
reaction resulting in the deposit. 

12. It is observed that the rate of deposition 
decreases with time when the total current of the 
beam is sufficient and if the beam does not strike 
the metal supporting screen. In the electron 
microscope an illuminating current of 3 to 5 
microamperes, focused to a 10-micron diameter 
spot and allowed to run for 20 minutes, reduces 
the rate of deposition to the vanishing point. 
Once the rate of deposition has been reduced to a 
negligible value under the above conditions, it 
does not return to its original value for many 
hours, even if the instrument is opened to the air 
several times in that period. 

13. If the conditions described in (12) are 
achieved and if part of the beam is allowed to 
impinge on the copper specimen screen, the rate 
of deposition of the material on the irradiated 
part of the membrane returns to the higher value 
for as long as the copper screen receives some 
bombardment. In this connection, another phe- 
nomena was often observed, namely, that in the 
case of a film mount a narrow clear strip on 
which there appeared to be no deposits and 
which followed the contour of the wire could be 
seen (Fig. 4). 

Ill. CONCLUSIONS 

All of the evidence given in the above observa- 

tions and discussions supports the hypothesis 


put forward by Watson‘ that the chemical proc- 
ess involved is the polymerization of a hydro- 


4]. H. L. Watson, J. App. Phys. 18, 153 (1947). 
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Fic. 4. A schematic representation of the factors giving 
rise to the clear uncontaminated border that is often ob- 
served near the edge of a grid wire. It is presumed that 
contaminating vapors are released from the grid wire 
(screen support) by the electron bombardment but that the 
released molecules travel in straight lines and are deposited 
only when they strike another area. The region of the 
membrane very close to the grid wire is sheltered from the 
bombardment of gas molecules, and hence no contamin- 


ation appears even though that region is under electron 
bombardment. 


carbon. It indicates further, however, that the 
deposit requires all three factors—the material in 
the gas phase, the surface, and the electron bom- 
bardment—leading to the suggestion that the 
surface under bombardment is activated so as to 
act as a catalyst for the reaction. While the ob- 
servations do not exclude the possibility that the 
material is formed in the gas phase, they do show 
that if this is so the material produced is not 
responsible for the contamination observed. 

The observations indicate that the gas or 
vapors involved diffuse into the system very 
slowly and at no time reach a high partial pres- 
sure. The probable origins are the metal walls of 
the instrument, the rubber gaskets, and the 
sealing or lubricating greases. There is a definite 
indication that it is released more rapidly by 
electron bombardment—probably through local 
heating. This is particularly true in the case of 
the specimen screen. The observation, 13, pro- 
vides the strongest evidence that the deposit 
occurs mainly on the area of the substrate simul- 
taneously bombarded by the gas molecules and 
the electron beam. 

There is good experimental evidence that the 
gases or vapors do not come from the oil diffusion 
pump. The same contamination problem has 
been encountered in instruments which have 
never used anything except mercury diffusion 
pumps with liquid air cold traps. Furthermore, 
it was found that when Octoil vapor was bom- 
barded in an experimental arrangement, the 
solid deposit obtained was white and consisted of 
phthalic anhydride. 


229 











While the evidence is not conclusive it does 
not support the suggestion put forward by Wat- 
son that there is a dependence on the conduc- 
tivity of the specimen. 


In a recent note Cosslett® suggests that the 
contamination is due to the ejection of absorbed 
matter and perhaps of metal itself from the grid 
wires under electron bombardment. He excludes 
organic impurities from the “absorbed matter” 
since cleaning and heating in vacuum prior to 
use brought about little improvement. The pre- 
sented observations, namely, that the contami- 
nation is largely carbon, has lower electron scat- 
tering than carbon, and is structure free to the 
limit of resolution of the microscope in_ this 
laboratory (10A), are almost in direct contradic- 
tion to his statements. The explanation for the 
discrepancy probably lies in the fact that the 
contamination has two sources (the object cham- 
ber space and the grid wires), a possibility which 
Cosslett does not seem to have considered. Since 
these can vary greatly in their relative strengths, 
observations in which the results are not sepa- 
rated have little meaning. The use of the fine 
probe in the microanalyzer system has made this 
separation possible in many of 
experiments. 


the above 


IV. ELIMINATION OF THE EFFECT IN THE 
ELECTRON MICROSCOPE 


These observations and conclusions point to a 
possible method of eliminating this very trouble- 
some phenomenon. 

In the probe type of instrument where the 
extremely fine beam seldom strikes the wires of 
the specimen screen, the rate of deposit can be 
reduced by several orders of magnitude® by the 

5V. E. Cosslett, J. App. Phys. 18, 844 (1947). 

6 N. Davidson and J. Hillier, “‘Single crystal electron dif- 
fraction by micro-crystalline materials,” J. App. Phys. 
18, 499 (1947). 
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simple introduction of an electron source in the 
specimen chamber. Spraying the chamber with 
an electron current of 50 ma at an accelerating 
potential of 1000 v for.1 minute cleans the system 
for several hours operation. Apparently it cleans 
the chamber space and also releases and re- 
deposits the gas near the surface of the walls of 
the chamber. In the microscope a similar scheme 
can be used. In this case, however, the deposit 
occurring as a result of the bombardment of the 
specimen still remains and must be 
avoided if the difficulty is to be completely elimi- 
nated. While this has not been investigated in 
this laboratory, it appears evident that the use 
of screen material from which hydrocarbons can 
be removed, probably by heating in a furnace, 
should provide an answer. The author has re- 
ceived reports that the deposit is appreciably less 
in the case of stainless steel screens that have 
been cleaned by heating. It would appear logical 
to investigate the use of specimen screen ma- 
terials such as platinum or tantalum for critical 
electron microscope work where extensive study 
of a single specimen is desirable. These could be 
cleaned and outgassed by heating in a furnace; a 
vacuum furnace might be preferable in extreme 


cases. 


screen 
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Possible Fluctuations in Electron Streams Due to Ions 


J. R. Prerce 
(Received August 11, 1947) 


The equations for low amplitude space-charge waves in an electron stream containing heavy 
positive ions are considered. At frequencies a little below the ion plasma frequency there 
should be a wave which increases in amplitude as it travels (negative attenuation). Ion motion 
tends to limit the rate of increase. In a beam filling a conducting cylinder and confined by a 
uniform magnetic field there should be an increasing wave at frequencies near the ion cyclotron 
frequency. Such waves might explain observed fluctuations in long electron beams. 





1. INTRODUCTION 


N the course of work at the Bell Telephone 
Laboratories on traveling wave tubes, Dr. 
L. M. Field and the writer observed fluctuations 
in electron streams in high vacua. Subsequently, 
Mr. J. A. Morton and the writer have found 
similar fluctuations in electron streams under a 
variety of conditions. The fluctuations can be 
observed as a fluctuating division of current be- 
tween electrodes due to fluctuations in the posi- 
tion of electron flow. 

In electron flow in a high vacuum, the elec- 
tronic space charge tends to form a potential 
minimum and this provides a mechanism by 
which positive ions can be trapped in the electron 
stream. The presence and motions of ions in the 
electron stream might afford one source of fluc- 
tuations. The purpose of this memorandum is 
not to present experimental data, but to give an 
analysis of the situation which shows that fluc- 
tuations should build up along an electron stream 
if ions are present. The analysis is carried out 
for a case in which the positive space charge of 
ions is just equal to the negative space charge of 
electrons, so there is no average charge density 
in the electron stream. The analysis is a small- 
signal or perturbation theory. 

First an essentially one-dimensional case is 
treated, in which neither d. c. nor a. c. quantities 
vary normal to the direction of electron motion. 
If the ions have no d. c. velocity, this leads to a 
wave which builds up infinitely rapidly with dis- 
tance at a frequency which may be called the ion 
plasma frequency. Such a wave might result in 
the amplification of small initial disturbances, 
such as those due to shot noise, to give the large 
fluctuations observed in some traveling wave 
tubes. Small ion velocities are shown to give a 
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finite rate of build-up. Collisions might result in 
some sort of ‘‘damping.”’ 

The case of a beam filling a conducting tube 
and an axial magnetic field is also considered. It 
is found that there is an increasing wave in this 
case for typical conditions. 


1.1 One-Dimensional Case, Ions at Rest 


We will consider first an essentially one-dimen- 
sional case, in which there are no electric fields 
in the x or y directions and a. c. quantities are 
functions of z only. We will use m.K.S units 
throughout. We will consider the case of a current 
of electrons of magnitude Jp/cm*, traveling in the 
+z direction with a mean velocity 


e i 
w= (2=10) ° 
m 
Here e/m is a positive quantity, the charge to 


mass ratio of the electron. The mean charge 
density of electrons will be 


(1) 


Peo = — Ip/up= — po- (2) 
We will assume a mean ion charge density 
Pi0 = — Ped = Pos | (3) 
composed of ions of charge to mass ratio 
q/M=(e/m)/N. (4) 


We will assume that the frequency is low 
enough so that a. c. magnetic fields can be 
neglected and the electric fields can be given well 
enough in terms of a scalar electric potential. We 
will consider small a. c. disturbances only, omit- 
ting all but first-order a. c. terms. This will 
linearize the equations and result in solutions in 
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which the a. c. quantities vary as 


exp(—I's+ jot). 


The a. c. quantities will then be represented by 
the constant factors identified below multiplied 
by exp(—I's+ jw). 


V =potential, 

p.-=electron charge density, 

pi =ion charge density, 
v.=electron velocity, 

vi =ion velocity, 

ge= electron convection current, 
qi =ion convection current. 


In terms of these quantities, we may write 
(Poisson’s equation) : 


1 
r?V = ——(pit+p,). (5) 


€ 


Definition of first-order convection current 


Ge = Peo — poe, (6) 
Ji = pov i- (7) 
Conservation of charge 
jwpe=Vq., (8) 
jwp: =Vqi. (9) 


First-order force equation 
jJwv.— Ul v, = — (e/m)TV, (10) 
(11) 


In combining these equations it is convenient 
to define quantities 


jv; = (e/m)TV/N. 


we = (po(e/m) /e)', 


wi = (po(g/M)/e)'=w,/N'. 


(12) 
(13) 


It is interesting to have some idea of the magni- 
tude of w, and w;. Suppose Jy is the current 
density in amperes/cm?. Then 


we = 1.83 XK 10'T,3/ Vo!. (12a) 


For the reasonable conditions of 0.1 ampere/cm? 
and 1500 volts 


w@, = 9.3 X 108 = (27) 150 10° 
for 
0-, N=2.9X10', N'=170 
w;=5.5X 10°= (27).87 K 10°. 
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Combining (5)—(13), we obtain 


we" wi" 
rv=( —_—___—+ rev. 
(w+j0'uo)? w? 


We might expect three types of solution, l=0, 
V=0, and 


(14) 


w," w;” 
{ =——____+—. (15) 
(w+ ju)? w? 
Those represented by (15) are the ones of inter- 
est. We can solve (15) for T° 


o wl fw" — 
-r=-j-4"|() -1| : 
Ug Ug (4) 
If w<w; this represents two waves traveling to 
the right with a phase velocity equal to that of 
the electron stream. One wave has a positive 
attenuation and the other a negative attenuation. 
We may imagine, for instance, that small ir- 
regularities in the electron stream at a point (say, 
a starting point) can be represented as having a 
component of the increasing wave. Then such 
irregularities can be the basis of large irregulari- 
ties further to the right, as the wave increases 
according to (16). The wave increases at an in- 
finite rate for w=w;, and hence we would expect 
very large disturbances near this frequency. 


(16) 


1.2 Effect of Small Ion Velocities 


The ions are presumably generated by electron 
collision, and we would expect a small ion drift in 
the direction of electron motion and a distribu- 
tion of ion velocities. 

Suppose we say that there is an ion charge 
density dpp associated with ions having z com- 
ponents of d. c. velocity lying in the range du 
at u. We will now have instead of (7), (9), 
and (11) 





dq: =dput+dpo;, (7a) 
jwdp;=Tdqi, (9a) 

1 e 
jov; = ulTv;=— —T V. (11a) 

m 

From these we obtain 
(e/m) r?V 

dpi= — (17) 


NN Wwo+jru)? 
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And, corresponding to (15) we have 


W," 1 u= dpo 
= tue f ——_———. (18) 
(w+ jlo)? po % ~~» (w+ ju)? 
We can progress a step further by assuming 
that none of the ions has a very large velocity, 
so that always 


Tru<w. (19) 


In this case we can rewrite (18), expanding the 
denominator of the integral 


1 f dpo 
po (w+ jP'u)? 


11 px v rT 
=— — (1-2-3 w--- dog 
w* po Fux w w” 


1 jv r? 
=— 1-2—-a—3—(ut)—---) (20) 
w* w w* 
Here we are describing the ion velocity distribu- 
tion in terms of the moments @, (u?), etc. If we 
use this expression, in connection with (18) we 
obtain 


We” 


1 =——_____ 


Tuo\? 
(sf) 
w® 
1'(u*) 


w,;” Ta 
(1 — j2——3—__-.. 


w~ w @* 


-). (21) 


For each additional moment included in describ- 
ing the ions, the equation becomes one degree 
higher in I’. This is natural, as inclusion of a, 
(u*), etc., requires some means for adjusting the 
ion motions at boundaries where the ions enter. 

The reasonable step is to include only the term 
in a Let us do this, rewriting (21) in terms of 
the parameters W and y 


~ 
w-(=) , (22) 
Wi 
Tuo 
sin (23) 
«; 


We obtain 


w.\? 7\? iy 
(~) = ( +i-) | w:- 1+ j2— =| (24) 
Wi W Uo W. 
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We are interested in frequencies such that W is 
very nearly equal to unity. 


@ 


—>|, (25) 
w; 
u 
—1l. (26) 
Uo 
Hence 
|y|>>1. (27) 
We may thus write approximately 
w\? —y¥y? “iy 
(*)'-2 Yn] ao 
@; Ww? uo W 
For a special case, W=1, this yields 
; we\'/2a\-! 
v-(9(=) (—) (29) 
Wi Uo 
2a w\'/2a\! 
. v= (!(=) (—) . (30) 
Uo w; Uo 


For singly charged ions we have from (4) 
M/m=N. 


If a, a, etc., are produced by collisions with elec- 
trons, we might expect &/up to be of the order of 
1/N, (u*)/u? to be of the order of 1/N?, etc. We 
have from (13) that (w,/w;)?=N. Assuming 


a 
—=1/N, (31) 
uo 

we obtain 

vy =(7)*(2)-3 N53, (32) 
2% . 
—y = (j)8(2)t3N-3, (33) 
Uo 


For 0-, N=2.9X 10‘ we have 
2a 


I—y| =5x10-2. 
|\Uo | 


Thus, (20) and (24) appear to be reasonably well 
justified. We have, then, using (31), 


—¥1=7.50X 104( —.866 — 7.5), 
—y2=7.50X 10*(+.866—j.5), 
—73=7.50X 104(j). 


71 represents a decreasing wave traveling to the 
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right, y2 an increasing wave traveling to the right, 
and 3; an unattenuated wave traveling to the left. 


1.3 Damping, as by Collisions 


The kinetic energy of an electron is for a given 
velocity, proportional to its mass. If we assumed 
some sort of frictional force (as, due to collisions) 
the effect could be taken into account in our 
analysis by assuming the mass or the charge-to- 
mass ratio to have a small complex component. 
It would seem offhand that the electrons would 
be more affected by collisions than the ions. 
However, we see that a complex value of e/m, 
meaning a complex value of w,, would not prevent 
lr from going to infinity. It requires damping of 
the ions to do this. Suppose, then, that there is a 
damping term, so that 


wi=w (1+ 76), 


Jwi = jw, — Sw. 


(34) 


If <1, this gives approximately 


wo wf /or\? _fay?y? 
roa jta't(*) 14202) 05 
Up Uhl \w w 
For instance, at w=, this gives 
: ®@ 1 We 1 We 
er 
uo 265) uo 25! uo 
and obviously [' does not go to ~ for any real 
value of w. 


(36) 


1.4 Cylindrical Beam with Magnetic Field 


The question now arises, what is the effect of 
confining the beam in a conducting tube and of 
having a uniform axial magnetic field. Suppose 
we assume an axial field of strength B in the z 
direction. Let us deal initially merely with the 
motions of electrons, dropping the subscript e. 
Let v., vy, v, be the Cartesian components of 
velocity. If we assume a variation in the z direc- 
tion, exp—T,, as before, we obtain, saving first- 
order terms, 


e fav 
(jw — Tuso.=—(—+Br, ), (37) 
m\ dx 
e /aV 
(jo—Tmso,=—(——Bo.), (38) 
m\ dy 
e 
(ja—Tuo)v,= ——T V. (39) 
m 
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We further have from the conservation of charge 


Ov, OV, OV, Op ; 
~ (4+ ) +00? — jwp, 











Ox Oy 02 02 
Ov, OV, 
(jw — Tuo) p= mo +—- rv). (40) 
Ox ody 
If we solve (37) and (38) we obtain 
¢ OV e aV 
—( jw — P'u)— +—ay— 
m Ox m dy 
v.=—_______—_—,,_ (41) 
(jw —T'uo)?+wo? 
OV e aV 
—(jJw — 'uo)— — —wo— 
m oy m Ox 
v,= <i tel (42) 
(jw —T'uo)*+ wo? 
where 
e 
wo=—B. (43) 
m 


Combining (39), (41), and (42) with (40) we 
obtain 


a°V 

e Ox? 

p= —-—Ppo : = - re ye . 
m \(wt+jl'uo)?—wo? (wt jlo)? 


o7V 


ay? mv 





(44) 


It has been convenient to work the problem 
out in rectangular coordinates. However, we wish 
to apply the result to an axially symmetrical 
system with a conducting wall. Suppose we 
choose a solution in which V varies as 


Jo(krje-™?, 
where 


ka=2.4. (45) 


The radius of the conducting tube being a, then 














eV aV 
——+—= —k?\. (46) 
Ox? dy? 
Poisson’s equation will be 
eV aV &V Pet pi 
pen nnn E aiemenmannn, 
Ox? dy? dz? € 
(47) 
Pet pi 
(—k?+Tr?) V= ———_. 
€ 
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It is easy to see that using (12) and (13) and 
letting 


e 
—B=wr, (48) 
m 
g. 
—B=woi=woi/N, (49) 


M 


we obtain, disregarding the solution V=0 (for 
which v, and v, need not be zero), 





— p2 2 
—k? +1? =] —— +——_] 
(w+ jl'uo)?—woe? (w+jl'uo)? 
— p2 rr 
+e] ———+—] (50) 
w?—wo;? w? 


If there is no magnetic field (we. =wo;=0), this 
reduces to (15) and 

—#+I?=0. (51) 

Thus, we have the waves given by (16) and also 

r= +k. (52) 


The waves given by (52), however, have a poten- 
tial identical with that in the tube in the absence 
of charge and are essentially ‘‘field waves’ as 
opposed to ‘‘space charge” waves. 

It is important to know the values of wo. and 


woi. If H is the magnetic field in gauss, 
woe = 17.8 10°H. (53) 


A focusing field of 100 gauss is common in a 
traveling wave tube. In this case 


Woe = 1780 X 10° = (27) 280 X 10°. 
For 0-, N=2.9X 10 and for this case 
wo; = 61,000 = (27) 10,000. 


(54) 


As (50) is rather complicated, it is helpful to 
use dimensionless units. 

We have earlier found important effects in the 
frequency range around w;. Let 


—=W, (55) 
Wi 
Tuo 
——mey, (56) 
W; 
kuo 
— =m. (57) 
WwW; 
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We then obtain 





w.\? 1 
Pais 
W; (W+ jy)? 


m? /* 


Woe 2 
(W+ jr)? (“) 
Wi 


It is important to know the magnitude of the 
factors m? and (w,/w;)?. For 0- as compared with 


electrons 
we\? 
(=) = N=2.9K10+. 
Wj 


Let us assume that r=0.24 cm. Then 





(58) 


(59) 


kr = .24k=2.4, 
k=10. 


Also, let us assume, up =c/10=3X 10° cm/sec. If 
we take the value of w; we assumed previously, 
w:;=5.5X10°, we find from (57) that 


m* = 3.0 X 107. 


The left-hand member of (58) is very small, while 
Y appears in the denominator of all the terms 
on the right. Thus, one solution is one in which 
y is very large. If this is so, (58) becomes 
approximately 


{ { 1 ) we\2)3 
y=u m? | 1 —————_-_——_ | -—- (~) 
oi W; 


ba m-(2¥ | 


: | ; (60) 


| W?} 








This can give either imaginary values of y, or 
real values representing either increasing or de- 
creasing waves. y is infinite and changes from 
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real to imaginary at 


W=1, 


WwW >= Wj. 


Thus, near this value of W our assumption that 
y is very large is justified. This is the ion plasma 
frequency which we previously estimated as 0.88 
megacycles. Also, y is infinite and changes from 
real to imaginary at 


w Wi 
=—=>—-, 
@i @; 

w@ >= Woi- 


This is the ion cyclotron frequency, which for 
100 gauss is 10 kilocycles. 

We see that increasing waves (—y real and 
positive) are to be expected in a beam surrounded 
by a conducting tube and confined by a magnetic 
field, and one of these is of a frequency near the 
ion plasma frequency while the other is near the 
ion cyclotron frequency. These are, of course, not 
the only values of y, for (58) is of the sixth 
degree in y. The other roots of the equation give 
unattenuated waves. 


1.5 Summary of Theory 


If we set up the linearized equations for an 
electron stream containing positive ions, we find 
that among the solutions there are waves which 
increase linearly with distance in the direction of 
electron flow. In both a one-dimensional case, or 
in the case of a beam surrounded by a conducting 
tube and confined by a magnetic field of any 
strength, such increasing waves occur at fre- 
quencies of the order of the ion plasma frequency, 
given by (13) and the ion cyclotron frequency 
given by (49). For physically reasonable condi- 
tions the first may be of the order of one mega- 
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cycle while the second may be of the order of 10 
kilocycles. 

When the boundary conditions are imposed at 
the end of the tube at which electrons enter, it 
must, in general, be found that the small fluctua- 
tions present give rise to an* increasing wave 
which will become much larger further along the 
stream. Thus, the existence of waves which grow 
in the direction of electron motion would seem 
to mean that disturbances at frequencies near 
the ion plasma frequency and the ion cyclotron 
frequency should always appear if the beam 
travels far enough. It is possible also that second- 
ary electrons traveling backward along the beam 
might provide feedback and result in oscillation 
as opposed to mere amplification of disturbances. 

If there is no damping and the ions have no 
d. c. velocities, the rate of increase of the waves 
becomes infinite at a certain frequency. It is 
shown that if the ions have a d. c. velocity the 
rate of increase does not become infinite. Pre- 
sumably, damping of the ion motion, as, by col- 
lisions, would make the ion plasma frequency 
complex, representing a decaying oscillation, and 
would also preclude an infinite rate of increase 
of amplitude with distance. 

Finally, non-linearities would soon limit the 
rate of increase along the tube. 


2. CONCLUSIONS 


Theory predicts that disturbances in an elec- 
tron stream containing ions will build up in the 
direction of electron motion. Experiment shows 
the existence of oscillations of about the right 
frequency and having roughly the right proper- 
ties to correspond with theory. Secondary elec- 
trons could carry disturbances back to the 
cathode end of the tube and result in oscillation 
rather than mere amplification of noise. 
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Application of the Laplace Transform in the Solution of Linear Integral Equations* 


LAWRENCE BayLor RosINson! 
Department of Physics, Howard University, Washington 1, D. C. 


(Received August 11, 1947) 


Integral equations are defined and references are given to the Laplace transform technique. 
Most of the operations used are illustrated through a solution of Abel’s integral equation. It is 
shown that the Liouville-Neumann series follows logically from an operational solution of 
Volterra’s integral of the second kind. Whittaker obtained solutions to the integral equation 
with the kernel expressed in the last three forms. His procedure was quite involved and it is 
shown here that the same solutions are obtained more directly through the operational calculus 


based on the Laplace transform. 





INTRODUCTION 


N integral equation? is one in which the 
unknown function appears under the sign of 
integration. The integral equations of greatest 
importance to the physicists are those known as 
Abel’s equations, 


t 
F=f (t—0)-*Y’(@)d@,, (0<b<1), (1) 
0 
and Volterra’s equation of the second kind, 


P(t) = Ft f (0) K (t—@)dé. (2) 


0 


These equations have been solved by various 
means.* The suitable solutions for computation 
purposes are those given by Whittaker. The 
procedures whereby the solutions are obtained 
are quite involved and a special method of attack 
is given for each type of problem. A general 
method of solving equations such as (1) and (2) 
is given by the operational calculus, based on the 
Laplace transform.* Koizumi used this method to 





* This article is taken from some of the lectures given by 
the author in a course in Mathematical Methods in Physics 
at Howard University. 

' Present address: The Physical Sciences, The College 
University of Chicago, Chicago, Illinois. 

2M. Bécher, An Introduction to the Study of Integral 
Equations (Cambridge University Press, Teddington, Eng- 
land, 1909); W. V. Lovitt, Linear ‘Integral Equations 
(McGraw- Hill Book Company, Inc., New York, 1924). 

3H. Margenau and G. M. Murphy, The Mathematics of 
Physics and Chemistry (D. Van Nostrand ny Inc., 
1943), pp. 503-525; E. T. Whittaker, Proc. . Soc. 
oan (A)94, 367 (1917); S. Koizumi, Phil. os 2, 432 

1931) 

4G. Doetsch, Theorie und Anwendung der Laplace- 
Transformation (Dover Publications, New York, 1943), pp. 
155-167. R. V. Churchill, Modern Operational Mathematics 
in Engineering (McGraw-Hill Book Company, Inc., New 
York 1944), pp. 36-44. 
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solve some integral equations and H. T. Davis® 
treated integral equations from other operational 
points of view. It is the purpose of this article to 
show that the use of the Laplace transform 
provides a more general method for solving some 
types of integral equations than is generally 
utilized, and that the solutions obtained by 
Whittaker are obtained by a ‘‘turn the crank 
procedure” through the operational calculus. 
This writer has employed the notations and 
methods of the operational calculus given by G. 
Doetsch in his Theorie und Anwendung der 
Laplace Transformation. For the reader unfamiliar 
with these methods, the clearly written and 
concise articles of L. A. Pipes® will provide a 
satisfactory background. Pipes defines the trans- 
form in a slightly, though not significantly, 
different manner from Doetsch, however. The 
operational technique has been used in most cases 
where the kern (also known as kernel and nucleus) 
of the equation K(¢—@) is given by a simple ex- 
pression in closed forms. Herein the method is 
discussed when the kern is expressed in other 
ways. 
ABEL’S INTEGRAL EQUATION 


Equation (1), Abel’s integral equation, has 
been solved by Whittaker and also by Koizumi, 
Doetsch, and Churchill using the operational 
calculus. A brief account of the solution is given 
here largely for purposes of defining the opera- 
tions contained herein. 


r= f (t—0)-°Y'(6)d@, (O<b<1), (1) 


8H. T. Davis, The Theory of Linear Operators (The 
Principia Press, Bloomington, Indiana, 1936), pp. 468-514. 
*L. A. Pipes, J. App. hys. 10, 172, 260, 30r (1939). 
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is to be solved for Y(t) where Y’(t) is dY(t)/dt. 
This equation may also be written as 


F(t) ={-* y(t), (3) 
where F,(t)*F2(t) is defined as 


F,(t)* F2(t) = f F,(t—6) F2(0)dé. (4) 


It is also shown that F,(t)*F2(t)=F.(t)*Fy(¢). 
Equation (4) defines the faltung of the functions 
F(t) and F,(t) ; this is also called the convolution 
integral. The Laplace transform of a function 
F(t) is by definition 


L{ F(t)} =s(5)= f e-*'F(t)dt, (5) 


0 


where the real part of s is greater than zero, or 
R(s)>0. F(t) must be sectionally or piece-wise 
continuous (stickweise stetig) and of exponential 
order so that this integral will converge as [> =~. 
The transform of the convolution is the product 
of the two individual transforms or 


L{ F,(t)*F2(t)} =fils)fe(s). (6) 
The transform of a derivative, as in Eq. (1) is 
shown to be 

L\ F'(t) } =sf(s)— F(+9), (7) 


where the second term on the right of (7) is read: 
the value of the function as ¢ approaches zero 
from the right. When Eq. (4) is transformed, the 
following equation results: 


f(s) =P(1—6)[1/s* JLy(s)—(¥(O)/s)]. (8) 


Equation (8) is then solved for y(s) and the 
inverse transform taken, whence the required 
solution is obtained. The inverse transform is 


L-\x(s)} =X (0) 


1 y+ ip 
=—— lim f exp(st)x(s)ds, (9) 
2Ti p+ J,_ig 


where s is a complex variable and the integration 
is performed in the s-plane. From (8) it is 
apparent that y(s) can be written as 


(0) s-* 


Y 
y(s) =——+ 
Ss 





—f(s). (10) 
r(i—d) 
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The inverse transform 








s-> te sinbr 
i = —_= t. (11) 
r(1—d) r(o)r(i—b) rg 
It is evident that 
d sinbr 1 
Y(t) = Y(0) +— —— —*Fi(), (12) 
dt -« tt 
or 
sinbx d é dé 
Y(t) = Y(0)+———- — F(@)———. (13) 
r dtd, (t—0)° 


As was afore-mentioned, the above solution is 
well known and is given here only as the means 
of introducing the definitions of the operations. 


VOLTERRA’S INTEGRAL EQUATION OF THE 
SECOND KIND 


1. The Liouville—Neumann Series 


Margenau and Murphy solve the Volterra 
integral equation of the second kind, 


(1) = F(t) +a (0) K (t—0)d8, (2) 


0 


by the Liouville-Neumann series: 


(?) < A", (t), (14) 


n=0 
where 


$,(t) = F(d), 


6,() = J K(t—6)0(0)d8, 


t 
$,() -{ K (t—0)4,(6)d8, (15) 


e.()= f K(t—6)®,_,(0)d0. 
0 


It is shown that the series (14) converges uni- 
formly for |\| <1/Mt where M is some positive 
number. Equations (15) may be also written as 


w= f KU-9 f° KU-0-~ 


tr—2 
x f K (tp2—te-1) F(ty1)dt,_1++-dt,d0. (16) 
a 
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Consequently, the solution to Eq. (1) may be 
written as 


&(t) = F(t) +0 f K (t—0)(6)d0 


+x f Kio f K (t,—t) F(t)dtd@+ - +++ 
0 0 


+n" f Ki-0 f K(ti—t)--- 
0 0 


tn—2 


»4 K (tn—-2—tn—1) 


X F(tn-1)dtn-1- + -dtd@+---. (17) 


Why one would seek the Liouville-Neumann 
series method as a solution to the equation is not 
immediately apparent, whereas the operational 
calculus yields precisely the same solution at 
once. 

Taking the transform of both members of 
eq. (2) yields 


(s) =f(s) +AP(s)R(s), (18) 
whence 
f(s) 
o(s) =. (19) 
1—Ak(s) 


When | A&(s)| <1, as is usually the case, Eq. (19) 
may be written as 


(s) =f(s)(1+AR(s) +A2R2(s) 
$---+NM(s)+-++), (20) 


or 


$(s) =f(s) 5 A"k(s). 


n=0 


(21) 


The solution (¢) is found by taking the inverse 
transform of both sides of the equation. The 
result is 


v= Fa f K(t—6) F(@)d6+---. (17) 


The solution obtained in this fashion is more 
apparent than those obtained without the use of 
the operational calculus. 


2. Kernel Expressed as Sum of Exponentials 


In physical problems, quite often the kernel 
can represent tabulated data at regular intervals. 
It is shown that when such is the case, the data 
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may have an analytical representation as a sum 
of exponentials, or 


K (t) =A ye"! +A 26%! + + -- +A ne! +->-, (22) 


where the constants are chosen so as to give the 
best representation of the numerical values.’ 
When (22) is substituted in (2), Eq. (19) is ob- 
tained where 





m A 
k(s) => — (23) 
n=1($—A,p) 
Equation (19) then has the form 
f(s) 
o(s) = (24) 





1—dv > (A,/s—an) 
n=1 
The inverse transform of $(s) calculated by the 
inversion integral gives the solution to the 
problem. Whittaker gives a solution to the 
problem by a very involved procedure, the 
method of which is extremely difficult to follow. 
March and Weaver,’ also without the use of the 
operational calculus, arrived at a solution to the 
problem by assuming that the solution would be 
given as a sum of exponentials and then evaluated 
the unknown constants. The inversion integral, 
applied to (24), gives a more direct solution. The 
solution is obtained with the aid of Cauchy’s 
residue theorem. 
f(s)e*'ds 


1 ytix 
oo=t f 
2n1 y-12 m Hs 
(1 a D ) (25) 
n=15—An 


™m 
=> p» exp(a,t), 


n=l 








where p, is the residue of the denominator at the 
poles a, and a, is a zero of the denominator. The 
zeros of similar equations have been studied by 
Raleigh.® It is well known that the residues are 
given by the ratio f(an)/q’(an), where q'(an) 
means the derivative of the denominator at an, 
and f(an) is the value of f(s) at a,. Equation (25) 


7E. T. Whittaker and G. Robinson, The Calculus of 
Observations, (Blackie and Son, Limited, London and 
Glasgow, 1937), pp. 369-371. 

8H. W. March and W. Weaver, Phys. Rev. 31, 1072 
(1928). 

® Lord Raleigh, Theory of Sound (MacMillan Company, 
Ltd., London, 1896), Volume 1, p. 334. 
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is identical with the solution obtained by The above equation may be rewritten as 
Whittaker and by March and Weaver. 





f(s) 

3. Kernel Expressed as a Polynomial ¢(s) = a ee ee en (28) 
When the nucleus is expressed as a polynomial, (- -»1) ( - rs) tee (--a.) 
s s s 
K(t) =a,t"+a,_1t" “t4. ° -+d9, (26) 
Be. (89) becenes where the \’s are the zeros of 
n! (n—1)! 

a f(s) _ (27) 1—a,—— —a,_1—_—_ : :- = 0. (29) 

n! (n—1)! ails il 

1—¢,—— —46,-1————_— °° ; 
ght s” An alternate form of (28) is 
s"f(s) 

$(s) = —_—____——. —- ——., (30) 


1 1 1 
Cayman a(t )(0-2)---(2-4) 
Ae Ae An 


Li FO (t)} +s"! F(+0) +s"? F'(+0)+--- 
$(s) =—— $$$, (31) 


-wefnale- 2-2-2) 


~ 


or 


When the numerator is of lower degree than the is known and, consequently, (32) becomes 


denominator (this is the only case discussed 
Fer) (\)i—! 


> > arive i > ba (n) r m 
here), the nth derivative of F(t), namely, F‘™ (2), &(t) =e © —+ 5 At, (33) 
vanishes and the numerator becomes a poly- n=i(r—mn)!(n—1)!  i=1 
nomial in s whose highest term is not greater than ; 
s*-!, The inversion integral of ¢(s) yields where the first term on the right corresponds to a 


factor (s—b)’. 
@(t) = Agt+A e"!'+A ce"'+---+A,e"", (32) 


where b;=1/), (i.e., the zeros of the denominator 4. Kernel Expressed as ite Series 


in (31)). The A’s are determined as in (25). When the kernel is expressed as 
Equation (32) was derived assuming that there : 
: + . . . Ki K,t" 
are no repeated roots. The solution of inversion K(t)=Ko+Kit+ oe a a 
integral with repeated roots in the denominator 2! n! 





the solution of the Volterra integral equation of the second kind may be written in a form similar 
to (33) or it may be expressed in a different form as shown below. -Equation (19) is 
f(s) 


o(s)= . | ae (35) 
1—K¢ /s—K,'/s*—K¢/s*§—---—K,_1'/s*—->: 





where K,’=\K;. Equation (35) may be rewritten in the form below: 


f(s) 
$(s) = (36) 
1—K,/p—K;'p?—Ky'p®— ---—Kes'p*—--- 





where p;=1/s;. The zeros of the infinite series in the denominator may be found by an extension of 
Graeffe’s root squaring method." 


” Reference 7, pp. 118-123. 
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When the factors are found, (36) becomes 


o(s) = 





or 


f(s) 
(37) 
(p—a1)(p—a2)---(p—an)-- 
Sls) 
ple odmaia (38) 


$(s) = 


(—1"C a] pie 


a; 


When conditions are such that each factor in the 
denominator can be expanded as (1—x)~', then 
(38) is written as 


f(s) 1 
aaa’ as 1" — = 
Ila: 
i=1 ; 
(ores 
a; a, 
Gn Qa" 
whence 
f(s) 1 
¢(s) =—____- — _ 
fa fhtse: 
Ie. 
X(1+A1pt+Aop?+---+Anp"+--:), (40) 


where the A’s represent the sum of the coefficients 
of like powers of p. Since p is 1/s, it is apparent 
that the inverse transform of ¢(s) in (40) is given 
by 


0) = BF+Bs f F(0)d6+ B, 


xf f F(\)ddd6-+ - - 


(41) 
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=) pre (lg)** ya’ - 
Ge 





An ‘alternate form of the solution i is obtained by 
writing (40) as 


6(s)=M( Jo) +- Ore 


5 


1 f(s) 
+ "4. 


s s* 


(42) 


Then 


o() =| F+ f K((-0) ods, (43) 


where 


i? t3 
K(t) = (144 es lees’ 


t” 
$e+d— te: ), (44) 


nN: 


Equation (44) is identical with the solution ob- 
tained by Whittaker. 


REMARKS 


Some of the examples have been given which 
show the utility of the operational calculus, based 
on the Laplace transform, in solving integral 
equations. The technique has the advantage over 
other proceduresin that the operational calculus 
gives the solution by an automatic method. 
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Response of Linear Resonant Systems to Excitation of a Frequency Varying Linearly 
with Time’ 


GuNNAR HOK 
Wesleyan University, Middletown, Connecticut 
(Received August 13, 1947) 


A general solution of this problem is obtained by means of Laplacian transforms. The 
resulting complex function is evaluated numerically, and universal response curves are pre- 
sented in order to facilitate the application of the solution to simple or complicated electrical 
networks as well as to other resonant systems represented by equivalent networks. 





ANY modern instruments are designed to 
plot on a cathode-ray screen the amplitude 
of a vibration or oscillation as a function of the 
frequency of excitation. Such instruments incor- 
porate narrow-band resonant systems that in 
operation are exposed to forces of varying fre- 
quency; usually the frequency varies linearly 
with time, synchronously with the left-to-right 
motion of the spot on the cathode-ray screen. 
The author first met this problem in connection 
with spectrum analyzers for study of the charac- 
teristics of microwave oscillators developed dur- 
ing the war. Later, the design of an instrument 
to portray the admittance spectra of piezoelectric 
resonators gave the impetus to an attempt to find 
a formal solution of this problem. 

Statements that such a solution has not been 
known have appeared in the current literature.? 
However, after this work had been essentially 
completed and presented at the Washington 
meeting of the Physical Society, the author’s 
attention was directed to an earlier paper on a 
similar problem in mechanical engineering.* Al- 
though for this reason complete novelty cannot 
be claimed, it is believed that a comprehensive 
study of this problem may be of value because of 
the frequent applications in present electronic 
measuring techniques. 

+ This investigation has been carried out in connection 
with Contract No. W28-003 sc-1556 with the U. S. Signal 
Corps. The fundamental solution of the problem was given 
in the second semi-annual report of July 9, 1946, under that 
contract. A brief paper on this subject was presented at 
the meeting of the American Physical Society in Washing- 
ton, D. C., May 1, 1947. 

2 E. M. Williams, “Radio-frequency spectrum analyzers,” 
Proc. I.R.E., Vol. 34, 18-22 (January 1946). 

3F. M. Lewis, “Vibration during acceleration through 
a critical speed,” Trans. of the Am. Soc. Mech. Eng. 
APM 54-24, 253 (1932). The author is indebted to Dr. 


H. Ekstein of the Armour Research Foundation for this 
and one later reference (7). 
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The following study is limited to such cases 
where the resonant system can be represented by 
a lumped-constant equivalent network (impe- 
dance operator Z in Fig. 1a). We consider only 
systems with relatively low losses, so that the 
resonance curve of each mode of vibration of the 
system never extends over more than a few per- 
cent of the frequency. This is not a necessary 
limitation; the mathematical procedure can be 
easily generalized to include lossy systems and 
large frequency variations. In a very large class 
of applications, however, this limitation offers 
the advantage that universal resonance curves or 
response curves plotted in dimensionless variables 
can be applied with very good approximation.‘ 

The frequency variation is not necessarily 
periodic; in fact, photographic recording is very 
often carried out with one single sweep of the 
frequency through the admittance band of the 
resonator. It is therefore logical to attack the 
subject as a problem in transient circuit analysis. 

It is convenient to represent the excitation by 
a complex vector 


E exp[je(t) ]=E expLjwott+e’)], — (1) 


of a constant amplitude E and an instantaneous 
angular frequency 


w=de/dt =wo+2et. (2) 


The differential equation for the simple circuit 
shown in Fig. 1b is then 
Pp" lds ; * idt ; 
el L—+Ri+ f —=E exp[j(wol+et?) ]. (3) 
dt ra 
Before proceeding with the solution it may be 
worth while to note an interesting equivalence. 


4F, E. Terman, Radio Engineer's Handbook (McGraw- 
Hill Book Company, Inc., New York, 1943), p. 136. 
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Let us change the time variable according to 
the equation 


t=7— (€/wo)r’. (4) 
Remembering that ¢ is always very small com- 
pared to wo, we can also write 
T>t+ (€/wo)l*, (5) 
and Eq. (3) is transformed to (6) 
L di 
—__—_—_§—:—+Ri 


1—2(€/wo)r dr 





* i(1—2(€/w)7) 
+f : dr=E exp[jwor]. (6) 


--00 


The equivalence of (3) and (6) means that the 
same solution will be obtained whether the reso- 
nance frequency of the system or the frequency 
of excitation is varied, at least if the relative 
increments of inductance and capacitance in the 
former alternative are equal. This point may 
seem trivial but the author has heard it seriously 
questioned and believes that it well justifies its 
place here. 

The solution of the general problem presented 
in Fig. 1a can be obtained conveniently by the 
method of Laplacian transforms. The differential 
equations are written down for all meshes of the 
network Z, and the transform is calculated for 
each term in all the equations under the assump- 
tion that the network is initially at equilibrium. 
The system is then solved for the transform of 
the current 

é 1 
i aki casipieRelineienie 


Z(p) Z(e) 
+00 
x f E exp[j(@oT +eT?) ]exp(—pT)dT. (7) 


If Z is a passive lumped-constant network, 
Z(p) is the quotient of two polynomials of p, and 
its reciprocal can be expanded according to the 
familiar methods of operational and transform 
calculus 


1 a) 3 ete) 
Z(p) h(p) —h'(p,)(p—p,) 


where ~1---p, are the roots of the equation 








(8) 
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h(p) =0. We consider only circuits for which this 
equation has no multiple roots.® 
Equation (7) thus becomes 


" (Py) 
_ a 
v=1 h'(py)(p— py) 





a0 
x f Eexp[j(woT+eT?)exp(—pT)dT. (9) 


The inverse transform is obtained as the sum 
of the Bromwich-Fourier integrals formed from 
the terms of the right member of Eq. (9). 








n 1 g(p,) f exp(pt) 
u(t) = — 
m= 2nj h'(p) J, p—p, 
+0 
xf E exp[j(woT + eT”) ] 
is Xexp(—pT)dTdp (10) 
n g( pv) sed . 9 
-£5. J E exp[j(woT +eT")] 
1 x aa 
‘ — f expel TI lar. an) 
21] Br p—P» 


The Bromwich-path integrals are equal to their 
respective residues exp[p,(t—T)] for t>T and 
zero for t<7. We can therefore write 


n (py) 
i)= 
EG) 





f Eowliertr] 
~ Xexp[p.(t—T)]dT. (12) 


At this point the dimensionless variable is 
introduced that will give the universal response 


z— 


@) 
Fic. 1. Circuit 
ZL diagrams. 
i 
£ei* 2 Cc 
é) 


’ How to cope with multiple roots in transform calculus 
is well known and does not require special attention here 
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curves for this type of excitation: 


y=at jp=se'(joo—p,)+jeT. (13) 


When y is used outside the integration sign, 
T is replaced by t¢. 
Equation (12) now assumes the following ap- 
pearance 
i(t) = E exp[j(wot+e?) ] > a) 
vt h’(p,) je: 


‘ 
Xexp(jy”) exp(—jy*)dy, (14) 
(0— ja) 
n g(py) 
i(t) = E exp[_j(wot+e?) ] >> —— —Fr(y). (15) 
= e'h’(p,) 


The lower limit of the integral in Eq. (14) is 
not strictly in accord with the change of variable 
in Eq. (12), but it is easily seen that any point on 
the quartcircle between (* —j-0) and (0O—j~) 
can be used as a lower limit without changing 
the value of the function. 

The solution is thus obtained in terms of the 
excitation vector and a function Fr(y). In order 
to apply Eq. (15) to any specific problem we 
must have tables or graphs of this function for 
all practical values of y. The remaining part of 
this paper is largely devoted to a study of this 





function 
. . Y ° « 
Fr(y) = —j explir) f exp(—jy’)d,. (16) 
0—ju 


The integral is related to the error integral and 
Fresnel’s integrals, which have unfortunately not 
been tabulated to any appreciable extent for 
complex values of the variable.*? As this function 
is likely to find extensive use in the future a con- 
venient name and symbol are desirable. The 
author suggests the name Fresnel function or 
fresneloid and the symbol Fr(y). 

In order not to make the discussion of the 
properties of the Fresnel function too abstract 
we shall introduce this function in the form of 
response curves for the simple circuit of Fig. 1b. 
Here 


Z(p) =(L/p)(p?+26—w1”), (17) 


pi R 1 R 3 
= —$+ joi= ——+j —-—). (18) 
bs 2b ~\LC 41? 


For simplicity, wo can be taken to be =a, the 
resonance frequency of the circuit. This means 
that we arbitrarily use negative as well as posi- 
tive values of time and stipulate that ‘=0 when 
the excitation frequency coincides with the reso- 
nance frequency of the system. 


Fic. 2. Universal response 
curves for a resonant system 
excited by a linearly vary- 
ing frequency (amplitude 
only). 


*E. Jahnke and F. Emde, Tables of Functions (Dover Publications, New York, 1943), pp. 23-40. 


7]. Barkley Rosser, ‘“Theory and applications of /(? exp(—x*)dx,"” OSRD 5861. 
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Thus 


v1 = }6e—'+ je't, (19) 


V2= 46€—!— j(woe!+ 32). (20) 


According to our assumptions, the imaginary 
component of yz approaches minus infinity and 
the corresponding integral zero. The solution is, 
therefore, 


1 
i(t)=E expLitud-+e") Fr(y), (21) 


4 


The definition of y in Eq. (13) has purposely 
been chosen such that a and 8 are proportional 
to and have the same sign as the resistance and 
the steady-state reactance of the circuit at the 
instantaneous frequency. 

If y is large (@ large, e small) the following ap- 
proximate solution results 


1 1 
i(t)=E exp j(wot+e) ]-—— -——— 
2j¢L —2 jy 


Eexp[ j(wot+el?)] Eexp[j(wot+et?) ] 


R44 jetL R+jX 


oS 





23) 


The solution is thus consistent with the obvi- 
ous requirement that for very slow frequency 
variation the current at any instant shall be equal 
to the expected steady-state current at that in- 
stantaneous frequency. 

Figure 2 presents the universal response curves 
in terms of dimensionless variables.’ In other 
words, it is the amplitude of the Fresnel function 
in terms of y=a+ j8. The real part, a, deter- 
mines the character of the curve while 6 gives 
the time scale of the curve. In any design problem 
involving the kind of excitation considered here 
a is the crucial parameter. According to Eq..(22), 
a is the logarithmic decrement divided by the 
square root of twice the rate of change of the 
angular frequency (the rate of change =2e ac- 
cording to Eq. (2)). The time scale is the square 
root of half the rate of change of the angular 
frequency. The group of curves of the current in 
the circuit, as shown in Fig. 2, would be observed 
if the rate of change of the frequency were kept 


* Some of the expansions used to calculate these curves 
are quoted in the Appendix to this paper. 
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Fig. 3. Observed response of a piezoelectric resonator to 
excitation of linearly varying frequency. 


constant and the resistance varied from one 
curve to the other. If the resistance were constant 
and the frequency variation 2e were changed, the 
actual response would have to be obtained by 
dividing both coordinates in the diagram by the 
value of e! for each curve. The order of the curves 
would then be reversed, the one with the largest 
a reaching the highest maximum. This is the 
way the oscillograms in Fig. 3 have been pro- 
duced. The resonant system for which these 
curves were recorded is the fundamental mode of 
vibration of an AT-cut quartz plate excited by a 
linearly varying frequency of three consecutively 
different rates of change. The top curve shows a 
perfectly symmetrical admittance curve obtained 
at slow sweep, while the other two represent 
a-values in the neighborhood of 0.5 and 0.2, 
respectively. 

The curves in Fig. 2 give information only 
about the amplitude of the current. For most 
calculations it is necessary also to know the phase 
relations between the current and the exciting 
voltage. Curves of this phase angle or of the in- 
phase and quadrature components of the current 
are very inconvenient to plot and.to use because 
the periodic components are much larger than 
in Fig. 2, even predominant for small values of a. 
It is possible, however, to find a function ®(7) 
that has no appreciable periodic components and 
is continuous for all values of a. Such a function 
is defined by the equations 


®(y)=Fr(y) for B<0, (24) 
®(y) = Fr(y) — (jn)! 
Xexp(jy’)[1—exp(js*)] for B>0. (25) 


It should be pointed out that it is not possible 
to separate the “‘forced”’ and ‘“‘free’’ oscillations 
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in terms of simple functions. The function ®(y) 
is introduced for practical, not for theoretical 
reasons. The curves of the real and imaginary 
components of ®(y) in Figs. 4 and 5 make it 
possible to calculate the in-phase and quadrature 
components of the current from Eqs. (24) and (25). 

The physical interpretation of the response of 
the resonant circuit is simplest for the ideal case 
of R=0. The reactance of the circuit suddenly 
changes sign when the frequency of excitation 
passes through the resonance frequency. This 
produces a transient oscillation of the latter fre- 
quency and of an amplitude equal to twice the 
current at that moment. Because of the absence 
of resistance this transient is not attenuated, and 
it produces beats with the current of the excita- 
tion frequency. 

In this case, y is a purely imaginary quantity, 
and the current of the excitation frequency can 
easily be calculated from tables of Fresnel’s in- 
tegrals. The curve for a=0 in Fig. 2 shows the 
total amplitude of the current, including the 
transient as well as the ‘‘forced”’ oscillation. In 
Figs. 4 and 5 the forced oscillations are repre- 
sented by the distance from the curves for a=0 
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Fic. 4. Graph of the real part of 
the function ¢(7). 


to the 6-axis for negative 8, and by the distance 
from these curves to the dashed line for positive 
8. The dashed line represents the amplitude of 
the transient when a=0. Both the real and the 
imaginary components of the forced oscillations 
are in this case antisymmetrical with respect to 
the point 8=0, #=0, and the transient oscilla- 
tion is: 


0 D 
exp(b)| f exp(j6*)d8-+ f exp is? | 
—oO 0 


(jx)! 
=2— 
2 





exp(pil). (26) 

When the circuit has resistance (a>0) the 
reactance does not change abruptly; conse- 
quently, the transient builds up gradually. The 
transient is, of course, exponentially attenuated. 
As is indicated by dashed curves in Fig. 2, the 
exponentials extrapolated back to 8=0 still have 
the virtual amplitude z}. As a@ is increased to 2 
and higher, the beats disappear and the curves 
become symmetrical admittance curves. Figures 
4 and 5 show how the current components under 
the same conditions approach the conductance 
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and susceptance curves of the simple resonant 
circuit. A better demonstration of the asymptotic 
approach of the response curves to the admit- 
tance curves for steady state is given by Figs. 6 
and 7. In Fig. 6 the real and imaginary com- 
ponents of the current at the steady-state half- 
power points are plotted in percent of their 
steady-state values. Similarly, in Fig. 7 the cur- 
rent components at the resonance frequency are 
plotted in percent of the steady-state resonance 
current. It is interesting to note that at the lower 
half-power point the.current leads the voltage by 
45° regardless of the a-value. A corresponding 
condition does not prevail at the higher half- 
power point, where both a free oscillation and a 
forced oscillation of different frequencies are 
present. 

The current at the lower half-power point is ob- 
tained by setting 8 = —a in Eq. (21). The result is 


gmp exp (6) f exp( — 6*)dé 
' (6=av2), 


2eée*L, 


(27) 


Fic. 5. Graph of the imaginary 
part of the function ¢(7). 


® Reference 6, pp. 24 and 32. 
” Reference 6, pp. 34-35. 


VOLUME 19, MARCH, 1948 


and at the upper half-power point (8 =a) 
E 


| 


2e§L 
5 

Xexp(— 6*) yn j f exp(@")da)) (28) 
0 


Tables of these integrals are available in the 
literature.® 


The current at resonance frequency (8=0) is 





E a) 
Ih= -j exp(ja*) f exp—(ja*)da. (29) 
eL a 

Since a is real, this expression can be used for 
numerical computations with the help of tables 
of Fresnel’s integrals.’ 

Figure 6 is of particular interest because it in- 
dicates how large a has to be for the response 
curve to be symmetrical. The half-power points 
are on the steepest part of the admittance curve 
and therefore give a very sensitive criterion on 





247 








Fic. 6. Curve 1: The current (each component) at the 
time when 8 = —a in percent of the corresponding steady- 
state value. Curve 2: The imaginary component of the 
current at the time 8=a in percent of the corresponding 
steady-state value. Curve 3: The real component under 
the same conditions as curve 2. 


symmetry. At a=2 the maximum deviation of 
the current from the ideal curve is nearly 9 per- 
cent. At a=3.5 the deviation is still 2 percent. 
Figure 7 shows that the quadrature component 
at the resonance frequency is an even more sen- 
sitive criterion. At a=3.5 this component is 
still 4 percent of the steady-state resonance 
current. 

A graph like Fig. 2, preferably with curves for 
a larger number of a-values, gives a considerable 
amount of information about the response of a 
single circuit. In Fig. 8 the “apparent band 
width” has been plotted versus a as extracted 
from Fig. 2. By apparent band width we mean 
the distance in frequency between points of an 
amplitude of .707 of the maximum. Similar 
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Fic. 7. The real and imaginary components of the 
current at the instant when the excitation frequency is 
equal to the resonance frequency. 


curves could be shown for the peak current and 
the displacement of the peak from steady-state 
resonance as functions of a. 

However, a graph like Fig. 2 does not give 
sufficient information when a network is con- 
sidered that contains more than one mesh so that 
not only the amplitude but also the phase of the 
current must be taken into account when the 
components are added up. Equation (26) and 
Figs. 4-5 furnish means to find the in-phase and 
quadrature components of the current. Because 
the function (7) is continuous and non-periodic 
it is easy to plot with fair accuracy, and reason- 
ably correct interpolations can be made between 
the curves. With sufficiently large families of 
curves of this kind, the current components can 
easily be found for any values of a and 8, so that 
the currents corresponding to various roots (1, 
po, p3, etc. can be calculated and combined to 
form the solution for a multi-mesh network. 

The curves given here are few and far between 
and computed with slide-rule accuracy only. The 
complete mapping of this function is a task for a 
large-scale calculating machine. However, the 
applications usually dictate the small deviations 
from the steady-state response curves that can 
be tolerated, and the problem is only to adjust 
the a-values so that these tolerances are met. 
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Empirical figures have of course been found and 
published several years ago.” 

Some numerical illustrations of the application 
to a single-mesh circuit may be of interest. At a 
frequency of 2 Mc per second and a frequency 
variation of 100 ke per second the maximum dis- 
placement of the peak towards higher frequen- 
cies, which occurs at zero resistance, is about 
If a shall re- 
main =2, the Q of the circuit must not exceed 
5000. A circuit with a Q of 11,000 will appear to 
have a band width of 207 cycles per second in- 
stead of 180; a circuit with a Q= 22,000 will seem 
to have 147-cycle band width instead of 45. The 
peak current will in these cases be reduced to 
about 95 and 88 percent, respectively, of the 
steady-state value. These figures have been 
chosen with piezoelectric resonators in mind, 
which explains the rather high Q’s. 

Equation (12) can be generalized to apply to 
all transient frequency modulation that can be 
represented by an excitation vector of the form 


1.5€!/27 or 134 cycles per second. 





Now 
n g(p») c' Be 
i(t) = J Eeslier+s7))] 
v=l h' (pv) —w 
xexp[p(t—T) dT, (32) 
or 
£(py) 
i(t) =F expL jot SO) X —— 
v=1h ‘(p,) 
Xexp[ — eof expl¢(T) |dT, (33) 
where 
¢(T) =(joo—p,)T+jf(T). (34) 


In this way one can quickly reach the core of 
the problem which is the evaluation of a complex 
function of a certain general type. 

It is possible to expand this general function by 
repeated integration by parts. However, the re- 
sulting series will usually have a limited region of 
convergence or semiconvergence. 


, 1 
ef e-d T=—___ 















—~ jJwo— pr 
6=E exp[ j{wott+f(t)} ]. (30) f'(t) f'(—2+f"O 
ne he —s 
It is helpful to assume, as before, that (joo— Pr) (joo — py) 
SF OI+3f' OF rT oe. (35) 
+ j—_—__ ——_—_———_— 
(df /dt)<Kwyo. (31) (jwo— p»)* 
| 
| 
P APPARENT BANOWIDTH 
4f 
VE 
se 
| 
Fic. 8. The apparent band 
width of a resonant circuit to 3} | 
excitation of linearly varying | 
frequency. e 
‘ | 
: a 
/ 4 
. ow 
| ra « 
{ ave- 
° 4 2 3 





VOLUME 19, MARCH, 1948 


249 











If f(t) is an even function of ¢ the response of a 
resonant circuit will have many common features 
with the previous problem where f(t) =. 

In the extreme case of zero resistance the in- 
tegral is again simplified and 


0 D 


f expLifT) MT = f explif(T) dT =A. (36) 
L 0 


We assume, as before, wo=w; and ‘“‘steady- 
state resonance” at =0. Then there is a transient 
=2A exp(pil) (37) 
generated at !=0, and one objective of design is 
usually to adjust the Q of the circuit so that this 
transient is negligible. Other objectives may con- 
cern the phase angle between excitation and re- 
sponse, which depends not only on the circuit 
constants and on the instantaneous frequency of 
the excitation but also on the rate of change of 
the latter. Numerical evaluation of the complex 
function or an experimental investigation would 
be required to obtain detailed information on 
these questions. 

The author acknowledges gratefully the help- 
ful interest of Dr. K. S. Van Dyke throughout 
the time this work has been in progress, as well 
as a stimulating discussion of the subject with Dr. 
H. Ekstein of the Armour Research Foundation. 


APPENDIX 


The basic expansion of the Fresnel function, 
convergent for all values of y, is 
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{Gn 

Fr(y) =exp(jy’) “— — 39 
7 x 77 | 
——+j—+—-:-- 
3 °2!5 317 


=exp(j”) | 3(x/2)'+r sind 


r3 r 
—— cos3é ——— sin5@+--- 
3 2!5 


+ j(3(x/2)!—r cosé 
r® g> . 
—— sin3é+— cos5é+---)}, 
3 2'5 


(38) 


where y =r exp(j@). 

For large values of 7 the number of terms that 
have to be included to get a satisfactory approxi- 
mation becomes prohibitive. A semiconvergent 
series in falling powers of y can then be used. 





1 i 3 3°55 3-5-7 
Fr(y)2—— j—-—+j——+——- 
y 4y* 875 167’ 327° 
1 1 
~— cosé —— sin3é—— cos5é+--:- 
2r 4r° 8r® 


1 1 
+f —— sin@é —— cos30 
2r 4r3 


3 
+— sin5é+-- ), (39) 
8r® 

Other expressions obtained by contour integra- 
tion, Taylor’s expansion, etc, have also been used, 
but the two given above have been found the 
most satisfactory. 
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Simple Relations for Calculating Certain Transient Responses 


W. J. CUNNINGHAM 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts* 


(Received August 15, 1947) 


The response of a linear transmission system to a step or an impulse is related directly to its 
steady-state response to a sinusoidal signal of variable frequency. Several very simple relations 
between these two types of response for low pass systems are collected here. These equations are 
empirical in nature, but generally may be expected to be in error no more than +25 percent. 
Such equations may be useful in preliminary design work where great accuracy is not required, 
or in checking the results of exact but tedious calculations. 





I. INTRODUCTION 


HE response of a linear four-terminal trans- 

mission system to a sudden change in its 
driving signal is of considerable interest in any 
over-all analysis of the system. This response 
may be found mathematically by any of several 
standard methods,' but if the system is not 
exceedingly simple, such a calculation is compli- 
cated and tedious. In order to simplify the 
analysis so far as possible, very simple forms for 
the driving signal are commonly used. The step 
function and its time derivative, the impulse 
function, are convenient for this purpose, since 
they are mathematically relatively simple but 
provide a severe test for the system. These two 
functions are shown in Fig. 1, a and c. The 
impulse function is approximated by a physical 
pulse of finite duration, but short compared with 
the reciprocal of the highest frequency trans- 
mitted by the system. 

There exist several very simple relations be- 
tween the steady-state response of a system toa 
sinusoidal driving signal of variable frequency 
and its response to a step or an impulse. These 
relations are usually strictly true only for an 
idealized system or a limited class of physical 
systems. However, they often apply to any un- 
known system with sufficient accuracy as to be 
very useful, particularly in preliminary design 
work where a quick answer is needed, even 
though a reasonable error may be present. The 
following discussion is a collection of such rela- 
tions together with a few pertinent remarks 
about each. 

* Now at Department of Electrical Engineering, Yale 
University, New Haven, Connecticut. 


1M. F. Gardner and J. L. Barnes, Transients in Linear 
Systems (John Wiley and Sons, Inc., New York, 1942), Vol. I. 
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The relations may be obtained in the following 
way. A standard system (either physical or 
idealized) having very simple steady-state trans- 
mission characteristics is used for a mathematical 
study. The response of this system to a step or 
an impulse is calculated by any rigorous method ; 
this process is not difficult because of the sim- 
plicity of the system. Significant relations be- 
tween the steady-state characteristics and the 
response to the step or impulse are now observed. 
These relations are taken as empirical equations 
which may apply to other systems with similar 
characteristics. 

Two systems with identical frequency charac- 
teristics must yield identical time responses. Two 
systems with similar frequency characteristics 
should yield time responses which are similar. 
For the latter case it is difficult to predict the 
difference to be expected in the time responses 
from a knowledge of the difference in the fre- 
quency characteristics. In an attempt to check 
this point, the empirical equations given here 
have been tested in a wide variety of typical 
systems, both through mathematical analysis 
and experimental measurements. In almost every 
case, the equations give information which is in 
error no more than +25 percent and usually less 
than this amount. It may be expected that errors 
occurring in the application of the equations to 
typical physical systems will lie within these 
limits. 

The equations are directly applicable to low 
pass systems, in most cases with the require- 
ment of minimum phase shift.? Often they may 


2H. W. Bode, Network Analysis and Feedback Amplifier 
Design (D. Van Nostrand Company, Inc., New York, 
1945), p. 242. 
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be applied to situations in which band pass 
systems are involved, by making use of the 
analogy which exists between low pass and band 
pass systems.’ The steady-state response curve 
of the low pass system should have only a single 
pronounced peak, although this peak may occur 
at zero, or any higher frequency. The equations 
are most accurate if the curves are relatively 
smooth. 


Il. DEFINITIONS OF SYMBOLS 


Since the various equations involve many of 
the same symbols, it is well to define them all 
at one time. Reference should be made to Figs. 1 
and 2, where the meanings of the symbols are 
illustrated. In most cases it is convenient to 
express the amplitude characteristic in decibels. 





Input signal (Fig. 1, a and c).—E,—value of 
input step, Eofo—area of input impulse. 

Steady-state response curve (Fig. 2).—A,—am- 
plification at zero frequency, fa—band width, 
three decibels below zero-frequency value (c.p.s.), 
fa—band width, six decibels below zero-fre- 
quency value (c.p.s.), fa—nominal cut-off fre- 
quency (c.p.s.), P—slope of chord connecting 
points at f=0 and f=f, (db/c.p.s.), &—slope 
(always negative) of tangent at point three 
decibels below maximum of curve (db/c.p.s.), 
8—phase shift through system (radians), w—an- 
gular frequency (radians/second). 

Response to step (Fig. 1b).—£,—final value, 
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- Fic. 1.a—Step function ; b—response of low pass system 
to step function; c—impulse function; d—response of low 
pass system to impulse function. 


7E. C. Cherry, “Pulse response, a new approach to 
a.c. electric network theory and measurement,” J. Inst. 
Elec. Eng. 92, III, 183 (1945). 
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s—maximum relative overshoot, f—build-up 
time (sec.), fg—time delay (sec.). 

Response to impulse (Fig. 1d).—é,—peak of 
response, fg—time delay (sec.), ¢-—duration at 
é,/2 (sec.), t;—duration measured at base of 
main peak (sec.), ¢,—duration measured at 0.012,,, 
taken to envelope if response oscillates (sec.). 
The units in which A,, /o, Z,, and é@, are meas- 
ured depend upon the type of system under 
consideration. For example, if Ko, Z,, and é, are 
all voltages, then A, is the voltage amplification 
of the system and is merely a ratio. 


Ill. TIME DELAY OF OUTPUT SIGNAL 


In any system having minimum phase shift, 
a signal begins to appear at the output just after 
a signal is applied at the input. The output signal 
may be exceedingly small for a time, after which 
it increases rapidly, giving rise to an apparent 
time delay. This time delay depends primarily 
upon the slope of the curve of phase shift through 
the system as a function of frequency. The time 
delay for the response to a step, or an impulse, 
is defined as shown in Fig. 1, b and d. The time 
delay is given by the equation‘ 


ta= |dB/dw}. (1) 


When the phase curve is not linear, it may be 
approximated with a straight line at pass band 
frequencies, and the slope of this line taken for 
use in Eq. (1). 


IV. RESPONSE TO AN IMPULSE 


If the input signal is a very short pulse 
(ideally approaching zero duration) the output 
signal will be a pulse, more or less broadened and 
possibly containing oscillations. 

The peak value of the output signal, é,, is 
given by the equation 


én=AnEoto(2fat7.4/|*|). (2) 


The basis for this equation is given in the ap- 
pendix. The equation may not be applicable if 
the steady-state response curve falls more slowly 
than 1/f? at frequencies well above the pass 


band. 


4E. A. Guillemin, Communication Networks (John Wiley 
and Sons, Inc., New York, 1935), Vol. II, Chap. XI, 
Sections 4 and 5. 
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If the steady-state curve falls as slowly as 1/f, 
a better approximation to the peak is 


ér,=AnEoto(2rfa). (3) 


This equation is the rule of thumb often quoted 
to the effect that the peak is proportional to 
the band width of the system. In case of doubt 
as to which of these equations, (2) or (3), will 
give the better results, Eq. (2) is probably the 
more accurate. 

The duration of the output signal may be 
defined in several ways. If measured at half the 
peak of the signal, the duration, ¢,, is given by 
the equation 


te=1/(2fa). (4) 


This equation is an extension of the relation, 
4fntn=1, between the nominal cut-off frequency, 
fn, and the nominal duration, ¢,, for real sym- 
metrical pulses. The nominal cut-off frequency is 
that which, when multiplied by the zero-fre- 
quency value of the amplitude curve (ratio form) 
yields the total area under the curve. The nominal 
duration is similarly defined. The frequency, fa, 
corresponds approximately to the nominal cut- 
off, and the duration, ¢., corresponds approxi- 
mately to twice the nominal duration. 

If only the main part of the output pulse, 
near zero time, is considered, the duration, ¢,, 
measured near the base of the pulse is given by 
the equation 

ty=1/fa. (5) 


This equation is an extension of Eq. (4), based 
on the assumption that the output pulse is 
roughly triangular in shape. 

The duration, ¢t,, measured at one percent of 
the peak of the output pulse is given by the 
equation 

t,=0.16| |. (6) 


The basis for this equation is given in the 
appendix. For an oscillatory response, ft, is 
measured to the envelope of the oscillations 
rather than to the response itself. 

If the damping of the system is large, the out- 
put signal will have no oscillations, and ¢; and ft, 
will be approximately the same. If the damping 
is small, the response will be oscillatory, and ¢, 


5H. A. Wheeler, “The interpretation of amplitude and 
phase distortion in terms of paired echoes,’’ Proc. I.R.E. 
27, 359 (1939). 
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Fic. 2. Important parameters in the shape of the 
steady-state transmission characteristics of a low pass 
system. 


will be larger than t;. The intensity of the oscilla- 
tions is indicated by the amount ¢, exceeds ty. 
Since the slope, P, is | P| =6/fa=6t;, and the 
slope, ®, is || =¢,/0.16=6¢,, oscillations are 
indicated if slope |®| exceeds slope | P|. 


V. RESPONSE TO A STEP 


If the input signal is a step, the output signal 
will rise to its final value more or less rapidly, 
always requiring a finite time for this build-up. 
The final value, EZ,, is given by 


E,=A,Eo. (7) 


The build-up time may be defined as the time 
required for the output signal first to reach its 
final value if it rises at a rate corresponding to 
that at half its final value. The build-up time, &, 
is given by the equation®? 


ty =1/(2fa). (8) 


A better approximation is to use f,, rather 
than fa, in Eq. (8), although f, is not easily found 
from the amplitude curve expressed in decibels. 

If the damping of the system is small, the 
response to the step will overshoot its final 
value, which will be reached only after a train 
of oscillations. Such behavior is predicted by 
the same criterion which indicates oscillations in 
the response to an impulse. The maximum rela- 
tive overshoot, s, is expressed as a fraction of the 
final value, and is given by the equation 


—21 —21 
$s =exp -——) =exo( -). (9) 
fal@—P| fa|\®| —6 


where |#| must be greater than | P|. 


’ See reference 4, p. 479. 
7™C. C. Eaglesfield, ‘Transition time and pass band,”’ 
Proc. I.R.E. 35, 166 (1947). 
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Fic. 3. a—Actual steady-state transmission character- 


istics of a low pass system; b—output signal due to im- 
pulse at input; c—output signal due to step at input. 
Circles indicate correct values. 


Two further observations may be made, based 
on the fact that the response to an impulse is 
the time derivative of the response to a step. 
The time for the maximum overshoot is the time 
of the first zero of the impulse response, and is 
approximately (tg+/;/2). The interval needed 
for the step response to begin its rise and then 
to settle very nearly to its final value is approxi- 
mately ¢,. 


VI. APPLICATION TO EXAMPLE 


As an example of their use, some of the fore- 
going equations are here applied to the problem 
of finding the response of a typical low pass 
system to an input step or impulse. The steady- 
state transmission characteristics of the system 
to be studied are shown in Fig. 3a. These are 
the actual characteristics of two sections of a 
circuit of the type shown in Fig. 4a. 

The necessary quantities may easily be evalu- 
ated from the curves of Fig. 3a, and are as 
follows: 


fa=O0.8f1, 
fa= 1.0f1, 


|| = 12.7/f1, 
|dB/dw| =0.5/f;. 


For convenience, these quantities are expressed 
in terms of a reference frequency, fi, which is 
here equal to fa. 

If an impulse is applied to the input of the 
system, important factors in the shape of the 
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output signal are found as follows: 


ta=0.5/fi, ‘r= 1/fi, 
t.=0.5/fi, t, =2.04/f, 
és 
a TY 
A2Eotof 


These factors are used to predict the shape of 
the output signal in Fig. 3b. Since the time f¢, is 
greater than the time ¢;, oscillations in the 
response are predicted. The envelope of the 
oscillations decays almost, exponentially, and 
may be assumed to have the form k; exp(— eft). 
Evaluation of the constants k; and ke by making 
use of the known points, ¢a and ¢,, gives for the 
envelope, 11.6 exp(—3.07f,t), which is plotted 
in the figure. The predicted response passes 
through the points located from the equations 
and lies inside the envelope. It oscillates with a 
period about twice the time f¢,;. 

So long as fy is at least twice fz, it is measured 
from zero time as in this example. However, if t; 
is less than twice ¢, (as will be true for systems 
with large time delay) ¢, and ¢t; will be approxi- 
mately equally distributed on either side of fg. 
In such a case, ¢, will be measured from the 
beginning of the output signal, and not from 
zero time. 

If a step is applied to the same system, the 
build-up time is &4=0.5/f;. The maximum rela- 
tive overshoot is s=0.043 and occurs approxi- 
mately at (tg+t;/2) =1/f;. The shape of the out- 
put signal is predicted in Fig. 3c. 

Points found from exact mathematical analysis 
are plotted as circles in Fig. 3, b and c. The 
agreement between the predicted curves and the 
correct points is fairly good. 
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APPENDIX 
Source of Equations 


Equations (2), (3), (6), (8), and (9) may be obtained 
from studies of the sort now to be described. The build-up 
time is found as follows. 
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At half the final value of the response to a step, the slope 
is E,,/ty, following the definition of t,. For a low pass system 
with zero phase shift, the response to a step is given by the 
Fourier integral® 


e(t)=2 f ” (Eo/juw)Ang(f) exp(jut)df. 


In this equation, (Eo/jw) is the frequency spectrum of the 
step, and A,g(f) is the amplitude characteristic of the sys- 


tem. The resulting response is symmetrical about zero 


time, under the assumptions made, so that the slope of the 
response at half its final value is 


(de(t)/dt):.o=2AnEo " g(f df =2A nEofn=En/to. 


Therefore, ts =1/(2f,), and if it is assumed that f, is ap- 
proximately the same as faz, Eq. (8) results. 

Consider now the circuit shown in Fig. 4a. This is a 
simple low pass filter which is driven by a pentode current 
source. The amplitude characteristic for this circuit may be 
written as 


Ea 
AE, 


1 

~ [(1—x?)?#+d)2x2}" 
where E, is the r.m.s. voltage at the input, E» is the r.m.s. 
voltage at the output, A; is the voltage amplification at 
zero frequency, (A1=gmR), gm is the transconductance of 
the tube, d; is the dissipation factor for the circuit 
d, = R/(Lw;)= R(C/L), x is a normalized frequency vari- 
able (x=w/w:), and w,=(1/LC)!. At the point (x=x,), 
where g is equal to 2~ times the zero frequency value, the 
slope of the curve of g vs. x is 





g= (A 1) 








(4) =(-1/4)[4+h)(h+[4+ PF, (A 2) 


x 
where 
h=(2—d,?). 
If d,<2', the curve of g has a peak at other than zero 


frequency. At the point (x=x,), where g is equal to 274 
times its maximum value at this peak, the slope of the 


curve is 
d 2(h4[4—h?}) 
(4) -=[Pee-") ay 
x/z2=z2p 4—h? 
The voltage ratio, g, may be expressed in decibels as 
G = 20 logiog, 
in which case the slope relation is 
dG 201 d 
es ae oe. (A 4) 


dx g dx 


The value of dG/dx, at x. for dj >2! and at x, for d; <2}, is 
given the symbol ¥. The response of this circuit to an 
impulse of area Eolo, when d;<2, may be written as 
eb _exp(—dywit/2 
Ai Edw (1—d,?/4)! 


where é is the instantaneous voltage at the output. The 


) sin[eoit(1—d)2/4)8], (A 5) 





8 See reference 4, Chapter XI. 


VOLUME 19, MARCH, 1948 


peak, @, occurs at the time, tp, given by 
tan[wit,(1 —d,?/4)!]=[(4/d,?) —1}!. 


Similar equations may be obtained for the cases of d;=2 
and d,>2. 

The quantities w:f, and WY may be determined from the 
preceding equations. The ratio wit,/W as calculated for a 
wide range of values of d, is plotted in Fig. 4b. It is evident 
that this ratio is surprisingly constant, and that to a good 
approximation it may be set equal to unity. A little alge- 
braic juggling then leads to 


tp=0.16|®|, (A 6) 


which is the same as Eq. (6). 
A second type of circuit which is easy to study.is the low 

pass analog of stagger-tuned circuits. The amplitude 

characteristic for such circuits is 

Ae? oe 
|AnEa| (1-+y?")n/2’ 





Zn (A 7) 
where E, is the r.m.s. output voltage, r is the number of 
circuits in the stagger-tuned unit, # is the number of units 
in the system, y is a normalized frequency variable 
(y=f/f-), and f, is the frequency at which g:=27 (i.e., 
n=1). For the present discussion it is not necessary to go 
into detail about the actual circuit, but merely to take 
Eq. (A 7) as a mathematical response to be studied. Maxi- 
mum value (unity) for g, occurs at y=0. At the point 
(y=ya), where g, equals 274, the slope of the curve is 


o a) _ we(Qve— 1 e-em 
dy/y=va 


Qati/n) 





(A 8) 


This quantity, y, may be expressed in decibels, and the 
symbol, I’, defined as 


2 (2) 
dy ¥Y=Va 


For the case of r=1, the response of the system to an 
impulse of area Eofo is 


(A 9) 


en a (wet )(e-) 
AnEdwe (n—1)! 





exp(—w-t), 


(A 10) 


where e, is the instantaneous voltage at the output. 
For the case of r=2 and n=1, the response to the 











Fic. 4. a—Low pass system for analysis ; b—curve show- 
ing the constancy of the ratio wit,/| | ; c—curves showing 
the region of application of the empirical equations. 
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Cor- Eq 
rect (A 13) 
fa welg én én 
r n welg le 7 I IU | AnEolowe AnEolwe 
1 i 4.6 1.0 —0.35 —4.3 1.06 1.0 0.59 
1 2 7.6 0.64 —0.64 -—7.9 0.96 0.37 0.35 
1 4 11.4 0.44 —1.03 —12.5 0.91 0.23 0.24 
i 6 14.2 0.35 —1.33 — 16.3 0.87 0.18 0.18 
i 10 18.6 0.27 —1,.77 —28.6 0.86 0.13 0.14 
2 1 8.1 1.0 —0.71 —8.6 0.94 0.46 0.45 
2 2 10.8 0.80 —1.03 —12.5 0.86 0.34 0.35 
2 3 13.2 0.71 —1.23 —15.1 0.86 0.29 0.31 
impulse is 
ey ° 
= 2) exp(—w-t/25) sin(wt/2!). (A 11) 
A 1E lowe 


Similar expressions apply for other values of r and n. 

The quantities wf, and I’, and their ratio, may be calcu- 
lated from the preceding equations and collected in Table I. 
From this table it is apparent that the ratio w.f,/\T| is 
always near unity in value, and, with some algebra, Eq. (6) 
may be obtained again. 

A study of Eq. (A 5) shows that if d; is very small or very 
large, the peak of the response is given by 


én 
——_=2sf, or 2zfb, (A 12) 
A,Edo ia f 
where f, or f, (as the case may be) is the frequency at 
which g is 2+ times its maximum value. This is the same as 
Eq. (3). 





A similar study of the data of Table I shows that 
Eq. (A 12) does not apply there with very good accuracy, 
and, further, that the following equation is better. 


a 0.3(fa /f-)+0. 1 / ly | e 


(A 13) 


Values found from this equation are tabulated and may be 
compared with the correct values obtained from equations 
of the type of Eq. (A 10) and (A 11). If expressed in terms 
more convenient for ordinary use, Eq. (A 13) becomes 


4 
Cn 


amar TANG. 
A,E do f 


(A 14) 
which is the same as Eq. (2). 

The values of é, for the low pass circuit as calculated from 
Eqs. (A 5), (A 12), and (A 14) are plotted in Fig. 4c. For 
values of d; in the order of unity (i.e., values commonly 
used in filter circuits) Eq. (A 14) is evidently more accu- 
rate; for other values of d,, Eq. (A 12) is better. For the case 
of r=1, n=1, in Table I, Eq. (A 12) is better than Eq. 
(A 14) for finding é,. By correlating the shape of the output 
signal with the applicable equation, it is observed that 
Eq. (A 12) is better when the output signal (or its envelope ) 
begins with a discontinuity ; otherwise Eq. (A 14) is better. 
In general, only those circuits whose steady-state charac- 
teristics fall as 1/f above the pass band yield a pulse re- 
sponse with a discontinuity, so that Eq. (A 12) is to be 
preferred only in those cases. 

Equation (9), for calculating the relative overshoot in the 
response toa step, may be set up through empirical methods 
similar to those just described. 
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On Creep and Relaxation. II.* 


B. Gross 
National Institute of Technology, Rio de Janeiro, Brazil 


(Received August 20, 1947) 


A phenomenological theory of the elastic aftereffect can be formulated with the aid of the 
principle of superp6sition. In a previous paper the theory of the transient effects, caused by the 
sudden application of a constant load or a constant deformation, has been developed. In the 
present paper the theory of the steady-state behavior under alternating load and deformation is 
given. Relations are established between the loss factor, the storage factor, the distribution 
functions, and the Laplace transforms of the creep—and the relaxation—function. 





N an experimental study of the behavior of the 
visco-elastic body one may choose two ways 
of attack; they are: (a) The measurement of the 
transient phenomena caused by the sudden appli- 
cation of a constant load (creep or strain re- 
tardation under constant stress) or a constant 
deformation (stress relaxation under constant 
strain). (b) The measurement of the steady-state 
behavior under alternating load or deformation. 
The transient is measured easily for times of 
about 1 minute or longer; the steady-state be- 
havior involves measurements with technical and 
higher frequencies, corresponding to periods 
smaller than 1 second. The two methods there- 
fore complete themselves mutually; the first one 
is suitable for the study of the processes with the 
long relaxation times and the second one for 
those with the short relaxation times. 
The ultimate aim of the measurement is the 
acquisition of information that may conduce to a 


general picture of visco-elastic behavior and the 


knowledge of the nature of the relaxation process. 
Its realization depends on additional evidence 
beyond the simple. experimental fact. This may 
be provided (a) by an interpretative theory of the 
fundamental relaxation process,' or (b) by a 
general phenomenological principle that allows 
for the coordination of the multiplicity of experi- 
mental facts under one point of view. Such a 
principle is the assumption of the linear super- 
position of effects, which, since Boltzmann, has 
guided research in this field. Many years of ex- 
perience have shown that Boltzmann’s principle 
is correct for a wide variety of materials at low 

* Publication assisted by the Committee on Inter- 
American Scientific Publication. 

‘Such a theory has been given for instance, by M. S. 
Green and A. V. Tobolsky, J. Chem. Phys. 14, 80 (1946). 
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and medium stress and strain; even when it 
breaks down, it still remains a first approximation 
to the true fact.* 

In a previous paper*® the phenomenological 
theory of the elastic after effect has been de- 
veloped, the treatment being mainly concerned 
with the transient effects. In the present paper 
these considerations shall be extended so as to 
include the steady-state behavior under alter- 
nating load and deformation. Mathematically 
speaking, while the former treatment has been by 
Laplace transforms, the present one is by Fourier 
transforms. Irreversible effects shall again not be 
taken into account. Now this is certainly of very 
little consequence, because, at least when loads 
or deformations near the yield values are ex- 
cluded, the time constant of the irreversible 
effects is many times greater than that of the 
reversible phenomena, which play a role in the 
range of frequencies that concern us here.‘ 

The transient solution for constant load is 
characterized by one function of direct physical 
significance, the creep function ; the steady-state 
solution for alternating load is characterized by 
two such functions, the loss factor, giving the 
energy dissipation per cycle, and the storage 
factor which is a measure of the reversible 


2 Recently, the principle of superposition has been ex- 
tended so as to include non-linear effects by H. Leaderman 
(‘Elastic and creep properties of filamentous materials,” 
J. Text. Found. Washington, D. C., 159 (1944)) and E. 
Guth, Paul E. Wack, and Robert L. Anthony, J. App. 
Phys. 17, 350 (1946). One of these generalizations applies to 
creep, the other one to relaxation. It is not yet clear how 
far these laws are compatible. Here we confine ourselves 
entirely to the linear problem. 

3 B. Gross, J. App. Phys. 18, 212 (1947). In the following 
this paper will be quoted as I. 

4A treatment of the transient effects which includes 
plastic deformations has recently been given by Charles 
Mack, J. App. Phys. 17, 1086 (1946). 
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storage of energy due to the time-dependent 
mechanism. An analogous situation exists for 
given deformation. Thus for one relation in the 
transient treatment there appear two relations in 
the steady-state treatment. In theory, all these 
functions and relations are equivalent ;5 in prac- 
tice it usually is just one way of treatment that 
fits adequately for a given problem. In particular, 
in the study of the interdependence of creep and 
relaxation phenomena the use of the Laplace 
transform is generally preferable to that of the 
Fourier transform. 

The theories of the elastic after effect, of 
dielectric absorption, and of linear electric cir- 
cuits are formally equivalent, only names have to 
be changed. Load corresponds to voltage, de- 
formation to quantity of electricity ; rate of creep 
and elastic modulus function correspond to 
indicial admittance and steady-state impedance. 
Thus results obtained in one branch of the 
general theory are usually significant for the 
other branches too. The steady-state solution in 
circuit theory has been worked out in great 
detail. Thus, many of the relations presented in 
this paper are contained in some form or other in 
the existing literature, but their systematic co- 
ordination, forming a consistent description of 
the various aspects of visco-elastic behavior, is 
not easily found, yet to our belief it is of interest. 

In measurements with constant loads or de- 
formations, oscillations do occur immediately 
after the application or suppression of the load, 
but usually observations start only after these 
oscillations have died away. Then the displace- 
ment proceeds so slowly that inertia reactions 
and viscous friction losses can be neglected and 
the experiment be considered as a static one. 
Under such conditions the measured values of 
load and deformation give direct indication of the 
values of stress and strain. It then becomes 
irrelevant, if the principle of superposition is 
stated in a form in which it correlates load and 
deformation (as has been preferred in 1), or stress 
and strain. In measurements with alternating 
joad and deformation, inertia reaction can be no 
longer generally neglected. Thus, for the sake of 
obtaining the relation between load and deforma- 


5’ Many correlations between these functions have been 
considered by Gemant, Physics, 7, 311 (1936); R. 
Simha, J. App. Phys. 13, 201 (1942); T. Alfrey and P. 
Doty, J. App. Phys. 16, 700 (1945). 
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tion, the dynamic equation of motion of the 
vibrating body must be established. It can be 
written down, when one knows (a) the mechanical 
properties of the system and (b) the stress-strain 
relation. It is this relation that is governed by the 
principle of superposition. In the present paper 
we do not consider the dynamic equation but 
study exclusively the stress-strain relations, which 
follow from the substitution of Hook’s law by the 
more complicated expression represented by the 
principle of superposition. For the sake of sim- 
plicity and to keep in line with the former paper, 
we employ the same terminology as in I, calling 
attention to the fact that our steady-state solu- 
tion does not contain the inertia terms. 


I. THE STEADY-STATE ADMITTANCE 


The deformation (strain) y(t) produced by a 
given load (stress) P(t) in a visco-elastic body is 
given by 


F t . 
yin=—| Pot f P()elt—ndr|. (1) 
E - 


E is the elastic modulus (instantaneous value), F 
a form factor, and g=dy/dt the rate of creep, y 
being the creep function.* The steady-state solu- 
tion for alternating load is obtained by putting 


P(t) =Poe™'. (2) 
Substitution in (1) gives 
y(t) =G(tw) P(t). (3) 


G(iw) is the steady-state admittance of the 
system: 


G(iw) = (F/E)[(1+L (iw) J, (4) 
where L(iw) is the Fourier transform of ¢/(t). 


Decomposition into real and imaginary com- 
ponent gives . 


L(iw) = A (w) —iB(w), (5) 


and A and B are, respectively, the cosine and the 


sine transform of ¢(t). 


A(w)= f coswr¢(r)dr, (6a) 
0 


B(w) = f sinwrg(r)dr. (6b) 


6 The deduction of this equation, also used in I, may be 
found in H. Leaderman, reference 2, p. 23; the above 
equation follows from Leaderman’s Eq. (8) by partial 
integration. 
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Both A and B have direct physical significance. 
The energy dissipation per cycle is given by 
BrP? F/E; the amount of energy, which twice in 
every cycle is alternately stored and returned, is 
(1+A)0.5P°F/E; only the term 0.5P,?F/E 
would be observed, if creep were absent. We call 
B the loss factor and A the storage factor for a 
given load. It follows from (6) that 


A(0)= f o(s)dr=¥(<), (7a) 


B(0) =B(«) =0. (7b) 


The zero-frequency value of the storage factor is 
therefore a measure of the total amount of creep. 
Il. THE STEADY-STATE IMPEDANCE 


The load (stress) P(t) necessary to produce a 
specified deformation (strain) y(t) is given by 


P(t) | d 8 
j=— - 7) g(t—7)dr |}. 
-l90-f veee-nar]. &) 


—O 


g(t) is the rate of relaxation.’ The steady-state 
solution for alternating deformation is gotten 
from 
y(t) =yoe'*!. (9) 
Substitution in (8) gives 
P(t) =Z(w) y(t), (10) 


where Z(iw) is the impedance function, 


E . 
Z(iw) =—(1—L (iw) }. (11) 
F 
L(iw) is the Fourier transform of g(t), 
L(iw) = A(w) —iB(w), (12) 


and A and B are again given by integrals of the 
type. A and B are, respectively, the storage 
factor and the loss factor for given deformation. 
The amount of elastic energy storage is given by 
(1—A)0.5yeE/F. It is worth while mentioning 
the asymmetry existing between this expression 
and the former expression (1+A)0.5P.°?F/E, 
which has been given above and refers to given 
load (stress). In the former case the contribution 
of the time-dependent mechanism increases the 


7 For the deduction of (8) see reference 6, p. 24, Eq. (11). 
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energy storage beyond the value, which would be 
observed without creep. In the present case the 
time-dependent effect diminishes the energy 
storage below the value, which would be ob- 
served without relaxation. We observe, finally, 
that the first mention of the integrals (6) has 
apparently been made by E. v. Schweidler.*® 


Complex Elastic Modulus 


In the foregoing the expressions steady-state 
admittance and impedance have been used with 
a view on the close analogy existing between the 
elastic and the electric relaxation theory. But 
again attention is called to the fact that these 
functions do not represent the total mechanical 
admittance or impedance of an oscillating system 
because they contain neither the mass reaction of 
the vibrating body nor possible friction losses. If 
one prefers, one may therefore define a complex 
elastic modulus function E*(—iw)=E,+iE2, 
where E,;=E(1—A) and E;=EB. This complex 
elastic modulus corresponds with the impedance 
function, and its reciprocal with the admittance 
function. 


Ill. RELATIONS BETWEEN LOSS FACTOR AND 
STORAGE FACTOR 


An inspection of the Eqs. (6) clearly reveals the 
fact, that A(w) and B(w) are mutually dependent. 
An inversion of the Fourier transform (6a) yields 
¢ in function of A; substitution of the result in 
(6b) gives B in function of A. There follows 
finally : 


2 7” w 
B(w) =— A (8) ———_48, (13a) 
= | ave a 





yo iw B 
A(w) = f B()- dp. (13b) 
To 


B?—w? 


The integrals are principal values. A deduction 
and discussion of these equations has recently 
been given by Gross*!° and, independently, by 
Cole." Very detailed information is given by 
Bode.” In view of the numerous literature exist- 


8 E. v. Schweidler, Ann. d. Physik 24, 711 (1905). 

* B. Gross, Phys. Rev. 59, 748 (1941). 

10 B. Gross, Am. Math. Monthly 50, 89 (1943), contains 
further bibliography. 

11 R. H. Cole, Phys. Rev. 60, 172 (1941). 

12H. W. Bode, Network Analysis and Feedback Amplifier 
Design (D. Van Nostrand Company, Inc., New York, 
1945). 


259 





ing about these relations, there is here no need for 
further discussion. We only observe that putting 
w =0 one obtains the relation 


2 7” BiB) 
A(0)= f dp, 
Tyo B 


(14) 


mentioned already by Kronig." It connects an 
integral over the loss factor with the zero-fre- 
quency value of the storage factor, or, according 
to Eq. (7a), with the total amount of creep. It is 
obvious, that the same relations hold between 


A and B. 


IV. RELATIONS BETWEEN IMPEDANCE AND 
ADMITTANCE FUNCTION 


The very definition of the impedance function 
makes it obvious that this function must be the 
reciprocal of the admittance function, and vice 
versa. Thus 


G(iw) = 1/Z (iw). (15) 


In the present treatment, the deduction of the 
expression for G has started with Eq. (1), and 
that of Z with Eq. (8). Thus, from our point of 
view, the relation between G and Z is more than a 
pure definition; it represents a condition of 
compatibility between the Eqs. (1) and (8) and 
results ultimately in a relation between the creep 
and the relaxation function. Indeed, substitution 
of (4) and (11) in (15) gives 


L(iw) +L (iw) L (iw) — L (iw) =0. (16) 


Since L(iw) is the Fourier transform of the rate of 
creep and L(iw) is the Fourier transform of the 
rate of relaxation, Eq. (16) permits the calcula- 
tion of one of these functions when the other one 
is given. The substitution iw=p converts the 
Fourier transforms L(iw) into the Laplace trans- 
forms L(p), and the relation (16), between the 
Fourier transforms, into the corresponding rela- 
tion between the Laplace transforms (1, Eq. (8)). 
This illustrates well the fact that the description 
in terms of steady-state admittance and impe- 
dance signifies treatment by Fourier transforms 
and that the description in terms of the transient 
effects signifies treatment by Laplace transforms. 
But while Z(p) is a real function, L(iw) is a 
complex one; therefore Eq. (16) implies the 


1% R. L. de Kronig, J. Opt. Soc. Am. 12, 547 (1926). 
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existence of two relations between the real and 
imaginary components of Z and L. After alge- 
braic transformations one gets 


_ A(1+A)+B? 

A=— —, (17a) 
(1+A)*+B? 

S B 

B= . (17b) 
(1+A)?+B? 


A similar set of equations is obtained for A and 
B in terms of A and B. Thus when A and B are 
known, A and B can be calculated, and vice 
versa. Furthermore, it is sufficient that either A or 
B is known, and A and B may still be calculated. 
This follows from the existence of the direct 
relations (13) between A and B. Here only the 
relations between B and B shall be considered in 
detail. For reasons of the following we make the 
substitutions 


B=V/s, (18a) 
B(\/w) =1b(w), B(/w)=rb(w), (18b) 
and get 

P b 

b = —_—_——— —— —-, (19a) 
" ds 7 
+f b(s)- + 2b? 
0 s—@ 

b 

b=- (19b) 


*. Ff ; . 
[ — { b(s)— +7°b? 
“0 S—-@ 

These relations allow for an interesting interpre- 
tation. In I the distribution functions for creep 
(p) and for relaxation (f) have been introduced. 
Now it is well known from the analogous electrical 
case, that the shape of the loss factor curve is 
similar to that of the corresponding distribution 
function. This similarity may have appeared to 
be somewhat fortuitous; now it becomes evident 
that it must exist of necessity. Indeed, a com- 
parison of the above relations with Eqs. I (27), 
which give the relations between p and p, shows 
the 2 sets of equations to be formally identical. It 
follows that every pair of functions b and 6 must 
serve as a pair of distribution functions, and 
conversely every pair of distribution functions p 
and f must fulfill the condition necessary for 
being considered as a pair of loss-factor functions. 
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V. RELATIONS BETWEEN THE LOSS FACTOR, THE 
STORAGE FACTOR, AND THE 
LAPLACE TRANSFORM 


For many problems, the knowledge of the 
Laplace transforms L(p) and L(p) is of interest. 
When both A and B are known, L(iw) =A —1B 
can be calculated and L(p) follows by the substi- 
tution iw=p. But again it can be shown that 
L(p) can still be calculated if either A or B is 
known. In the following a very simple deduction 
shall be given. For more detailed considerations 
we refer to a paper by Nijenhuis and Stumpers.'* 
Multiplying Eq. (6a) on both sides by 1/(w?+ ’), 
integrating over w from 0 to ~, and changing on 
the right side the order of integration, one gets 


* A(w)dw . * COSWT 
f nt f otal f ——de. (20) 
0 w'+p? 0 0 w+p 


The cosine integral is known, giving 


* COSWT —"e 
J dw = —- 15 (21) 
0 w*+ p? 2p 





Substitution of (21) in (20) shows, that the right- 
side term is identical with wL(p)/2p, and 
therefore 


D 





2 ¢” 
L(p) =- f A(w) dw, (22a) 
"9 





w+ p? 
se w 

Lip)=— f B(w) dw, (22b) 
T#9 w+ p? 


the second equation being obtained in a way 
analogous to that given above. A set of analogous 
equations connects A and B with L. 

These equations may find various applications. 
It has already been mentioned, that in the short 
time interval the measurement of the storage 
factor and loss factor is performed with much less 
difficulty than that of the creep or relaxation 
function. The result of the measurements with 
alternating load or deformation may then con- 
veniently be used for computation of the creep or 
relaxation function in the short time interval. An 

“4 W. Nijenhuis and F. L. Stumpers, Physica 8, 289 
ts See F. Hurwitz and R. Courant, Funktionentheorie 
(Julius Springer Verlag, Berlin, 1929), p. 312. The integral 


(21) is obtained by integrating the integral given by 
Hurwitz and Courant. 
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inversion of the Fourier transforms (6) gives 
directly g(t). But when the analytic expressions 
for A or B are not simple, then it may become 
more advantageous to calculate first the Laplace 
transform L(p) by means of one of the Eqs. (22). 
The inversion of the Laplace transform can be 
carried out by modern operational methods and 
become less difficult than the Fourier inversion. 

The Eqs. (22) can supplant the direct relations 
(13) in computations of the storage factor from 
loss-factor measurements or vice versa. To this 
purpose we observe that these equations obvi- 
ously are equivalent to 


L(pe*'*"") =A (p)¥1B(p), (23) 

or, written differently 
A(p) =ReL(pe***), (24a) 
B(p) = +ImL (pe*‘*!2). (24b) 


(23) is really identical with Eq. (5), with which 
we have introduced the factors A and B, only for 
+7 now we have substituted e* '*’?. This substitu- 
tion becomes necessary because in making trans- 
formations between L(p) and L(iw) we really 
consider one function, L(z), of a complex variable 
z; taken along the positive real axis, L(z) gives 
what here has been called the Laplace transform 
and taken along the imaginary axis it gives what 
has been called the Fourier transform. In the 
complex plane L(z) is found to have ramification 
points. To avoid ambiguities one has to chose one 
determination of the function. Along the positive 
real axis, L(z) must be real. Thus one is led to 
chose the phase of z in the interval between —7 
and x. An example for such an application of 
Eqs. (22) is given at the end of the paper. 

The technique of measurement frequently 
allows for the simultaneous determination of loss 
factor and storage factor. One may then perhaps 
doubt the usefulness of formulae containing only 
one of these functions as compared with the much 
simpler relations containing their sum or differ- 
ence. Now the result of the measurement is a set 
of numerical data which may be represented 
analytically by suitable expressions A(w) and 
B(w). But if these expressions shall constitute the 
starting point for further mathematical develop- 
ments, the functions A(w) and B(w) cannot be 
chosen independently of each other because then 
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Pets 








they would not satisfy the relations (13), the sum 
and difference A+B not being a function of the 
argument ww and not constituting a function 
L(iw). It follows, that one is free in the choice of 
one function only, and for this reason one needs 


the “‘single’’ relations of the theory. 
£ \ 


VI. DETERMINATION OF THE DISTRIBUTION 
FUNCTIONS 


In I, the distribution functions p(p)dp for creep 
and ~(p)dp for relaxation have been formally 
introduced by the relations 


g(t) -{ e~”'p(p)dp, 


(25a) 


a(t) = f e~?*a(p)dp. 


“oO 


(25b) 


A method has been given for the determination 
of p and f from creep and relaxation data. From 
the present point of view, A(w) and B(w) are the 
given function. Thus the problem arises how to 
calculate the distribution functions from the loss- 
factor and storage-factor data. This problem has 
first been considered and solved by Fuoss and 
Kirkwood.'® Here we shall give a different deduc- 
tion based on the foregoing developments. For 
this purpose multiply both sides of Eq. (25a) by 
coswf and integrate over ¢ from 0 to «. Changing 
the order of integration on the right side, it 
follows that 


D 


f y(t) coswidt 


-{ p(pdp | e~?' coswidt. (26) 


The left-hand integral is equal to A(w), the inte- 
gral over the exponential function is known. 
Making a similar integration with sint, there is 


a 


p 
A(w) = f ai (27a) 
w*+ p? 


0 


Biw)= f oD) —dp.7 —— (27b) 
0 +7 


16 R. Fuoss and J. G. Kirkwood, J. Am. Chem. Soc. 63, 
385 (1941). 

‘7 Equations of this type have been already derived by 
E. v. Schweidler, reference 8. Compare W. A. Yager, 
Physics 7, 434 (1936), Eqs. (5b) and (6b). 
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A comparison of these equations with the Eqs. 
(22) shows the formal identity of (27a) and (22b), 
and (27b) and (22a). But in the former case, the 
functions under the integral were known and the 
left-hand functions unknown; now the contrary 
happens, Eqs. (27) representing integral equa- 
tions for p(p). The method of solution is indicated 
by the analogy existing between the sets (22) and 
(27). It has been shown, that (22) is equivalent 
with the relations (23) and (24), in which the 
connection between A and B on one side and 
L(tp) on the other side is given. Thus we now 
suppose the existence of 2 ‘functions A(z) and 
B(z) of a complex variable z, such that along the 
positive real axis A(z) ouinchdes with the given 
A(w) and B(z) with the given B(w). Then in 
analogy with (24) one gets 


2 2 

p(w) =— ImA (we-**/?) = —— ImA (wet **/?), (28a) 
rg 7 
2 2 

p(w) =— ReB(we-'*/*?) =— ReB(wet'*’*).'5 — (28b) 
T T 


The solution obviously implies the principle of 
analytic continuation. 
The functions p and # are not normalized ; they 


_are related to the normalized distribution func- 


tion, F(r)dr, of retardation times of strain and, 
F(r)dr, of relaxation times of stress by the 
Eqs. I (14): 


BF(r)dr=—p(p)dp/p, p=1/r, 


and an analogous relation for F and p. The 
normalization constant 8 can be determined 
from the condition that the integral over F must 
be equal to 1 (I, Eq. (12)). But frequently 
another method of determination may be found 
convenient. Integrating (25a) from 0 to » and 
taking into account (29) there results 


(29) 


p= | o@dt=v(~)=A(0). (30) 


An analogous relation holds for 8. 


‘8 This solution is easily seen to be identical to the 
formulae of Fuoss-Kirkwood, reference 16. It can be 
arrived at in many different ways; f.i., by the substitution 
p?=s the.Eqs. (27) can be transformed into Stieltjes’ 
integral equation, the solution of which is known; c.f., 
I (16). 
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VI. APPLICATIONS 


We now shall apply some of the formulae given 
above to the discussion of a problem which has 
appeared primarily in dielectric theory but which, 
to our belief, is of interest also for the treatment 
of visco-elastic behavior. Fuoss and Kirkwood, in 
their paper frequently referred to,'® have repre- 
sented their loss-factor data by the law 


1 
B(w) =—_—_—_. (31) 
w*+w% 

The coincidence between measured and calcu- 
lated values is good. A little later, Cole and Cole*® 
formulated another law, that represents equally 
well many experimental data of loss-factor and 
storage-factor measurements and therefrom de- 
duced an expression for the corresponding func- 
tion g(t).2! The measurements of Fuoss and 
Kirkwood have been made with high polymers; 
in consequence of the close analogy existing be- 
tween the electric and the elastic relaxation 
effects for these materials, their results are 
significant not only for the electrical effect where 
they originally have been established, but also 
for the elastic effect. Therefore it seems to us 
worth while trying to get the function g(t), that 
belongs to the expression (31). Its knowledge 
may facilitate a comparison between the theories 
of Fuoss and Kirkwood and Cole and Cole. 
Unfortunately a treatment based on (31) is more 
complicated than one based on Cole’s expression. 
So we shall confine ourselves to the particular 
case 


:; (32) 


which admits of an easy solution in terms of 
known transcendents. We shall start with the 
calculation of the Laplace transform, because 
once this function is given, the other character- 
istic functions are obtained relatively easily. 
The Laplace transform L(p) belonging to (31) 
and (32) is gotten from (22b) ; integration gives 


2 pr} 
Lip)=—|1+( ) (o—1)| (33) 
1+p? 2 


' Reference 16, Eq. (43). For the sake of rendering as 
simple as possible the following expressions we shall omit 
the factor H(0) and, instead of w/wm, write simply w. 

Pe K. S. Cole and R. H. Cole, J. Chem. Phys. 9, 341 

41). 

"1K. S. Cole and R. H. Cole, J. Chem. Phys. 10, 98 
(1942). 


Nie 
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The Fourier transform L(iw) is obtained from 
(33) by the substitution p=we'"’?; there results 





1 
L(tw) = 
1 2 


—w- 


1 
[2—(w)w+1)]—i———._ (34) 
wt+w 


® 


The storage factor is given by the real compo- 
nent of L(iw) ; 


1 
AoA rns (35) 


The loss factor is given by the coefficient of the 
imaginary component of L(iw) ; 


1 
B(w) =——— (36) 


w+w 


This expression coincides with the given function 
(31); the coincidence constitutes a proof of the 
correctness of the foregoing because there can 
exist only one function L(iw), the imaginary part 
of which gives (31). 

The distribution functions are most easily 
gotten from the Laplace transform by the method 
given in I. Along the negative real axis we have 


. 2 {es} 
Lipe-*)=——|1+i(~) a+p)|. (37) 
i+p? 2 


The imaginary component, according to I, Eq. 





108. y¥—, 














tog t—— 


Fic. 1. I—Eq. (42) (Fuoss-Kirkwood). 
II—Eq. (45) (Cole-Cole). 
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(23) gives the distribution function for creep 


(2p)' 1+p 
p(p)dp =——— ——dp; (38) 
r 1+ )° 


the real component, according to I, Eqs. (21a) 
and (27a), gives for the distribution function for 
relaxation 


p(p)dp 
p(p)dp =-— —. 


(39) 
(14+2/ A+ p? )? +r p*(p) 


The distribution functions F and F are obtained 
making the substitution (29), but there is no need 
for explicitely writing down the results ; one may 
only add that the normalization constant 8, 
according to (30), is equal to 2, and B=2. 

The creep function g(t) is obtained by inverting 
the Laplace transform (33). Consider first the 
last term (p/2)!/(1+ 7). Since 1/(1+?) is the 
transform of sint and (2p)~} is the transform of 
(2nt)~', there is, by Borels theorem,” 


‘ sin(t—r) 


1 d 
— = f dr. (40) 
(2p)'(1+p?) dt Jo (2xr)! 


The integral can easily be expressed in terms of 





Fresnel’s integrals 
‘ sint 
sw=f ——-dt, (41a) 
0 (2nt)! 
* cost 
C(t) = dt.” (41b) 
0 (2xt)! 


It is interesting to note, that the expression (40) 
also appears in the theory of the transmission 
line; in this connection it has been discussed by 
McLachlan** with the aid of the complex in- 
version theorem. The second term in (33), ob- 
tained from the foregoing one by multiplication 
with ~, obviously is equal to the derivative of 
(40), the identification of » with d/dt being 
justified because the expression to which it applies 


% See R. V. Churchill, Modern Operational Mathematics 
in Engineering (McGraw-Hill Book Company, Inc., New 
York, 1944), p. 37. 

% Tables and expansions for the Fresnel integrals are 
given by E. Jahnke and F. _ Funktionentafeln (B. G. 
Teubner, Leipzig, 1923), p 

«N. W. sdclechion’ Ce ye Variable and Operational 
Calculus (Cambridge University Press, Teddington, Eng- 
land, 1939), p. 92. 
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disappears for ¢=0. Therefore the final result is 


2 
saaeees il sintL1— C(t) — S(t) ] 


g(t) = 
(2nt)! 
+2 costLS(#)—C(t)]. (42) 
For small values of ¢ there is 
t0: v(t) =2/(xt)}. (43a) 


For very great values of ¢ the asymptotic ex- 
pansion formula for the Fresnel integrals* gives 


g(t) = (43b) 


The behavior of the function (42) at very small 
and at very great values of ¢ is up to a constant 
identical with that of the corresponding Cole 
function (I (40), (41)). It can be shown, that this 
is not restricted to the value a=}, but happens 
for all values of a. Figure 1 gives in a logarithmic 
scale the function g(t)/2 given by (42). 

In I, it has been shown, how the generalization 
of the simple power law can lead to a very general 
type of creep function which is identical with the 
function indicated by Cole. It can be expressed 
in terms of the Mittag-Leffler function E£,,(—2). 
For m=}3 the Mittag-Leffler function reduces to 
the expression”® 


E,(—2) =exp(f)[1—(4) ], 


(t) is the error integral; this gives for the 
creep function 


l— @ ; 


t(2mt)}. 


(44) 


where 


1 
g(t) =—— —e'[1—4( V2) ]. 


(45) 
(xt)! 


This function has also been plotted in Fig. 1 
Clearly the difference between the two functions 
is small. 

For a of the form a=n/m, where n and m are 
integers, the integral (22b) can still be evaluated, 
but with increasing values of and m the 
calculus becomes increasingly tiresome. The 
creep function can without difficulty be obtained 
in series form, which allows for numerical 
computation. 

I am very much indebted to the Director of the 
Institute, Professor E. L. da Fouseca Costa, who 
made this study possible. 


25F, M. de Oliveira Castro, Zeits. f. Physik 114, 116 
(1939), Eq. (37). 
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A Determination of the Elastic Constants for Beta-Quartz 


E. W. Kammer, T. E. PARDUE,* AND H. F. FRisseL** 
Naval Research Laboratory, Washington, D. C. 
(Received September 2, 1947) 


The complete set of adiabatic elastic constants for beta-quartz have been obtained for a 
series of temperatures ranging from 580°C to 800°C. The elastic moduli were computed from 
the elastic constants and are compared with the adiabatic and isothermal values which have 


been determined by others. 





INTRODUCTION 


, experimental method for determining the 
adiabatic elastic constants for beta-quartz 
has been described by Kammer and Atanasoff! 
and values for the constants were reported at a 
single temperature of 600°C. The data contained 
in the present report are an extension of their 
work using the same experimental method to 
cover most of the temperature range in which 
beta-quartz is known to exist. 


GENERAL CONSIDERATIONS 


The general theory of wave propagation in 
homogeneous aeolotropic media as developed by 
Christoffel and others (for example, see Cady’) 
achieves a desirable simplicity when applied to 
thickness vibrations in plane-parallel crystal 
plates of infinite area. In order to employ this 
method of analysis on the vibrations generated in 
small plates of alpha-quartz, Atanasoff and Hart? 
demonstrated experimentally that the perturbing 
effect of the finite plate edges on the resonance 
frequencies, f,, diminishes as the number of nodal 
planes between two opposite surfaces of a 
vibrating plate are increased, where f, is the 
observed frequency and m the order of the 
harmonic. When the order n of the harmonic was 
made sufficiently large, the ratio f,/n approached 
an asymptotic value which was found to depend 
only on the thickness of *the plate, the mode of 
vibration, and the elastic properties of the 


* Now at the University of North Carolina. 
. a“ Cornell Aeronautical Laboratory, Buffalo, New 
ork. 
1E. W. Kammer and J. V. Atanasoff, Phys. Rev. 62, 
395 (1942). 
2 W. G. Cady, Piezoelectricity (McGraw-Hill Book Com- 
pany, Inc., New York, 1946), pp. 104-106. 
asdiy. Atanasoff and P. J. Hart, Phys. Rev. 59, 85 
941). 
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medium. Under these conditions the plate be- 
haves essentially as a two-dimensionally infinite 
plate, and it is permissible to apply the Christoffel 
method. 

The theoretical treatment used in this report is 
identical with that outlined in the work of 
Atanasoff and Hart* on alpha-quartz modified as 
explained in the paper of reference 1 so as to 
apply to the hexagonal holoaxial symmetry (Dz 
Schoenflies) present in beta-quartz.‘5 This crystal 
class has a single sixfold axis and six twofold axes 
perpendicular to it so that there are five inde- 
pendent elastic constants. The matrix in Voigt’s® 
notation has the following form: 


Cu Ci C13 O O 0 
Ci2 Cir C13 0 0 0 
Cis Ci3 C33 0 0 0 
0 0 0 C44 0 0 
0-0 O O Ca 0 


0 0 0 0 0 4(€41—C12). 
The secular equation for the plate orientations 
used in this study becomes: 


| 
| 


Cu’ — K? 0 Cis. 
0 Ces’ — K? 0 =(). (1) 


Co’ 0 C55 — K? 


In this expression K? = 4pd?(f,,/n)*, where p is the 
density and d is the thickness of the plate. Here 
a double-subscript notation for the elastic con- 
stants has replaced the four-subscript notation 
used in reference 3. The primes indicate that the 
constants are referred to the rotated axes. The 
rotated axes, x,'x2'x3’, are related to the principal 


*W. Bragg and R. E. Gibbs, Proc. Roy. Soc. London 
A109, 414 (1925). 

5R. W. G. Wyckoff, Am. J. Sci. 11, 112 (1926). 

®W. Voigt, Lehrbuch der Kristallphysik (Teubner, 
Leipzig, 1928), pp. 585-586 and 830-831. 
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TasLe |. Dimensions and orientations of the quartz 
plates at room temperature. The angle, theta, is measured 
between the x; and x,’ axes, being a rotation about a 
common x. axis for both coordinate systems except for 
the Y-cut. The x,’ axis is in the thickness direction through 
the plate. Lateral dimensions of all the plates were nearly 
alike, averaging about 2.3 centimeters. 





oO 


4 30° 45° 60° (Z*cut) Y-cut 


Thicknessincm 0.4549 0.4549 0.4549 0.4333 0.4445 





x1 Xe X3 
| X1 ay 0 43 
wl o lal o |, 
Xa ne —dA13 QO | d33=dy1 


where @;, etc. are the direction cosines. Equation 
(1) is limited to plate orientations for which the 
x2’ axis is parallel to the x2 direction. The x,’ axis 
is along the thickness direction. The primed 
constants can then be expressed in terms of the 
principal elastic constants by the following 
equations: 


Cy = Ay‘. +a 13°C33 + 20117013°C13 +4011°013°C55, (2) 


Coe = 413°Cqa t+ O11°Coe, (3) 
Coe. = 4171 3°C11 +1101 3°C33 
— 244171 3°13 + (dir? —@13")*C55, (4) 
Cie = Co = — 041°01 3611 +1101 3°C33 
+ (41°13 — 11013") C13 
+2 (a4;°@13 — 11413") C55. (5) 


In the above equations the equality between C44 
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Fic. 1. The ¢;; adiabatic elastic constant for beta-quartz 
as a function of temperature. 
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and cs, has not been used since c4q4 and C55 are 
coupled differently by means of the piezoelectric 
constants to the applied fields. In this way it is 
possible to identify. the modes which are excited 
if the direction of the field is known, as will be 
explained below. Also cs is measured as a single 
constant and the equality ces =3(¢11—C¢12) is not 
used until later. 

When the expressions for the c’ constants are 
substituted into Eq. (1) and expanded, the rela- 
tion between the elastic constants, direction 
cosines, and K? becomes: 


[ais'Cast+11°Co6 — K*] 
x [ai7ai3e13"7 + 2011°13°C55C13 
= 411° 3°C11€33 + (Gi PK? —11'C55) C11 


+ (d13?K? — 13455) C33 — (K* — C55) K? ] =0. 


, 


(6) 


There are five independent unknown elastic con- 
stants in this equation and three predicted modes 
of vibration indicated by the three possible values 
for K*. The process of determining the unknown 
constants in Eq. (6) is begun by selecting one 
plate with the direction cosines a,,=0, and 
a,3;=1. This has been called a Z-cut. Under these 
conditions terms with ca, and cs5 are factorable, 
giving 
(C44— K*) (C55 — K*) (C33 — K*) =0. (7) 
The two piezoelectric constants for beta- 
quartz,’ ¢14= —é25, permit modes involving the 
C44 elastic constant to be excited by a field in a 
direction x;, and, correspondingly, the field must 
have a component in the direction x2. when 
driving modes dependent upon the ¢;5 elastic 
constant. Hence it is not possible to excite the 
mode for ¢33; in Eq. (7) by means of coupling 


ELASTIC CONSTANT cs 





625. 650 675 700 725. ~~ 750 775 800 
TEMPERATURE -DEGREES CENTIGRADE 

Fic. 2. The c33 adiabatic elastic constant for beta-quartz 

as a function of temperature. 





7™W. Voigt, Lehrbuch der Kristallphysik (Teubner, Leip- 
zig, 1928), pp. 830-831. 
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through a piezoelectric effect. The remaining two 
modes can be driven by this means. Additional 
verification of the value for c4, was obtained from 
a Y-cut. For a plate of this orientation the x,’ 
axis is directed parallel to a principal x2 axis and 
the secular equation reduces to: 


(Ces — K?) (€11 — K”) (C44 — K?) = 0. (8) 


The value of c44 now being known, it is possible 
to determine cs from the linear factor of Eq. (6). 
However, it is necessary to select the proper mode 
of vibration for computation of the K? which 
occurs in this factor. This can be done by choosing 
the direction of the applied alternating field so as 
to utilize the appropriate piezoelectric constant 
for coupling with cq. The frequency of vibration 
associated with the factor containing Cg. in 
Eq. (6), will appear when the field is applied so 
as to have a component along the x; principal 
axis. The constant ¢s5,; which occurs in the 
quadratic factor requires that the field be applied 
along the x2 principal axis in order to excite the 
modes associated with this factor. 

The terms in the quadratic factor of Eq. (6) 
can be rearranged as is shown below: 


11°43" | 13° + 2€55C13 — C11€33} 
+ (17°K?—ay1'¢55) C11 
+ (a13?K? — a13'C55) C33 = (K? — C55) K°. (9) 


lf three or more independent equations of this 
form can be found it is possible to eliminate the 
terms in the braces and solve for cy; and ¢33. The 
elastic constant ¢;3 can then be obtained from any 
one of the original quadratic factors by use of the 
quadratic formula or numerical approximation 
after substituting the values for cy, and ¢33. The 
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Fic. 4. The cy. adiabatic elastic constant for beta-quartz 
as a function of temperature. 


390 





365 


& 
° 





ELASTIC CONSTANT % 
w 
~_ 
wo 





570 600 620. 640 660 680 700 720 740 +760 780 800 
TEMPERATURE —DEGREES CENTIGRADE 


Fic. 5. The ¢4, adiabatic elastic constant for beta-quartz 
as a function of temperature. 


elastic constant, C¢¢, is related to cy; and C1. by the 
equation : 


Coe = 3 (C11 — C12). (10) 


Since ces and ¢,; are now available, ci2 can be 
computed. This completes the solution of Eq. (6). 


EXPERIMENTAL 


The apparatus used in the present experiment 
was improved in several respects over that 
employed for the earlier work.' A General Radio 
Type LR frequency standard with interpolation 
oscillator gave greater precision and speed in the 
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Fic. 6. The ces adiabatic elastic constant for beta-quartz 
as a function of temperature. 
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Fic. 7. The sy; adiabatic elastic modulus for beta-quartz 
as a function of temperature. Isothermal values determined 
by Perrier and de Mandrot are also shown. 


measurement of the frequencies of the plates. 
Changes were made in the construction of the 
furnace to permit closer control of the tempera- 
ture. The electrode assembly holding the quartz 
plate in the furnace was redesigned so as to have 
greater resistance between the electrodes at high 
temperatures. This increased the sensitivity of 
the detector and permitted measurements to be 
made on harmonic resonances of higher order 
than previously attainable. 

Five plates were used to obtain the frequency 
data. These plates were prepared by the National 
Bureau of Standards. The orientations were 
measured with an x-ray goniometer. The quartz 
material was examined for flaws and twinning by 
optical methods developed at the Bureau and the 
final surfaces were etched. No conducting ma- 
terial was deposited on any of the surfaces. The 
dimensions and orientations are given in Table I. 

The uncertainty in the orientation angle @ as 
given in Table I| is about 2 minutes of arc. These 
specifications apply to the plates at room temper- 
ature. The orientation and thickness change 
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Fic. 8, The s33 adiabatic elastic modulus for beta-quartz 
as a function of temperature. Isothermal values determined 
by Perrier and de Mandrot are also shown. 
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slightly when the plates are elevated to the tem- 
peratures necessary for beta-quartz. Corrections 
for these changes in orientation and thickness 
were made in the same manner as is described in 
reference 1. They are based on the thermal 
expansion coefficients measured by Jay.* The 
maximum correction in orientation angle is en- 
countered with the 45° plate and is about 12 
minutes of arc 

The density of beta-quartz used in the calcula- 
tions for the elastic constants was estimated from 
the density of alpha-quartz at room temperature, 
2.648 g/cm*, and the expansion coefficients 
vy = 0.01733 and y,=0.00998. Thus: 


Ps = pa/(1+7y)?(1+7:). (11) 


This yields a value of 2.533 g/cm’ at 600°C. Ac- 
cording to the expansion coefficient data obtained 
by Jay and also the direct measurements of Day, 
Sosman, and Hostetter,’ the density of beta- 
quartz may be considered to remain constant for 
the temperature range under consideration in 
this experiment without introducing serious error. 
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Fic, 9. The si3 adiabatic elastic modulus for beta-quartz 
as a function of temperature. 


Beta-quartz is abnormal in that the density in- 
creases slightly with increasing temperature. 
RESULTS 

Frequencies for ten modes of vibration were 
measured for the five plates. One mode from each 
of the Z- and Y-cuts was used to determine C44. 
The 30°, 45°, and 60° plates each furnished one 
mode from which cs could be computed with the 
C44 constant being known. Five more modes were 
observed in the same three plates to furnish an 
over-determination of ¢1;, C33, and C3. 

8 A. H. Jay, Proc. Roy. Soc. London A142, 237 (1933). 


9A. L. Day, R. B. Sosman, and J. C. Hostetter, Am. J. 
Sci. 37, 16 (1914). 


JOURNAL OF APPLIED PHYSICS 

















A numerical check on the consistency of the 
experimental results was made by evaluating ¢1; 
from each of the five equations after ci; and ¢3; 
were known. At 580°, four equations yielded 
values of ¢,3 consistent to within 0.5 percent while 
the fifth gave a value about 20 percent higher 
than the mean of the other four. For temperatures 
above 600°C, the variation of any value from the 
average at the temperature in question was less 
than 0.5 percent. The change in frequency of 
vibration of the plates was large for small tem- 
perature variations near the transition point. 
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Fic. 10. The siz adiabatic elastic modulus for beta-quartz 


as a function of temperature. 


This increases the difficulty in ascertaining the 
frequencies in this region and perhaps largely 
accounts for the major inconsistency in the 
results for ¢i3 at 580°C. 

The values of the elastic constants are presented 
in graphical form in Figs. 1 through 6. Since the 
crystal plates were not coated with conducting 
material and were excited in an air gap, the 
measurements were made under “‘open circuit” 
electrical conditions and strictly speaking should 
be denoted by c?. Because of the small piezoelec- 
tric coupling in beta-quartz, the difference be- 
tween the “‘open circuit” (constant charge) con- 
stants c? and the “short circuit’’ (constant field) 
constants c¥ is less than experimental error so the 
superscripts are omitted in this paper. 

During the measurement of the frequencies of 
vibration for each mode no attempt was made to 
set the furnace at the temperatures for which 
values of the constants are shown. Instead, the 
furnace was permitted to come to equilibrium at 
a series of arbitrary temperatures. Interpolation 
on a graph of these frequencies permitted the 
selection of values for use in computing the K? 
terms for each mode at the same temperature. 

It is possible to compute the elastic moduli 
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Fic. 11. The sy; adiabatic elastic modulus for beta-quartz 
as a function of temperature. 


from the elastic constants, as shown, for example, 
by Wooster.!® The values for the moduli obtained 
in this way are shown in Figs. 7 through 11. 

Perrier and de Mandrot" measured by iso- 
thermal methods the s1; and s33 moduli as well as 
a modulus at 50° with respect to the optic axis for 
beta-quartz. Their values for s;; and s33 are 
plotted without correction to adiabatic conditions 
beside the adiabatic values in Fig. 7 and Fig. 8 
respectively. 

The elastic constant c;2, shown in Fig. 4, be- 
comes negative near the transition temperature. 
In this connection reference is made to the data 
in Table V of Atanasoff and Hart.* They have 
tabulated data for ¢1, C66, C44, and C14 for alpha- 
quartz up to the transition temperature. Using 
Eq. (10), ci2 can be computed from their values 
given for the constants ¢,; and cg. Table II was 
obtained in this manner and shows the way ¢12 
changes with temperature. Apparently the elastic 
constant, Cis, is negative during the transition 
and changes sign in the beta-state at a tempera- 
ture a little above 580°C. It is also of interest to 
note that all of the elastic constants for beta- 
quartz become stiffer as the temperature in- 
creases. 

Lawson” has measured the adiabatic Young’s 
modulus for beta-quartz along a direction at 45° 
to the optic axis. This was done by setting up 
longitudinal vibrations in a rod-shaped specimen. 
The modulus was computed from the relation : 


1/33’ =4f?L2p(1 +k), (12) 


1W. A. Wooster, A Text Book on Crystal Physics 
(Cambridge University Press, Teddington, England, 1938), 
pp. 235-237. 

11 A, Perrier and B. de Mandrot, Memoires de la Société 
Vaudoise des Sciences Naturelles (Imprimeries Réunies 
S. A., Lausanne, 1923), No. 7, pp. 333-363. 

122 A. W. Lawson, Phys. Rev. 59, 608 (1941). 
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TABLE II. Values for cys as nd Harts from the data for 
Cu and ces of Atanasoff and Hart (x= en - 240. ) 


cir X10 cu X10” 

Temp. °C dynes /cm? Temp. °C _dy nes /cm? 
0 6.70 500 - 13.8 
100 4.90 540 —19.3 
200 2.60 550 —21.8 
300 —0.70 562 —25.8 
400 —5.40 573 — 32.9 


where f is the frequency, L the length of the rod, 
p the density, and & is a Rayleigh correction to 
compensate for the depression in frequency due 
to the loading effect of the finite rod cross section. 
The magnitude of the correction term k for 600°C 
was estimated to be 0.003. The density for quartz 
at 600°C used by Lawson was 2.518 g/cm‘*. If this 
be changed to 2.533 g/cm* to permit comparison 
with the data being presented here, Lawson’s 
determination of the s3;’ modulus is reduced from 
1.067 x 10-"* cm?/dyne to 1.061 « 10-" cm?/dyne. 

The s33’ adiabatic modulus at 45° to the optic 
axis can be computed from the elastic moduli,” 
shown in Figs. (7) to (11), by 
relation : 


means of the 


ae 2\2 > 
S33 = S11(1 — a3")? + 533034 


+ (2513 +544)(1—a;3*)a;*. (13) 


Since the direction cosine, a3, is equal to v2/2 for 
the 45° orientation, expression (13) reduces to: 


$33 = 4 (Sir +S33+S44+2513). (14) 


1 W. Voigt, reference 5, p. 740. 
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The values to be used for the moduli on the right- 
hand side of Eq. (14) at 600°C are: 


si1= 0.941 ’ 
533> 1.062, 


54g = 2.773, 
$13= —(.262, 


each multiplied by the factor 10~-" cm?/dyne. 
When these constants are substituted into Eq. 
(14), S33’ becomes 1.062 10-" cm?/dyne, to be 
compared with Lawson’s modified result af 
1.061 X 10-" cm?/dyne. At higher temperatures, 
the 45° modulus, as computed from the elastic 
moduli, deviates in a uniform manner from the 
45° modulus as computed from Lawson's data, 
the former being about 5 percent lower at 800°C. 
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Dark-Field Electron Microscopy 


II. Studies of colloidal carbon 


C. E. HALy 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received September 8, 1947) 


Several commercial carbon blacks were examined with dark-field illumination in the electron 
microscope. Anomalous intensities in the dark-field images are interpreted as resulting from 
coherent scattering by parallel layer groups of graphitic planes, which were previously identified 
by Biscoe and Warren using x-ray diffraction. The dark-field micrographs indicate that in a non- 
reinforcing and a semireinforcing black, the parallel layer groups near the surface of the roughly 
spherical particles are oriented with their representative vector approximately normal to the 
surface. No definite evidence was obtainable for such an arrangement in a reinforcing channel 
black with average particle diameter about 300A. After prolonged heat treatment to the point 
where graphitization occurs the channel black exhibits numerous diffraction images in dark- 
field, indicating the presence of crystals with dimensions in the range 20 to 90A. 











INTRODUCTION 


ARBON black produced by incomplete com- 
bustion of gaseous hydrocarbons is widely 
used as a pigment in inks and lacquers and is of 
particular importance as an ingredient in rubber 
compounding. There are many grades, differing 
from one another in ultimate particle size and 
other properties. Average particle size may be 
from 50A up to several thousand A. Biscoe and 
Warren! have shown through x-ray analysis that 
carbon black is not finely divided graphite. The 
structure is one of true graphite layers arranged 
roughly parallel and equidistant, but otherwise 
completely random as to orientation about an 
axis normal to the layers. The interplanar spacing 
is somewhat larger than that for graphite. Di- 
mensions of the parallel layer groups in different 
blacks were found to be from 12.1 to 26.6A for 
the axis perpendicular to the graphite planes 
(L.) and from 18.6 to 47.5A parallel to the planes 
(L,). These results indicate that a carbon particle 
is a mesomorphic solid, neither crystalline nor 
amorphous, for which class of substances they 
suggest the term “‘turbostratic.’”’ When carbon 
black is heated to about 1000°C and higher, the 
graphite layers grow and thermal agitation allows 
orientation of the graphite planes over limited 
regions to produce small graphite crystals. At 
the same time there is a decrease in the interplanar 


1 J. Biscoe and B. E. Warren, J. App. Phys. 13, 364 
(1942). 
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spacing. From measurement of the x-ray line 
breadths, Biscoe and Warren have estimated the 
average dimensions of these graphite crystallites 
to be up to 40.0A for ZL, and up to 65.2A for La. 

Further evidence concerning the structure of 
carbon black comes from the calorimetric de- 
termination of heats of adsorption reported by 
Beebe, Biscoe, Smith, and Wendell.? Experiments 
on the same black with and without devolatiliza- 
tion, indicate that it is the state of the underlying 
solid rather than the presence of a chemisorbed 
layer which determines the magnitude of the 
binding energy for the adsorbed molecules of 
nitrogen. The binding energy for a monolayer is 
considerably higher than the heat of liquefaction. 
In the initial stages of adsorption, the binding 
energy is abnormally high, indicating that the 
first fraction of a monolayer is bound at active 
centers. When the black is heat-treated to the 
point where graphitization occurs as in the x-ray 
experiments, there is a marked decrease in the 
heat of adsorption in the initial stage, indicating 
disappearance of the active centers as a more 
ordered structure is produced. Although the role 
of the active centers in rubber reinforcement 
remains to be determined, the surface activity as 
determined by heat measurements closely paral- 
lels the rubber-reinforcing properties of the black.’ 


2 R. A. Beebe, J. Biscoe, W. R. Smith, and C. B. Wendell, 
J. Am. Chem. Soc. 69, 95 (1947). 

3 W. R. Smith, F. S. Thornhill, and R. I. Bray, Ind. Eng. 
Chem. 33, 1304 (1941). 
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The electron microscope is widely used to 
measure the particle size of carbon blacks and to 
observe their state of aggregation. The literature 
contains many papers on the subject of which the 
work of Ladd and Wiegand‘ and Watson® may be 
cited as examples. Although the ultimate particles 
are in most cases roughly spherical, Hofmann and 











Fic. 1. Dark-field images with tilted beam of carbon 
particles (P33), objective current decreasing from (a) to (c). 
Arrow indicates direction of aperture. 


*W. A. Ladd and W. B. Wiegand, Rubber Age 57, 299 
(1945). 
5 J. H. L. Watson, J. App. Phys. 17, 121 (1946). 
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Héper® have shown that after prolonged heat 
treatment at 3000°C the particles develop angu- 
larities in outline. In a recent note’ the author 
included a dark-field micrograph in which images 
of carbon particles exhibited regions of enhanced 
intensity at their edges. The reason for this ap- 
pearance was not obvious, but it was suggested 
that the enhanced intensities were in some way 
related to the crystal structure near the particle 
surface. It appeared that the dark-field method 
was capable of showing structural features that 
are not visible in bright-field. These phenomena 
have been studied further with techniques de- 
scribed in Part I and are the subject of the 
present report. 


OBSERVATIONS OF CARBON PARTICLES IN 
DARK-FIELD 


The bright spots which appear at the edges of 
dark-field images of certain relatively large car- 
bon particles,’ exhibit peculiarities in the manner 
of their occurrence and in their appearance in a 
through-focus sequence. They occur only on 
edges directly toward or directly away from the 
aperture (the region of bright-field illumination) 
and never at edges whose tangent is in the aper- 
ture direction. In some instances, enhanced spots 
have been observed apparently in this latter 
position, but further study has established that 
these come from smaller particles, sometimes 
poorly resolved, at the edge of a larger one. Spots 
away from the aperture which are produced by 
pencils bent under the particle, are not in focus at 
the same objective setting as those toward the 
aperture, which are produced by pencils reflected 
away from the particle. This anomaly is shown in 
the three micrographs reproduced in Fig. 1. The 
material, consisting of relatively large particles, 
is known as P33 and is classified as non-reinforcing 
in rubber. The objective current was decreased in 
steps from (a) to (c) and the aperture direction is 
indicated by an arrow. In (a) the spots are 
visible on edges away from the aperture, but the 
edges toward the aperture appear as a broad, 
bright line. As the objective current is decreased 
from (a) to (c) the images away from the aperture 
go out of focus and those toward the aperture 
come into focus. For these micrographs, the 


*U. Hofmann and W. Héper, Naturwiss. 32, 225 (1944). 
7C. E. Hall, J. App. Phys. 18, 588 (1947). 
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incident beam was tilted with respect to the lens 
axis to minimize the displacement of the en- 
hanced images that normally occurs when the 
objective current is altered. The fact that the 
two sets of images are not in focus at the same 
objective setting indicates either that observed 
structures are at different levels, or that the re- 
spective imaging pencils make quite different 
angles with the lens axis. 

In Fig. 2 is reproduced a typical micrograph 
for which the objective current was decreased 
slightly below the value for best focus so that the 
object plane is above the object. In this case the 
incident beam is closely parallel to the lens axis. 
The scattered intensity toward the aperture is 
displaced into the particle, while the enhanced 
spots sway from the aperture are displaced away 
from the particles. The displacements indicate 
that the scattered pencils make relatively large 
angles with the lens axis, and the coherency of the 
displaced images indicates that the pencils are 
confined to narrow relative aperture. These are 
just the characteristics of pencils produced by 
Bragg reflections. It is therefore concluded that 
the enhanced intensities at the particle edges are 
owing to the presence of crystalline structure. 
However, the enhanced images have diameters 
in the range 100 to 900A and the images are 
quite diffuse compared to dark-field images of 
single crystals in this range of dimensions. This 
particular carbon (black was studied by Biscoe 
and Warren! who give the values, L.= 16.8 and 
La=28.0A. Obviously £ single crystals of these 
dimensions could not account for images of the 
size and intensity observed at the edges. The 
diffraction images must consist of the super- 
imposed images of assemblies of parallel layer 
groups. The fact that the diffraction images are 
deflected approximately as a unit indicates that 
within an assembly the individual crystallites 
are oriented with their reflecting planes nearly 
parallel to one another. 

That enhanced intensities never occur on edges 
whose tangent is the aperture direction can mean 
one of two possibilities, (a) diffracted rays are 
prevented from leaving the particles with an 
appreciable component tangential to the image 
periphery by refraction or some other factor, or 
(b) reflecting planes responsible for the principal 
scattered intensity do not occur near the surface 
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Fic. 2. Dark-field images of carbon particles (P33), 
objective current lower than for best focus. Arrow indicates 
direction of aperture. 


except with their representative vector approxi- 
mately normal to the surface. Possibility (a) is 
discounted because there appears to be no reason 
why, if the particles contain a random assembly of 
crystallites, rays should not leave at all azimuths. 
Possibility (b) is therefore proposed as the 
explanation. To identify the planes most prob- 
ably responsible for the diffraction images, we 
refer to the x-ray diagram of Biscoe and Warren.! 
Their Fig. 1(a) shows that almost all the diffracted 
intensity is contained in the rather broad (002) 
ring. Since relative electron diffraction intensities 
follow those in x-ray diagrams rather closely, it is 
assumed that the diffraction spots in the dark- 
field electron microscope images result mainly 
from the (002) spacing. This would mean that the 
(002) vector is, near the surface, always approxi- 
mately normal to the surface. The conclusion 
then follows that in carbon particles of this 
particular sample the parallel layer groups origi- 
nally determined from x-ray analysis are, near 
the surface, all oriented with the graphite layers 
approximately parallel to the surface of the 
particles, but are otherwise randomly disposed 
with respect to one another so that there is lack 
of coherency between the individual groups. 

The P33 carbon black shown in Figs. 1 and 2 
has an average particle diameter of about 1400A 
and is classified as non-reinforcing in rubber. 
Several other types of black were examined, one 
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of which is classified as semireinforcing and is 
known under the trade name of Sterling R. The 
average diameter in this case is about 800A. 
Dark-field micrographs exhibit phenomena simi- 
lar to those described above and typical micro- 
graphs of a group of particles are reproduced in 
Fig. 3. For (a) the objective current is at, or a 
little higher than, that for best focus and the 
scattered intensity from edges away from the 
aperture is imaged as diffuse spots which taper 
inward. In (b) the intensity from edges away 
from the aperture is deflected away from the 
particles to form brush-like streaks, while the 
intensity on edges toward the aperture appears as 
spots. (Different exposure times were used for the 
upper and lower parts of the print in (b) in order 
to show the intensity distribution to better 
advantage.) This effect is like that shown in 
Fig. 2, but the resolution of the deflected in- 








Fic. 3. Dark-field images of carbon particles (Sterling R), 
(a) objective current close to best focus, (b) objective 
current lower than for best focus. In (b), upper and lower 
portions of the print have been given different exposure 
times. Arrow indicates Girection of aperture. 
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tensity into narrow streaks indicates inhomo- 
geneities in the surface of 100A and less. That fine 
structure was not observed in the deflected in- 
tensity as shown in Fig. 2, may mean either that 
the structure for P33 is too small to be resolved, 
or that there is too much overlapping in the 
image. 

A few comments should be made on the curious 
alterations that take place in the image with 
alterations in objective current as illustrated in 
Figs. 1, 2, and 3. It is to be noted that the 
distribution of intensity seen in the image is a 
cross section, real or virtual, of the scattered 
beam at a plane determined by the objective 
current. There is a notable difference between the 
two sets of images, those on edges toward the 
aperture and those at edges away from the 
aperture. In an image which is in best focus, 
edges away from the aperture exhibit diffuse 
spots which taper inward, while the structures at 
edges toward the aperture tend to fuse together 
along the periphery as seen in Figs. 1(a) and 3(a). 
When the objective current is decreased from 
best focus, the typical behavior is as shown in 
Figs. 2 and 3(b). (Figure 1 is not typical of this 
condition since the incident beam was tilted.) At 
edges away from the aperture the scattered in- 
tensity fans out normal to the periphery in fine 
streaks which may or may not be resolved and at 
edges toward the aperture the intensity occurs as 
diffuse spots which are in the nature of a virtual 
crossover above the object plane (Fig. 3(b)). 
These differences indicate that the scattering 
mechanism is quite different at the opposite 
edges. Although it is not possible to give a 
quantitative explanation of these phenomena, 
some of the factors involved may be mentioned. 
If, ashas been concluded, the parallel layer groups 
near the surface are oriented with their graphite 
planes approximately parallel to the surface, the 
scattered beam is reflected at edges away from 
the aperture as if from a concave surface, while 
the scattered beam at edges toward the aperture 
is reflected as from a convex surface. Also, in 
considering the possible factors contributing to 
the difference in appearance and behavior of the 
two sets of images, the effect of refraction at the 
particle surface must be considered. As was 
pointed out in Part I, when rays are incident at 
small angles to a surface, the refraction deflection 


JOURNAL OF APPLIED PHYSICS 




















Fic. 4. Reinforcing channel black (Spheron Grade 6), (a) bright-field, (b) dark-field. 


may be of the same order of magnitude as that 
given by the simple Bragg law. Furthermore, 
refraction effects would be quite different ac- 
cording to whether the scattered pencils were 
bent under the particle or reflected away from 
the particle. In view of these differences in the 
scattering mechanism, it is not surprising that 
images at opposite edges appear different at a 
given objective setting, or that the cross sections 
of the scattered undergo anomalous 
changes in intensity and shape with variations in 
objective current. Qualitatively, rays can be 
traced through particles near the edges to ac- 
count for the observed anomalies in the images, 
but without more specific information on the 
structure of the particles or the magnitude of 
refraction effects, a rigorous treatment of the 
problem is not possible at present. 

Probably the most interesting carbon blacks 
are those classified as reinforcing for rubber. 
Particle diameters are about 300A for this type. 
Unfortunately, the particles occur in aggregates 
which are difficult to separate and this, together 
with the small particle size, makes them difficult 
to study in dark-field. An example, shown in 
Fig. 4, is known commercially as Spheron Grade 
6. In the bright-field micrograph may be seen the 
usual clusters of particles about 300A in diame- 
ter. In dark-field the individual particles are 
poorly resolved. Superimposed on the images of 
the clumps are weak spots, difficult to observe 


beams 
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because of the diffuse intensity, but which are 
probably to be associated with assemblies of the 
parallel layer groups because they exhibit slight 
displacement with changes in objective current. 
Evidence, however, is lacking for any high degree 
of order as has been described for the larger 
blacks. As far as the electron microscope can 
show at present, the parallel layer groups are too 
small to be resolved individually and are randomly 
disposed. The x-ray results for this material are, 
L.=12.7 and L,=20.0A. 

A sample of Spheron 6 which had been 
graphitized in commercial production was ex- 
amined in dark-field in order to determine whether 
the changes observed in the x-ray work! could be 
detected with the electron microscope. Micro- 
graphs of the heat-treated sample (known com- 
mercially as Graphon) are shown in Figi 5. The 
material had been heated in commercial pro- 
duction for about 24 hours at about 2500°C. In 
bright-field, (a), the particles appear to be 
sintered together and some show internal struc- 
ture. In dark-field, (b), the aggregates exhibit 
numerous, sharp diffraction images with mini- 
mum dimensions of about 50A. The relation 
between image diameter and crystal dimensions 
is presumably like that estimated in Part I, 
Fig. 14. The largest diffraction images in Fig. 5 
are about 90A and where they are intense, they 
probably represent crystals of about this size. 
The weak images probably represent crystals of 
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the minimum detectable dimensions which were 
estimated in Part | to be about 20A. Although 
there is some uncertainty in estimating crystal 
dimensions from such small dark-field images, it 
is concluded that the graphite crystals are 
present with dimensions between 20 and 90A. 
Biscoe and Warren give L, =40.0A and L,=65.2A 
for a sample heated at 2800°C for two hours, 
which was the longest heating time used at com- 
parable temperatures. Of the x-ray samples, this 
one most closely approximates the conditions of 
preparation of the commercial sample used in the 
electron microscope work. Since it would be ex- 
pected that crystals would be present on either 
side of the average x-ray values, it is apparent 
that the electron microscope observations are 
quite consistent with what would be expected 
from the x-ray results. 

Another carbon black examined in dark-field 
was Standard Black (acetylene 
black). Watson® has studied this material in 


Shawinigan 


bright-field and has observed, attached to the 


spherical particles, rod-like elements which are 
apparently platelets oriented with their faces ap- 
proximately parallel to the incident beam. Similar 
structures have been observed in this work, but 
only with this particular black. The rods are not 
very numerous and the bulk of the material ap- 
pears to consist of roughly spherical particles of 
about 500A average diameter. In bright-field the 





particles exhibit considerable internal structure. 
In dark-field, numerous faint streaks about 75A 
wide, appear superimposed on the images of the 
particle aggregates and project around the edges. 
This structure has been seen only with this 
particular black and seems to indicate the pres- 
ence of single, relatively large parallel layer 
groups. The streaks are like those which fan out 
from the edges of the Sterling black described 
above, but in the Shawinigan black they occur 
around the edges and throughout the aggregates. 
It is concluded that the parallel layer groups are 
relatively large in this black and occur in less 
orderly arrangement near the surface than in P33 
or Sterling. The conclusion that the parallel layer 
groups are unusually large is consistent with 
Watson’s observation of platelets in bright-field 
and it is also consistent with the observations of 
Biscoe and Warren, who calculate that the 
average dimensions of the parallel layer groups in 
this material are almost twice as large as those in 
the other blacks examined in the x-ray study. 


CONCLUSIONS 


The anomalous intensities appearing in dark- 
field micrographs of carbon particles are in- 
terpreted as originating in assemblies of the paral- 
lel layer groups originally determined through 
x-ray analysis. From the assumption that most of 
the coherently scattered intensity is in the (002) 





Fic. 5. The same type of black as shown in Fig. 4 after heat-treatment (Graphon), (a) bright-field, 
(b) dark-field. 
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reflection, it is concluded that for a sample of non- 
reinforcing black and another classified as semi- 
reinforcing, the parallel layer groups are oriented, 
near the particle surface, with the graphite planes 
approximately parallel to the surface. Evidence 
for structure of this type is lacking for a rein- 
forcing channel black of much smaller particle 
size. Whether the absence of such evidence for 
the reinforcing black is owing to lack of resolution 
or contrast, or whether it represents a significant 
difference in the structure, it is not possible at 
present to say. After prolonged heating at ele- 
vated temperatures, a sample of channel black 
exhibited numerous sharp diffraction images indi- 
cating the presence of graphite crystallites with 
estimated dimensions between 20 and 90A. The 
absence of this type of image in the unheated 
samples is confirmation of the mesomorphic 
nature of colloidal carbon as determined from 
x-ray analysis. Observations on a sample of 
acetylene black indicate that in this material the 
parallel layer groups are relatively large. Evi- 
dence for preferential orientation is much less 
marked for this material than for the non- 
reinforcing and semireinforcing samples. 

It is concluded that the size and relative 
orientation of the parallel layer groups vary from 
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SAMPLES. 
Sample Manufacturer 
P33 Thermatomic Carbon Company, New York 
Sterling R Godfrey L. Cabot, Inc., Boston 
Spheron Grade 6 Godfrey L. Cabot, Inc., Boston 
Graphon Godfrey L. Cabot, Inc., Boston 


Shawinigan Standard Shawinigan Products Corporation, New York 








one type of black to another. No correlation has 
been made between the observed structures and 
the properties of carbon blacks in rubber com- 
pounding and other applications, but the results 
are strongly suggestive that some correlations 
may be possible with the acquisition of more and 
improved data. Although the resolution in dark- 
field is generally inferior to that in bright-field, 
the dark-field method appears to be useful in 
providing information concerning the crystalline 
structure of small colloidal particles. 
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Mass Emissivity of Powder Gases in Solid Fuel Rockets’ 


S. S. PENNER* 
Allegany Ballistics Laboratory, Cumberland, Maryland 
(Received October 3, 1947) 


An analysis has been carried out to determine an upper limit for the effective mass emissivity 
of the radiating gases in a solid fuel rocket motor. The calculations depend upon a comparison 
between experimentally observed performance data relating to a rocket motor and theoretical 
results based on a theory for the effect of radiation on the performance of solid fuel rocket 
motors. A simplified treatment is developed for correlating the rate of burning of the powder 
with the temperature rise of a powder grain caused by the absorption of radiant energy. 


INTRODUCTION 


SOLID fuel rocket motor is powered by a 
charge of combustible material. The power 
output of the motor depends upon the rate of 
generation of gaseous combustion products from 
the burning solid, on the temperature of the 
combustion products, the dimensions of the 
rocket chamber, etc. One of the factors which 
determines the rate of gas generation and, there- 
fore, the performance of the rocket motor is the 
temperature of the powder before burning. This, 
in turn, is related to the absorption of radiant 
energy from the hot gases surrounding the solid 
fuel.” 

The surface of a burning powder grain recedes 
in a direction normal to itself. As the grain burns, 
heated gases are evolved which may or may not 
undergo chemical reactions in the gas phase. 
Radiant energy is assumed to be emitted uni- 
formly in all directions from the gaseous envelope 





! This is the first of three articles abstracted from OSRD 
report No. 5251 which was written in 1945 at the Allegany 
Ballistics Laboratory. The work was carried out under 
OSRD Contract OEMsr-273 with the George Washington 
University, Washington, D. C. The original paper was 
reissued as Part IV of OSRD report No. 5817. 

A summary of the original paper is given in the Bibli- 
ography of Scientific and Industrial Reports, Office of 
Technical Service, Department of Commerce, Vol. 3, No. 4, 
Oct. 1946. The PB number is 34764. 

The radiation program at the Allegany Ballistics Labora- 
tory was initiated by the Director of Research, Dr. R. E. 

Gibson. Dr. W. H. Avery, to whom the author is indebted 
for advice throughout the course of the work, directed the 
radiation studies. 

The author wishes to acknowledge the aid of the Edi- 
‘torial Staff of the Jet Propulsion Laboratory, California 
Institute of Technology, in preparing this paper for 
publication. 

* Present address: Jet Propulsion Laboratory, California 
Institute of Technology, Pasadena, California. 

2W. H. Avery, OSRD report No. 3880, July 1944; J. 
Beek, Jr., W. H. Avery, M. J. Dresher, F. T. McClure, 
and S. S. Penner, OSRD report No. 5817, 1946. 
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surrounding the burning powder. 
this energy 


A portion of 
is absorbed by the gases separating 
the emitting layer from the powder surface. The 
light which is not reflected from the powder 
surface is absorbed by the powder grain. The 
decrease in intensity per unit distance depends 
on the absorption coefficients of the powder grain 
and on the wave-length distribution of the inci- 
dent radiation. The burning rate of the powder 
grain, in turn, depends on the temperature rise 
accompanying the absorption of radiation. 

The preceding physical picture requires several 
restrictive assumptions before useful calculations 
can be carried through, viz., 

(1) The gaseous envelope radiates like a gray- 
body, i.e., the emissivity of the gases is the same 
for all wave-lengths. This implies that the tem- 
perature dependence of the radiation intensity is 
the same as that of a blackbody. It is recognized 
that this simplification cannot correspond to fact. 
However, it is necessitated by the absence of 
more detailed information. 

(2) All of the radiant energy striking the 
powder surface is assumed to be incident in a 
direction normal to the surface of the powder 
grain. More general treatments without this sim- 
plifying assumption are quite complex and have 
been carried through in only a few selected 
cases.*: 4 

A complete treatment of the effect of radiation 


*B. L. Crawford, Jr. and R. G. Parr, University of 
Minnesota Report No. 27, February 1945. This memo- 
randum is available in the files of Division 3, Section H, 
NDRC, 

‘J. Beek, Jr., Armor and Ordnance Report No. A-287, 
OSRD No. 4033, Aug. 1944. An abstract of this report 
appears in the Bibliography of Scientific and Industrial 
Reports 2, 883 (1946), PB 34774. 
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on the performance of solid fuel rockets requires 
study of the following phenomena: 

(1) Absorption of radiation by powder grains, 
necessitating determination of absorption coeffi- 
cients as a function of wave-length. The required 
experimental work was carried out by F. Daniels 
and co-workers’ at the University of Wisconsin 
and by B. L. Crawford, Jr., et al.6 at the Uni- 
versity of Minnesota. 

(2) Emission of radiation from the hot gaseous 
products of combustion of a solid fuel rocket 
motor. The radiation emitted from the gaseous 
envelope has been studied experimentally by 
R. S. Craig? who employed an experimental 
technique developed by F. Kracek and W. S. 
Benedict. The present treatment obviates the 
necessity for experimental emissivity determina- 
tions by allowing calculation of an effective emis- 
sivity from observed performance data. 

(3) Calculation of the temperature rise associ- 
ated with the absorption of radiant energy and 
determination of the change in burning rate 
associated with the temperature rise. Avery’s 
theory for calculation of temperature rise and 
its effect on performance has been used in the 
present paper. 


I, AVERAGE TEMPERATURE RISE IN A POWDER 
GRAIN, CAUSED BY THE ABSORPTION 
OF RADIANT ENERGY 


Avery’ has shown how to calculate the average 
temperature rise in a powder grain resulting 
from the absorption of radiant energy. An ex- 
tended treatment, utilizing Avery’s method, is 
given here for a solid cylindrical powder grain 
which is restricted by the use of inert material 
so that it can burn from one end in only one 
direction, namely, parallel to the cylinder axis. 
Application of the theory developed in the follow- 
ing sections to experimentally determined per- 
formance data permits direct evaluation of the 
mass emissivity which makes theoretical calcu- 
lations consistent with experimental observations. 

If « represents the emissivity in the wave- 


°S.S. Penner and F. Daniels in NDRC Report No. 
A-485, OSRD Report No. 6559, Dec. 1945. 

°B. L. Crawford, Jr. and R. Gledhill in OSRD Report 
No. 6374, Dec. 1945. 

7R.S. Craig, OSRD Report No. 5832, 1945. 

* F. C. Kracek and W. S. Benedict, Armor and Ordnance 
Report No. A-252, OSRD No. 3291, 1944. 
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length region between \ and A+d) of the hot 


gases, produced as the result of powder burning, 
then 


log(1— «) = —kygp, (1) 


where k, represents the mass emissivity for the 
corresponding wave-length region, p is the den- 
sity of the powder gases, and g stands for the 
path length in a given direction throughout 
which emitters of radiant energy are uniformly 
distributed between the energy receiving powder 
surface and the rocket chamber wall. 

The wave-length dependence of k) is neglected 
in all of the following calculations in accord with 
the assumption of gray body emission. With this 
simplification Eq. (1) may be written as 


Ih=Io,(1 —e-**2) (2) 


after replacing e by 4,/Jo,, where Jo, represents 
the energy radiated by a blackbody emitter at 
the temperature of the powder gas per sq. cm of 
area per sec. over an angle of 2m steradians in 
the wave-length region under consideration. 

A cylindrical rocket chamber of length L+Lo, 
containing a powder grain of length L, is shown 
in Fig. 1. The powder grain may be divided into 
small sections which are not necessarily of the 
same width: AL;, ALs, ---AL,. The average rate 
of burning of the slab of width AL, is designated 
as r; and corresponds to an average gas density 
pi, etc. Thus, the slab of width AZ; burns with an 
average rate r; at an average gas density p;. 

It is desired to determine the temperature rise 
of the section of width a as the result of the ab- 
sorption of radiant energy. This section is located 
at a distance } from the original surface of the 
powder grain. 

The distance from the right-hand surface of 
the slab of width a to the burning powder surface 
at any time is designated as X. While the powder 
burns over a distance ALi, 
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Fic. 1. Cylindrical motor chamber containing a 
powder grain which burns from one end only. 
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Fic. 2. Progressive pressure vs. time curve. 


while burning occurs over the distance AL», 
X = (b—AL;) — rot; 
while burning occurs over the distance AL ;, 
X =b—(ALi+AL2+--++AL;-1)—rjt, etc. (3) 


Only radiant energy incident in a direction 
normal to the burning powder surface is con- 
sidered, as a first approximation.** With this 
simplification Eq. (2) becomes 


Ih = Io, (1 — e~*?1“2) (4) 


while the powder burns through the distance 
AL;. Equations similar to Eq. (4) hold for the 
incident intensities of radiation while the powder 
burns through other selected distances. 

The decrease in radiation intensity per unit 
length of powder is strongly dependent on the 
wave-length distribution of the incident light 
since the absorption coefficients of the powder 
vary from one wave-length to another. 

The average absorption coefficient of the 
powder for light in the wave-length region be- 
tween A and A+dA is designated as K, in the 
following, and defined by the relation 


DT! = he-*, (5) 


where J,’ is the radiant energy in the wave-length 
region between \ and A+dA transmitted through 
a length x of powder when the incident radiant 
energy is J, for the same wave-length region. 

The preceding physical picture leads to Eq. (6) 
for the radiant energy Ea absorbed by the 
powder grain in the wave-length region between 

\ and A+dd (see Appendix for details). 


** Neglect of radiant energy incident in a direction other 
than the normal to the powder surface is justified by the 
large absorption coefficient of the propellant charge. The 
contribution of radiant energy from directions other than 
the normal to the powder surface is relatively small for a 
strongly absorbing powder at some distance below the 
powder surface. 
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The temperature rise, A7\, produced by the 
absorption of Ea, cal./sq.cm/sec. is given by the 
relation 
AT, = Fay /Cp,a, 


where C represents the heat capacity in cal./ 
degree/gram of the powder, and p, is the powder 
density. When a is very small the approximation 
may be made that 


(1—e-*a*) /a~K,. 
The temperature rise of the powder grain associ- 
ated with the absorption of radiant energy in the 


wave-length region between \ and \+dA is, there- 
fore, given by Eq. (7): 


exp | —K,(o— , AL;)} 
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Fic. 3. Temperature vs. time curve. 


For accurate calculations the temperature rise 
must be computed for each individual wave- 
length interval AN=AX,, for which an average 
value K,=Ky, can"be assigned to the absorption 
coefficient of the powder. It is then necessary to 
sum over m in order to obtain the total tem- 
perature rise. ; 

Equations (6) and (7) require more detailed 
information than is usually available. For this 
reason it is necessary to assume an average burn- 
ing rate at an average pressure for reasonably 
long distances. For each of these distances Eq. (7) 
reduces to Eq. (8): 
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Fic. 4. Radiation temperature rise at various distances for 
K=2 cm", O.5u<A<1.0p, Lo=7.00 cm. 
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where 7 represents the average burning rate in 
cm/sec. and # is the average gas density in 
grams/cc. , 


II. APPROXIMATE METHOD FOR THE CALCU- 
LATION OF PROGRESSIVENESS OF 
PRESSURE-TIME CURVES 


Experimental techniques have been evolved 
which permit the determination of the pressure 
in a rocket chamber as a function of time.® A 
pressure vs. time curve of the type shown in 
Fig. 2 is said to be progressive. One of the factors 
which determines the degree of progressiveness 
is the depth of the hot radiating gas layer which 
increases continuously in a motor chamber of the 
type shown in Fig. 1. 

In order to correlate the effect of radiation 
with progressiveness of the pressure vs. time 
curve, it is necessary to develop a relation be- 
tween instantaneous powder temperature and 
instantaneous burning rate. 

The product 7AT can be calculated from Eq. 
(8) for a definite powder temperature with fixed 
path length and constant pressure. The actual 
temperature rise is obtained by dividing 7AT by 7 
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Fic. 5. Radiation temperature rise at various distances for 
K=10 cm™, 1.0n<A<2y, Lo=7.00 cm. 


®W. H. Avery, R. E. Hunt, and M. N. Donin, OSRD 
report No. 5827, January 1946; W. H. Avery, R. E. Hunt, 
and L. D. Sachs, OSRD report No. 5824, March 1946. 
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Fic. 6. Radiation temperature rise at various distances for 
3cm?<K<60 cm, 2u<A<10pn, Lo= 7.00 cm. 


where an approximate average value is used for 7. 
Dividing b, the initial distance of a from the 
burning powder surface, by 7 gives the time at 
which burning occurs at the right-hand plane of 
the section of width a (cf. Fig. 1). Addition of 
the initial powder temperature to the values of 
AT gives the necessary data for the construction 
of Fig. 3. 

The ordinate of Fig. 3 represents the powder 
temperature at which burning occurs and the 
abscissa the time at which the indicated powder 
temperature is reached. 

Let 7;-1 represent the average temperature 
during the time interval At;_;, and 7; the average 
temperature during the time interval At;. Assum- 
ing the validity of the relation® 


r=c’P"/(T,—T), (9) 
_ it follows that 
dr/r=ndP/P+dT/(T,—-T). (10) 


Here c’, n, and 7, are constants characteristic of 
the powder under investigation, P represents the 
pressure, and T is the powder temperature. Dur- 
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ing the early stages of burning'? 


Spr = CpA P+d(pV) /dt, (11) 


where S=surface area of the burning powder, 
A,=throat area of the rocket nozzle, Cp =dis- 
charge coefficient, V =chamber volume, and r, p, 
and P have their previous meaning. After the 
very rapid initial temperature rise, which occurs 
immediately after ignition of the rocket powder, 
the value of d(pV)/dt is small compared with 
CpyA.P. Therefore, 


r=CpA-P7 Sp, (12) 
or 
dr/r=dP/P. (13) 
Combining Eqs. (10) and (13), 
dP =(P/(1—n) \[d7/(T,-—T) }. (14) 


Let dT ;_;,; represent the average temperature 
rise occurring between the Af;'th and the Af;_,’rst 
second, 1.e., aT 5-1, ;=7;-Ty-1; if dP ;_1, ;=P; 
—P;_, is the corresponding pressure rise, then, 

j g 


from Eq. (14), 
dP j-1, ;=((Pj+Pj-1)/2(1—n)] 





[dT 5-1, (/(T1-3T5— 375-1) J. (15) 
But P;=P;1+dP;-_,,; and, therefore, 
P-1 dT j-1, ; 
1—n T,—3T; $7 -1 
dP; ia (16) 
1 dT j-1, ; 





—_ ; 
2(i—n) T™1—37;—3T5-1 


A pressure vs. time curve can be constructed 
from Eq. (16) provided one value of P= Pp» at a 
definite time =f) is known. Since Eq. (16) is 
based on the assumption that nearly steady state 
conditions exist, an experimental pressure value, 
after the rapid initial pressure rise, is chosen as 
starting point for the calculations. The remaining 
pressure vs. time curve may then be constructed 
by use of Eq. (16). Calculation of 7; requires the 
selection of a definite value for the mass emis- 
sivity of the powder gases. If the assumed value 
for the mass emissivity does not lead to a 
pressure-time curve which agrees with the experi- 
mental curve, a different mass emissivity is 


10 Rocket Fundamentals, OSRD Report No. 3711, May 
1944; PB Report No. 50835, Chapter 3 by R. B. Kershner. 
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chosen, and the calculations are repeated. The 
effective mass emissivity is the value at which 
theoretical and experimental pressure-time curves 
agree with one another. 


III. CALCULATION OF MASS EMISSIVITY BY TRIAL 
AND ERROR 


Calculations of mass emissivity from the pro- 
gressive pressure-time curves observed experi- 
mentally are made on the assumption that all of 
the observed increase in pressure with time is 
caused by radiation. Other factors than increased 
radiant energy incident on the powder surface 
may be responsible for causing progressiveness 
of the pressure-time curve. For example, the 
vaporization of the inert material which prevents 
burning of the powder grain in a direction normal 
to the cylinder axis may be responsible for a 
continuous increase in pressure during burning. 
Nevertheless, a calculation which assumes that 
the only important factor is radiation is of inter- 
est since the result of this study gives an upper 
limit for the mass emissivity. 

Some of the data required for calculation are 
summarized in a series of graphs. The dimension- 
less factor 7ATCp,/Jo, is plotted against initial 
distance from the powder surface for various 
wave-length regions for which average absorption 
coefficients have been assigned to representative 
double base rocket powders*** (see Figs. 4 to 6). 
Experimentally determined absorption coeffi- 
cients for powders which do not contain added 
darkening agents are given for the required wave- 
length regions in each of the figures. Absorption 
coefficient data are taken from the work of 
Daniels,® Crawford,® and their co-workers. 

The term 7AT Cp,/ 0 is effectively independent 
of A, for K,b>10, as can be seen from Eq. (8). 
For this reason the entire wave-length interval 
from 0.50 to 10.0 microns, which is of interest in 
radiation calculations on rockets, can be sum- 
marized for all wave-lengths by the plot shown in 
Fig. 7 for a darkened grain for which 12< K <56. 
In Fig. 7, 6 has been chosen to be greater than 
or equal to 1 cm. 

An experimental pressure vs. time curve is 
shown in Fig. 8 for a double base powder which 





*** A double base rocket powder consists of approxi- 
mately 60 percent nitrocellulose and 40 percent nitro- 
glycerin, compounded with various other ingredients. 
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did not contain added darkening agent such as 
finely divided carbon. The average pressure dur- 
ing burning was P=209 p.s.i. and the average 
rate of burning *=0.447 cm/sec. Also shown in 
Fig. 8 are the calculated pressure-time curves for 
a powder gas temperature T = 2400°K at kp = 0.02 
cm and kp=0.05 cm, and for a powder gas 
temperature T= 3000°K at kp=0.02 cm—. 

It can be seen from Fig. 8 that reasonably good 
agreement exists between experimental and theo- 
retical data for ko=0.05 cm= at a flame tem- 
perature of 2400°K. At this temperature and a 
pressure of 209 p.s.i., p has the value 0.00204 
g/cc. Therefore k=24.5 cm?/g for the given 
conditions. This value of the mass emissivity is 
in good agreement with approximate experi- 
mental findings for powders of the type in- 
vestigated.’ 

The theoretical pressure vs. time curves do not 
flatten out in the same manner as the experi- 
mental curve, as can be seen from Fig. 8. This 
may be caused by the fact that the emitting 
layer of hot gases does not increase uniformly in 
depth during burning. 

The method of determining mass emissivity 
from the progressiveness of the pressure vs. time 
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Fic. 7. Radiation temperature rise at various distances for 
12 cm ?<K<56 cm™, 0.54<4<10p, Lo=7.00 cm. 
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Fic. 8. Experimental and calculated pressure vs. time 
curves_for a restricted rocket¥propellant (Lo=8.8 cm., 
P=209 tbs./in.2, #=0.447 cm/sec., initial powder tem- 
perature = 62°C). 


curve requires knowledge of the flame tempera- 
ture. The calculated mass emissivities do not 
respond to slight variations in flame temperature. 
Nevertheless, it is necessary to use a value for 
the flame temperature within 100 degrees of the 
correct value in order to obtain calculated mass 
emissivities which are of the correct order of 
magnitude. 


APPENDIX TEMPERATURE RISE IN A CYLINDRI- 
CAL POWDER GRAIN BURNING FROM ONE 
END ONLY, CAUSED BY THE ABSORP- 
TION OF RADIANT ENERGY 


The temperature rise in a cylindrical powder 
grain, caused by the absorption of radiant energy, 
may be determined by application of Eggs. 
(3)—(S). 

When the powder burns a distance AL; (see 
Fig. 1) the total radiant energy received by a 
is given by 


tfinal 
Eu=In f [1—exp { —kpi(Lo+rit)} | 
0 


Xexp{—K,X}dt. (A-1) 
The transmitted energy is 
tfinal 
Euy= Im f [1—exp{ —kp:(Lo+rit)} ] 
; xexp{ —K,(X +a) dt. 


The limits of integration in Eqs. (A-1) and (A-2) 
are given by the condition that ‘=0 when X =), 
t=AL;/r; when X=b—AL,. In the section of 
thickness AL;, X =b—r;t. Therefore, the energy 
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(A-2) 


absorbed by the section of thickness a while the 
powder burns through the distance AL; is 
given by 


E, =Io(1—exp{ — Kya} )e~*» 


ALi/ri 
xf [1—exp{ —kpi(Lo+rit)} ] 
0 


Xexp| Kyrit}dt. (A-3) 


Similarly, when the powder burns through the 
section of thickness AL.» the incident energy on 
the section of thickness a is given by the relation 


Eie= Io 


tfinal 
xf [1—exp{ —Rpe[ (Lo+AL1) +ret |} | 
4 


Li/ri 


Xexp{—K,X}dt. (A-4) 


The transmitted energy is therefore 


Ewu=In 


xf [1—exp{ —Rpe[ (Lot+AL;) +ret |} ] 
A 


Li/ri 


Xexp{|—K,(X+a)}dt, (A-5) 


and the absorbed energy becomes 


FE. = Io, (1 —e~*r2) 


x f [1x —kpol (Lot+ALi) +ret ]} J 
Xexp{ —K,X }dt 


subject to the boundary conditions t=AZ;/r: 
when X =b—-AL); t=AL,/r1:+AL2/re when X=b 
—AL,—AlL2, and X=b—ALi—ret at any time 
while the section of thickness AL» is burning. The 
energy absorbed by the section of thickness a 
while the powder burns through the distance AL2 
is, therefore, seen to be 


E.=Io,(1—exp{ —K,a}) exp{ —K,(b—AL,)} 
ALi/rit+AL2/re 
x [1—exp{ —kpe(Lo+ALi+ref)} ] 


ALi/nri 


Xexp|Kyret}dt. (A-6) 


Repeating the preceding argument, it is found 
that the energy E; is absorbed by the section of 
thickness a while the powder burns through the 
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distance AL;: 


E;= Io,(1 —e-**) 


Xexp{ —K,(b—ALi—---—AL;-1)} 
Obi, ,, ALi 
Xx f [1—exp{ —kp,(Lot+ALi 
ATA op Mitt 


+---+AL;1+7et)} | exp| Kyrjt}di 


j-1 
= Io,(1 —e-*)*) exp{| —K,(b—> AL,)} 


i=1 





j-1 
x f [1—exp{ —kefLlet+ 5 AL: +r4)}] 
i=l 


ina (7H 


Xexp|Kyrjt}dt. (A-7) 


In order to obtain the total energy of radiation 
absorbed by a before the flame front reaches a 
it is necessary to sum over all of the j’s from 
j=1 to j=n. The result is 


Ea = > E;= Io, (1 —e-*y) 


7=1 


x exp | ~Ky(b-¥ AL;)} 


j=1 i=l 


ALi 


1 7 


j-l1 
Xx f [1—exp{ —kp;(Lo+ > ALi+r7t)} J 


i=l 


iM~ 


3-1 
ALi 

-— 

Pr a! 


x exp { Kyr3t } dl. (A-8) 


Equation (A-8) reduces to Eq. (6) after in- 
tegration. 





Viscous Slip Along Grain Boundaries and Diffusion of Zinc in Alpha-Brass* 


T’1inc-Sur K& 
Institute for the Study of Metals, The University of Chicago, Chicago, Illinois 
(Received October 6, 1947) 


The viscous behavior of the grain boundaries in 70-30 alpha-brass was demonstrated by 
anelastic measurements. It was shown that the grain boundaries cannot sustain permanently a 
shear stress and have a coefficient of viscosity decreasing with an increase of temperature. 
Using the same alpha-brass specimen, the activation energy associated with the stress-induced 
preferential orientation of pairs of zinc atoms in alpha-brass (which is a diffusion process) was 
determined. This activation energy was found to be identical, within experimental error, to 
the activation energy associated with the grain boundary slip in alpha-brass. This finding indi- 
cates that the grain boundary slip in alpha-brass is a diffusion process and the diffusion mecha- 
nism is similar to the volume diffusion in alpha-brass. 


I. INTRODUCTION 


HE mechanical behavior of grain bound- 

aries in metals has previously been studied 
by considering the grain boundary as an entity 
without inquiring into its true nature.’* Thus it 
has been demonstrated that the grain boundary 
behaves in a viscous manner in the sense that it 
cannot sustain permanently a shear stress and 
has a coefficient of viscosity decreasing with an 








* This research has been supported by ONR (Contract 
No. N-6ori-20-IV). 

'T. S. Ké, Phys. Rev. 71, 533 (1947). 

2 T. S. Ké, Phys. Rev. 72, 41 (1947). 
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increase of temperature.! The concept of a vis- 
cous grain boundary was further elucidated by 
later experiments? where it was shown that when 
a periodic stress is applied to a polycrystalline 
metal, the manifestations of the grain boundary 
viscosity can be expressed as functions of the 
grain size (G.S.) the frequency of vibration, f, 
and the absolute temperature, T, only through 
the parameter (G.S.) Xf Xexp(H/RT), where R 
is the gas constant and H the heat of activation. 

The viscous behavior of grain boundaries 
demonstrated above gives no definite information 
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regarding the nature or the structure of the 
boundary region between adjacent crystals. It 
should be emphasized that the concept of a 
viscous grain boundary is consistent with both 
the two principal theories regarding the nature of 
grain boundaries, e.g., the amorphous cement or 
irregular intercrystal layer theory and the transi- 
tional lattice theory. It can be seen that any 
region of disturbed crystallinity, when considered 
as an entity, may behave in a viscous manner. So 
far, most of the work was done on high purity 
aluminum. The purpose of this paper is to 
describe some findings on alpha-brass, which, in 
addition to the further demonstration of the 
viscous behavior of grain boundaries in metals, 
might shed light on the mechanism of the 
viscous slip along grain boundaries. 


Il. GRAIN BOUNDARY VISCOSITY IN 
ALPHA-BRASS 


If the grain boundaries behave as if they were 
viscous and have a coefficient of viscosity de- 
creasing with an increase of temperature, then as 
the temperature of a polycrystalline specimen is 
raised, the internal friction associated with the 
viscous slip along grain boundaries (a convenient 
measure of this internal friction may be taken as 
the product of the viscous slip along the bound- 
ary and the resistance to slip) must go over a 
maximum value. The reason for this can be 
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Fic. 1. Variation of internal friction with temperature in 
polycrystalline and single-crystal alpha-brass. (Frequency 
of vibration is about half a cycle per second.) 
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illustrated as follows. At low temperatures, the 
viscous slip along the grain boundaries is negli- 
gibly small, and therefore the internal friction is 
small. At high temperatures, the resistance to 
slip along the grain boundaries is negligibly 
small, so again the internal friction is small. Only 
in an intermediate temperature range, when both 
the viscous slip and the resistance to slip are 
appreciable, will the internal friction be appreci- 
able. Such a maximum has been observed in the 
case of aluminum and magnesium.' For alpha- 
brass, internal friction measurements have been 
taken by Zener and his associates from room 
temperature up to 450°C at a frequency of 
transverse vibration of between 600 and 900 
cycles per second.’ It has been found that the 
internal friction increases with the grain size of 
the specimen. However, a maximum was never 
attained. According to the proposition that an 
increase of the frequency of vibration will shift 
the internal friction curve (versus temperature) to 
higher temperatures,” apparently still higher 
temperatures were needed to reach the maximum 
unless a lower frequency of vibration is used. 

In the present experiment, the internal friction 
(logarithmic decrement divided by 7) and the 
shear modulus were determined using a simple 
torsional pendulum. The frequency of vibration 
was about one cycle per second so that the 
measurements could be made by visual observa- 
tion. The apparatus used was practically identical 
with that described before.' The alpha-brass used 
was similar to that used by Burke and Shiau in 
their study on the grain growth in alpha-brass.+ 
It was prepared by melting cathode copper and 
commercially pure zinc together in an induction 
heated graphite crucible. The copper content of 
the final alloy was 70.14 percent, and the per- 
centages of the impurities as determined by 
spectrographic analysis were: Fe, 0.005; Pb, 
0.002 ; Ag, 0.001 ; Cd, 0.001 ; As, 0.005 ; Sb, 0.002. 
A three-inch piece of this specimen, having a 
diameter of one-quarter inch was cold-swaged 
and drawn to a final diameter of 0.035 inch with 
several intermediate anneals at 650°C. The final 
cold-work was about 70 percent reduction in area 


3C. Zener, D. van Winkle, and H. Nielson, Trans. 
A.I.M.E. 147, 96 (1942). 

4]. E. Burke and Y. G. Shiau, Metals Tech. 14, No. 6 
(September, 1947). 
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(70 percent R.A.). A wire 12 inches long was 
mounted up in the furnace and all the measure- 
ments on internal friction and shear modulus 
were taken in an atmosphere of argon having a 
purity of 99.8 percent. The temperature of the 
wire specimen was controlled and measured 
within +1°C with a Tag Celectray indicating 
temperature controller. The torsional stress ap- 
plied to the specimen was very small and the 
maximum shearing strain on the surface of the 
wire was of the order of magnitude of 10~-5. Under 
such a small stress, the internal friction and shear 
modulus were found to be independent of stress. 


A. Variation of Internal Friction 
with Temperature 


The alpha-brass wire described above was 
annealed at 600°C for five hours in an argon 
atmosphere, and it has been found that the 
recrystallized grains have an average grain diame- 
ter of about 0.03 cm. Measurements on internal 
friction and shear modulus were taken from room 
temperature up to 525°C and then coming down 
again to room temperature. The variation of 
internal friction and shear modulus with temper- 
ature are shown, respectively, in Fig. 1 and Fig. 4. 
It can be seen from Fig. 1 that the expected 
internal friction maximum occurs at a tempera- 
ture of 415°C and at this temperature the natural 
frequency of vibration is 0.482 cycle per second. 
In order to determine the heat of activation 
associated with the grain boundary slip, two 
different frequencies of vibration were used and 
the corresponding internal friction curves were 
plotted against the reciprocal of absolute temper- 
ature (1/T) in Fig. 2. It can be shown that the 
internal friction obtained in measurements with 
different frequencies of vibration should be 
identical if the frequency, f, and the temperature 
of measurement, 7, be related through the heat 
of activation, H, according to the relation 


d \Inf/d(1/T) = —H7/R, 
or 


H=R W(f2/f1)/(1/T1—-1/72), 


when two frequencies, fz and f;, are involved. 
The horizontal shift in 1/7 required to bring the 
two internal friction curves in Fig. 2 to superpose 
on each other is 0.055 X 10-*. As the average ratio 
of these two frequencies in the temperature range 
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Fic. 2. Effect of frequency of vibration on internal fric- 
tion associated with grain boundary slip in alpha-brass. 
Frequency of vibration: ©, about 0.5 c.p.s.; X, about 1.5 
c.p.s. (see Fig. 3). 


studied is, from Fig. 3, 3.10, this gives a heat of 
activation of 17 =41,000 calories per mole, which 
should be accurate to within five percent. 

It can be seen from Fig. 1 that this internal 
friction peak does not occur in alpha-brass single 
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polycrystalline alpha-brass. 





yn oN 
o - 
i 


a oe oe oe | 











on 


(FREQUENCY of VIBRATION)" 
2 


n n i L 
100 200 300 400 500 
TEMPERATURE of MEASUREMENT (°C) 


Fic. 5. Variation of shear modulus with temperature in 
alpha-brass crystal. 


crystal. This indicates that the internal fric- 
tion observed can be ascribed only to the grain 
boundaries in alpha-brass. However, it is to be 
noticed that there is a smaller internal friction 
peak around 300°C in both polycrystalline and 
single-crystal alpha-brass. These will be dis- 
cussed in some detail in the next section. 

From the data given in Figs. 1 and 3, we can 
estimate the thermal variation of the viscosity 
coefficient associated with the grain boundaries 
in alpha-brass. It has previously been shown that 
the coefficient of viscosity can be estimated from 
the formula! 


nr = (Gu) rrrd/(GS.), 


where Gy is the unrelaxed shear modulus, r the 
relaxation time, d the effective thickness of the 
grain boundary, and (G.S.) the average grain 
diameter of the specimen. Following similar pro- 
‘cedures as in the case of aluminum, we get 


rr =8.20X 10-" exp(20,500/T), 


utilizing the facts that the internal friction is 
- maximum at 688°K when the frequency of vibra- 
tion is 0.482 cycle per second and the heat of 
activation is 41,000 calories per mole. 

The unrelaxed shear modulus is of the order of 
magnitude of 10" dynes per cm? for all metals 
and d may be taken as of the order of magnitude 
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of 10-§ cm. Since the average grain diameter for 
the alpha-brass specimen is about 0.03 cm, we 
have 


nr» 10° exp(20,500/T7). 


At the melting point of alpha-brass, which is 
940°C, the coefficient of viscosity is, according to 
the above relation, of the order of magnitude of 
10-* poise, which is of the same order as the 
coefficient of viscosity of most molten metals. 


B. Temperature Variation of Shear Modulus 


When an over-all stress, however small, is 
applied to a polycrystalline specimen, the shear 
stress across all grain boundaries, because of 
their viscous behavior, will gradually relax. The 
locking effect of the grain edges and corners will 
insure that the over-all stress relaxation will be of 
limited extent for a fixed over-all strain. The 
theoretical maximum shear stress relaxation in 
metals is,® 


1 —Gr/Gu =1—2(7+50)/5(7 —4e), 


where Gz is the shear modulus of a polycrystalline 
specimen in the case where the grain boundaries 
are viscous and the relaxation of shearing stress 
across the grain boundaries has been complete, 
Gy the shear modulus where no slip occurs across 
the grain boundaries, and o the Poisson’s ratio. 

The temperature variation of the shear modulus 
(which is proportional to the square of the 
natural frequency of vibration) in polycrystalline 
and single-crystal alpha-brass are, respectively, 
shown in Figs. 4 and 5. The straight line in Fig. 4 
was converted from Fig. 5. Let G(T) be the shear 
modulus in polycrystalline alpha-brass at tem- 
perature T, then we have 


G(T)/Gu = (fp/f.)*, . 


where f, and f, are, respectively, the frequency of 
vibration in polycrystalline and single-crystal 
alpha-brass at the temperature 7. The ratio 
G(T)/Gu is plotted as a function of the tempera- 
ture of measurement in Fig. 6. It can be seen that 
the grain boundary relaxation has not been com- 
pleted at the highest temperature of measurement 
but we can infer that Gr/Gy is smaller than but 
close to 0.73. The over-all stress relaxation across 


5 C. Zener, Phys. Rev. 60, 906 (1941) ; see also reference 1. 
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the grain boundaries in the alpha-brass specimen 
is thus greater than, but close to 27 percent. This 
lies within the range of theoretical value of from 
24 to 45 percent when o covers the range from 3 
to } for metals. 


Ill. ANELASTIC RELAXATION IN ALPHA-BRASS 
ARISING FROM ATOMIC DIFFUSION 


It has been mentioned above that there is a 
smaller internal friction peak around 300°C in 
both polycrystalline and single-crystal alpha- 
brass (Fig. 1). The frequency of vibration used 
was of the order of one cycle per second. It is 
believed that this internal friction peak has the 
same physical origin as that observed by Zener in 
70-30 alpha-brass crystal. The optimum tem- 
perature observed by Zener was around 420°C 
using a frequency of vibration of 620 cycles 
per second. The origin of this internal friction 
has been recently determined by Zener as the 
stress-induced preferential orientation of pairs 
of zinc atoms in alpha-brass.? This mechanism 
is formally similar to that arising from the 
stress-induced preferential distribution of the 
interstitial solute atoms in alpha-iron (solute 
atoms: carbon or nitrogen)* and in tantalum 
(solute atoms: carbon or oxygen).°® 

In a copper crystal which has a face-centered 
cubic structure, there are twelve closest neighbors 
for each individual atom and these neighbors are 
so arranged as to form a cubic symmetry. When 
zinc atoms are dissolved in copper to form a 
substitutional solid solution or alpha-brass, an 
elastic distortion of the lattice will take place 
since the zinc atom is larger than the copper atom 
it replaced. If another zinc atom happens to be- 
come one of the twelve closest neighbors of a zinc 
atom, then distortion of the lattice will no longer 
have cubic symmetry and the lattice will be 
extended more along the axis of the pair of zinc 
atoms than along any other directions. In a 
stress-free alpha-brass specimen the pair axes will 
be distributed at random along various per- 
missible crystallographic directions. But when a 
tensile stress is applied to the specimen, the zinc 
pairs will tend to be oriented preferentially along 





6 C. Zener, Trans. A.I.M.E. 152, 122 (1943). 

7C. Zener, Phys. Rev. 71, 34 (1947). 

8 J]. Snoek, Physica 8, 711 (1941). 

* T. S. K@, communicated to the Physical Review. 
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those permissible axes which make the least angle 
with the direction of the applied stress. Some time 
is required for the establishment of this prefer- 
ential orientation of the zinc pairs and, conse- 
quently, strain lags behind the applied stress. 
This will cause all kinds of anelastic effects 
including internal friction. 

It can be seen that if the above picture is cor- 
rect, then the activation energy associated with 
the observed internal friction should be identical 
with the activation energy for the diffusion of 
zinc atoms in alpha-brass. An estimation of this 
activation energy has been made by Zener* from 
the slope of the internal friction curve plotted 
against 1/T and he obtained a value of 33,600 
calories per mole, which is much smaller than the 
value of 41,700 calories per mole for the diffusion 
of zinc (29.08 percent) in alpha-brass.’® A careful 
determination of the activation energy using the 
standard procedure of measuring internal friction 
with various frequencies of vibration is thus 
desirable. 

The alpha-brass single crystals used in the 
present experiment were kindly prepared by Mr. 
Harold Hirsch of the metallurgy section of this 
institute. They were grown from the melt in a 
graphite mould. The crystals so prepared, #5” in 
diameter and about 6” in length, are quite 
homogeneous in composition throughout the 
whole length, but the internal friction back- 
ground of such a crystal is unusually high. Subse- 
quent annealing at 800°C has the effect of re- 
ducing partly the internal friction background, 
but it was still several times higher than that of 
the polycrystalline specimen (see Fig. 1). It was 
thus decided to use polycrystalline specimens 
having very large grain size with the hope that 
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Fic. 6. Relaxation of shear modulus in polycrystalline 
alpha-brass. 


10 A, E. Arkel, Metallwirtschaft 7, 656 (1928); see also 
R. F. Mehl, Trans. A.I.M.E. 122, 11 (1936). 
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Fic. 7. Effect of frequency of vibration on internal 
friction associated with the preferential orientation of zinc 
atom pairs in alpha brass. Frequency of vibration: ©, 
about 0.6 c.p.s.; X, about 2.3 c.p.s. 


the internal friction due to the grain boundary 
slip will be shifted sufficiently to higher tempera- 
tures and be separated from the internal friction 
due to the preferential orientation. The alpha- 
brass wire used in the study of grain boundary 
relaxation (Section I) was sealed in a thin quartz 
tube filled with }j-atmosphere of argon (99.8 
percent pure). This was annealed at 800°C for 5 
hours. After such a treatment the grains in the 
wire specimen were much larger than the diame- 
ter of the wire and most of them were extended 
clear across the specimen. The concentration of 
zinc in the specimen after such a treatment was 
found by chemical analysis to be 29 percent, 
indicating that little zinc was lost by this treat- 
ment. The internal friction curve of this wire is 
shown in Fig. 7, in which the internal friction 
concerned is fairly well separated from the in- 
ternal friction due to grain boundary slip. Using 
two frequencies of vibration having a ratio of 
_ 3.89 (average value over the temperature range 
studied), it was found that the heat of activation 
associated with the preferential orientation of 
zinc atom-pairs in alpha-brass is 40,000 calories 
per mole within 5 percent. This agrees well, 
within experimental error, with the value of the 
activation energy for the diffusion of zinc in 
alpha-brass mentioned above. 


290 








IV. GRAIN BOUNDARY SLIP AND 
ATOMIC DIFFUSION 

The activation energy associated with the 
grain boundary slip in 70-30 alpha-brass has been 
determined in Section I as 41,000 calories per 
mole. Using the same specimen, we have also 
shown in Section II that the activation energy 
associated with the preferential orientation of 
pairs of zinc atoms in alpha-brass, which is a 
diffusion phenomenon, is 40,000 calories per 
mole. Both values are equal, within experimental 
error, to the activation energy for the diffusion of 
zinc atoms in alpha-brass ¢29.08 percent zinc). It 
is believed that the experimental accuracy in the 
determination of the activation energies de- 
scribed in Sections | and II are accurate to 5 
percent and, therefore, such a close agreement 
deserves special attention. This would indicate 
that the grain boundary slip in alpha-brass as 
manifested in the anelastic measurements is a 
process of atomic diffusion. The identity of the 
activation energies for grain boundary slip and 
for volume diffusion is somewhat unexpected, be- 
cause it is a conventional belief that grain 
boundary diffusion is faster than volume diffu- 
sion and the activation energy for the latter is 
higher, for there the atoms are held in greater 
restraint and thus a greater energy of activation 
is required to set them free."' However, this 
finding is consistent with Elam’s” observation 
that zinc does not penetrate preferentially along 
the grain boundaries of copper at a temperature 
below 500°C. 

In conclusion, it is a pleasure to the author to 
acknowledge his indebtedness to Professors 
Clarence M. Zener and Cyril S. Smith of this 
Institute for their this 
problem. 

Note added in proof: The identity of the activa- 
tion energies for grain boundary slip and for 
volume diffusion has also been found in the case 
of alpha-iron (T. S. Ké, communicated to Metals 
Technology). A discussion on the possible in- 
ference as to the structure of grain boundaries in 
metals is lately given by this author.” 

1 R. F. Mehl, Trans. A.I.M.E. 122, 11 (1936). 


2C. F. Elam, J. Inst. Metals 43, 217 (1930). 
8 T.S. Ké, Phys. Rev. 73, 267 (1948). 
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The Plastic Flow of Thick-Walled Tubes with Large Strains’ 


C. W. MacGrecor,? L. F. Corrin, Jr.,3 anv J. C. FisHer‘ 
(Received October 20, 1947) 


The theory for the plastic flow of thick-walled tubes is discussed for a material having an 
arbitrary true stress-strain curve, loaded by any combination of internal pressure, external 
pressure, and end load, where the strains may become large necessitating the use of finite true 


strain definitions. 


Comparisons are made of the finite strain theory, the infinitesimal strain theory, and ex- 
perimental results for the special case of a closed end thick-walled tube of aluminum. The com- 
parison shows that the infinitesimal strain theory is considerably in error even at moderate 
strain values. The finite strain theory appears to be satisfactory even for strains somewhat 
beyond the value for the local bulging of thick-walled tubes. 


FINITE strain theory for the plastic flow 

of thick-walled tubes to large deformations 
was developed during the War by Bridgman.® An 
independent solution of the same problem was 
obtained later by the present authors.* A com- 
parison of the two solutions now reveals that they 
are identical save for minor differences of termi- 
nology and method of treatment. This duplica- 
tion of effort resulting from the limited flow of 
scientific information at the time has only re- 
cently been discovered. 

The present paper gives experimental data 
from the expansion of a thick-walled tube under 
internal pressure indicating the superiority of the 
finite strain theory over the infinitesimal strain 
theory with which it is compared. In view of the 
limited availability of the original references it 
was felt that an introductory review of the 
derivation of the finite strain theory would be of 
value. 


! This paper is the result of a portion of a research pro- 
gram carried out at the Massachusetts Institute of Tech- 
nology for the Bureau of Ordnance, Navy Department. 
Approval for the release of this information has been re- 
ceived from the Chief of the Bureau of Ordnance. State- 
ments expressed are those of the authors and should in no 
way be construed as the opinion of the United States Navy. 

* Professor of Applied Mechanics, Department of Me- 
chanical Engineering, Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts. 

3 Assistant Professor of Mechanical Engineering, De- 
partment of Mechanical Engineering, Massachusetts Insti- 
tute of Technology, Cambridge, Massachusetts. 

‘Instructor, Department of Mechanical Engineering, 
Massachusetts Institute of Technology, Cambridge, Mas- 
sachusetts. 

5 P. W. Bridgman, Watertown Arsenal Laboratory Re- 
port Number WAL 111/7-3, October 8, 1943. 

6C. W. MacGregor, L. F. Coffin, Jr., and J. C. Fisher, 
“Project G Report No. 1,” submitted to the Chief of the 
Bureau of Ordnance, Bureau of Ordnance, Navy Depart- 
ment, Washington, D. C., February 10, 1945. 
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REVIEW OF THEORY 


In problems concerned with the deformation of 
materials in the elastic and in the early plastic 
ranges, the magnitudes of the displacements 
encountered are small compared to the dimen- 
sions of the deforming element. For such prob- 
lems it is sufficiently accurate to employ the 
infinitesimal strain concept wherein the deforma- 
tions are referred to the original dimensions of the 
element rather than to the instantaneous values. 
As the magnitude of: the plastic deformation 
increases relative to the dimensions of the ele- 
ment,’ a condition is soon reached where an 
analysis of the state of stress and strain requires 
the application of the finite strain concept. 

Such a situation is presented in the problem of 
the finite plastic deformation of a thick-walled 
tube subjected simultaneously to an internal 
pressure, an external pressure, and an end load.* 
The need for a solution of this problem is three- 
fold: (a) to aid the designer of industrial equip- 
ment wherever use is to be made of the added 
strength available through the heavy cold- 
working of tubes; (b) to serve as a guide for the 
pressure versus expansion relationship during the 
cold working operation; and (c) to provide a 
means for determining the pressure at which 
failure by local bulging will occur in thick-walled 
pressure vessels. 

The radial equilibrium of an element of ma- 


* The cases of the partially yielded tube and the tube 
just completely yielded have been treated under the 
assumption of infinitesimal strains (see references 7 and 
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Fic. 1. Expansion of tube to very large deformations. 


terial in an expanded tube demands that 
de, 
(r+-u)——=01-<,,"* (1) 
d(r+u) 
where r+u is the radial distance to the element 
in the expanded tube. 

The cross section of a tube with original inside 
and outside radii Ry and R, respectively, ex- 
panded to internal and external radii Ro+u» and 
Ri+4, respectively, is shown in Fig. 1. An arbi- 
trary element at the radius r has moved to the 
radius r+u. The true radial and tangential 
strains become 


d(r+u) ; (14— 
€- = ee oe Oe ae 
dr d 


(2) 
(r+u)do 


r 
Pm 
¢,=log————___= le. 1 +") 
rd r 


After the elimination of u in Eqs. (2), the equa- 
tion for the compatability of the strains is 





de: 
r—=etr-“t—1, (3) 
dr 
Equation (1) can be rewritten as 
du 
i+— 
do, dr 
r = (¢,—0,)———_ = (o,—a,)e"" 7s (4) 
r u 
i+-— 
r 


if use is made of the relationships given in 
* Eqs. (2). Thus Eqs. (3) and (4) are the equations 
of strain compatability and of equilibrium re- 
spectively for the large deformations of a thick- 
walled tube. 


* 


** See Nomenclature for definition of variables. 
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It is further assumed that (a) the material 
obeys the distortion energy law of yielding, (b) 
the volume remains constant, (c) the ratios of 
principal stress differences to the corresponding 
principal strain differences are equal,*** and 
(d) the axial strain is constant over the cross 
section. These are the commonly accepted laws 
of plasticity as applied to problems of this type.” ® 

The distortion energy law of plasticity states 
that the effective stress 


o = (v2/2)[(o1—02)? + (¢2—03)?+(¢3—01)?]' (5) 


required to maintain plastic flow in an element of 
material is a function of the effective strain 


e= (V2/3)[(€1— €2)?+(€2—€3)?+(€s—€1)?]! (6) 


associated with that element, where o1, a2, o3; 
€1, €2, €; are the principal stresses and strains, re- 
spectively. This relationship can be stated 
symbolically as 


o=fi(e). (7) 


The constants in front of the radicals in Eqs. (5) 
and (6) were so chosen that for the tension test 
Eq. (7) reduces to the true stress-strain relation- 
ship o,= f(e,). 

Since €1, €2, and e; represent true plastic strains, 
the constancy of volume can be expressed by the 
equation 


eiteste;=0 (8) 


independently of the magnitude of plastic de- 
formation when there is no rotation of the princi- 
pal strain axes. 
From assumption (c) and Eq. (8) it follows 
that 
Ee,-=D([e-—4(o.+<¢.) | 
Ee=D[or—3(0.+0,) ]¢, (9) 
Ee.=D(o.—}(0-+0:) | 


where D is a function of the radius. 
Eqs. (5), (6), and (9) it follows that 


Ee 


From 


*** This is a special case of the more fundamental as- 
sumption that the ratios of principal stress differences to 
the corresponding principal strain rate differences are equal 
during plastic flow. 

7A. Nadai, Plasticity (McGraw-Hill Book Company 
Inc., New York, 1931). 

8C. W. MacGregor, L. F. Coffin, Jr., and J. C, Fisher, 
Partially Plastic Thick-Walled Tubes, to be published 
shortly. 
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Inserting this value of D in Eqs. (9) gives 


€ 
a 
o 


cH 
. 


es ee (10) 
oC 





€ 
e,=—[o:—43(¢,+<:) | | 
o J 
Let e,=k=constant, satisfying the fourth as- 
sumption (d) above. Then from Eqs. (10) 


ok 
o:,=—+3(o,+0,). (11) 
€ 


Substituting this value of ¢, into E 


effective stress, 
s=(1-=), (12) 


where for a a. tube the upper sign 
corresponds to an excess of internal over external 
pressure, and the lower to an excess of external 
over internal pressure. 


q. (5) for the 


Ce oe = 


The value of o,—o, from Eq. (12) can be 
substituted in the equation of equilibrium (4) 
which becomes 


do, 20 k?\3 
a. 
dr v3 é 


Integrating this equation the — stress is 


o.= [do-+00= —pit 





(13) 


do, 


r=R, 


2 ph 
-—p.— f (1-=) rd 


where the upper sign applies to an excess of 
internal pressure and the negative external pres- 
sure — fp; appears as the constant of integration. 

Expressions for a, and o, are obtained from 
Eqs. (11) and (12) as 


—_ -=)] 


=}(o,+01) 


(14) 


a (15) 





4 


t The variable of integration could equally well be taken 
as r+u, in which case Eq. (14) would appear as 


om peed, fo" of 1 Eye. 


e r+u (14a) 
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Considering now the plastic strains, it follows 
from Eq. (8) that 


€- = —k— «&%. 


Substituting this value of e, into the compatability 
Eq. (3) 


de: 
r—=e~ht2e) 1, 
dr 
Integrating, 
c 
é,=4 toe | +=] (16) 
r- 
where ¢ is a constant of integration. 
Defining 
c 
= log ++], 
r2 
the strains become 
E=-— k = ¥3 
€:= 32 (17) 
e:=k 


The value of the effective strain can be com- 
puted from the values of €,, €:, and ¢, in Eqs. (17) 
according to the definition of ¢ in Eq. (6). Thus 

e= (k?+43(k+2)?)?. (18) 

The expressions for the stresses and strains 

contained in Eqs. (7), (14), (15), (17), and (18) 


are summarized below: 


2 Ri k?\3 dr ) 
nowt dpe 
v3 é r 


r 


2 k?\3 dr 
o,.= — 2 F— f o(: --) etre 
v3 a r 


(19) 





€ € 
¢,-= —k—}2 
€,=42 
é,=k 
=f(e) J 
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where 


e=(k?+ a(k-+s)*)4) 
toe | +] || 


u 

—=et—1 

r 
and where c is constant for a given deformation 
of the tube. 

The integrals in Eqs. (19) can be evaluated 
readily by numerical methods when the relation- 
ship «= f(¢) has been determined from a tension 
test and when the loading system is specified. In 
general, several trial solutions are needed in order 
to determine which pair of values for ¢ and k 
satisfies the boundary conditions. 

A particular loading condition of practical 
interest occurs when the tube has closed ends, 
and is subjected to no end load in addition to that 
due to the internal pressure p>» and the external 
pressure /;. It will be shown that for this loading 
the axial strain is e,=k=0. 

When k=0 Eqs. (19) and (20) reduce to 


Ri gd(r*) | 
o,= “ars aS . 











rete | 
Rt gd(r*) 261 
v= —pF— f a 
34, r+e v3! 
1 Rigd(r*?) o | 
o,= —pPiF— +—r (21) 
v3 re+e v3 | 
€¢-= —}z | 
€,=42 | 
¢,-=0 | 
o = f(e) 
_ where 
2 
€«=— 
v3 
c]| 
£=log 1 +2 - (22) 
yr? 
u 
—=e'—1 
rT J 





The expression e**t has been replaced by its 
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(20) 


ea) a a er 


value 


in determining the above expressions. 

From Eqs. (21) it is seen that 2€,=z. Substi- 
tuting the value of e, from Eqs. (2) and the value 
of zs from Eqs. (22) in this expression 


BY 


u\? c 
oe. 1 +") = oe. 1 +=) 
r a 


from which 
(r+u)?=r?+c (23) 
and 
d(r+u)?=d(r)?. (24) 
By means of these relationships the stress 
Eqs. (21) can be reduced to 


2 Ritui gd(r+u) ) 
o,= —piF— f —— | 


V3 Yu r+u 
ut od(r+u) 2e 
ope f —+—p>. (25) 
ee r+u v3 


2 Ry+ 
G:>= -—Ar~ ~ f 
v3 


se ar v3 


The end load acting on the tube 
evaluated as 


can now be 


Ritu 


F= 2no(r+u)d(r+u). 


Ro+uo 


(26) 


Let 


p=r+u, ( 


then the expression 


bho 
~ 


F—rpo*potmpr"pr 


becomes ’ 
_ Tpo’ Pot rps" pi 
pl 2 Plg oC 
te f"|- pre “apse 
po v3 


p Pp v3 


. Plg 
+| -pe— J “dp |rou-+x0'P 
V3 Joo p 
4r 7"( fo 
= - Bil “dp| dp 
3 p 


oe pl 2rp.” fa 
—_— opdp+ — f —dp. 
v3 p 


(28) 
po v3 0 p 


JOURNAL OF APPLIED PHYSICS 








) 





Evaluating by parts, the first integral in Eq. (28) 
becomes 


Plg 2 Plg pl 
Sis “ip}ode= —™ ff “det { opdp. 
p p0 


Substitution of this value in Eq. (28) results in 
F— po potrmp:*pi=0 
F=n2(Rot+uo)*po—m(Rit1)*p1 (29) 


or 


corresponding to a closed end tube, and justifying 
the original assumption that e,=0 for such a 
tube.T 

For a material with a constant yield stress 
a=a0, Eqs. (25) can be integrated directly, 
vielding 


2c0 Ritu, | 
o,= — pi +F— log.——_ | 
v3 r+u 
Ritnm 
o,;= prt —“(1 Hoe. ~~) 
r+u 
Ritu 
o:= ~pt—(1 — log, ~~) , (30) 
r+u 
cg "Gy 
¢,=9 
u 
—ae*t—] 
r 





4 


Using the upper signs, the expressions for the 
stresses are identical with those given by Nadai’ 
for the case of internal pressure; except that 
Ri+u, and r+u replace R; and r, respectively, 
and the external pressure ~; is included. This 
indicates quite reasonably that Nadai’s stress 
expressions are valid for large deformations of a 
material exhibiting a constant yield stress if the 
actual instantaneous radii are used.ff 


+ Both Bridgman’s tests on the collapsing of closed end 
steel tubes under external pressure and the tests described 
below for the expansion of a closed end aluminum tube 
check experimentally the value «,=0. 

tt For the more common case occurring in practice for 
large strains, namely, the arbitrary stress-strain curve, it is 
not sufficient to replace R; and r by Ritu, and r+-4, re- 
spectively, in NAdai’s infinitesimal strain solution to obtain 
agreement with the present theory. Agreement is obtained 
only when o=ao, since for ¢=gp the strain distribution is 
irrelevant; when o varies with e, however, the infinitesimal 
strain theory gives incorrect values of the stresses as it 
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COMPARISON WITH INFINITESIMAL STRAIN 
THEORY AND EXPERIMENT 


It is of interest to compare the predictions of 
the finite strain theory with those of a theory 
based upon the assumption of infinitesimal strain. 
Such a comparison has been made for the pres- 
sure versus external tangential strain curve of a 
closed end 17SO aluminum tube with wall ratio 
R;/Ro=2.0 and outside diameter of 5 inches. A 
typical stress-strain curve for the aluminum is 
shown in Fig. 2. 

The curve predicted by the new theory was 
computed by numerical integration of the expres- 
sion in Eqs. (21) for o,, assuming p;=0. The 
integration was carried out for several values of c, 
and the corresponding pairs of values of external 
tangential strain 


Cc 
E11 =} log-( 1+) 
2 


and bore pressure Pp = —o;o are plotted in Fig. 3. 
The curve predicted by the infinitesimal strain 





STRESS-STRAIN CURVE 
FOR 17SO ALUMINUM 
SPECIMEN CUT IN 
TANGENTIAL DIRECTION 


TRUE STRESS IN THOUSANDS OF PSI 
~ 


TRUE STRAIN 
Fic. 2. 


bases them on approximate values of « in the relationship 
o=f(e). This difficulty is not removed by introducing the 
actual instantaneous radii. 
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theory was computed from the equations 


BORE PRESSURE IN THOUSANDS OF P.S.I. 
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1 pe 
o= —p— f 
v3 J, 
1 fe 
v= —pF— f 
; v3 v7, 
1 Ri 
n= —F— f 
v3 7, 
C} 
é3-— 
r? 
C1 
é=-— 
r? 
e,=0 
2 C1 
= . 
v3 r? 
a= fie) 











od(r*) 
r- 
od(r*) 2¢ 
=~ 
r? v3 
od(r*) o@ 
=— 
r? v3 
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FOR ALUMINUM TUBE 
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which correspond to the infinitesimal strain as- 
sumption.’ In this case as before ~; was con- 
sidered zero. A similar method of numerical 
integration was employed for the determination 
of the bore pressure corresponding to several 
values of c,, and the pairs of po, €1 values so 
found are indicated in Fig. 3. 

Figure 3 also gives the experimentally de- 
termined curve of pressure versus external tan- 
gential strain for the aluminum tube. It is evident 
that the infinitesimal strain theory indicates an 
indefinitely increasing bore pressure, while the 
large strain theory and the experimental curve 
both indicate a maximum pressure corresponding 
to a strain of about 0.04, after which the tube 
bulges locally and the pressure falls off. At strains 
larger than about 0.01 the infinitesimal strain 
theory departs considerably from the large strain 
theory and from experiment. 

No great significance can be attached to the 
smaller disagreement between the predictions of 
the large strain theory and the actual behavior of 
the aluminum tube. Extensive study of the ma- 
terial of the tube shows it to be highly anisotropic, 
with substantial differences in the stress-strain 
curves for axially and tangentially oriented speci- 
mens. It is not likely that the simple relationship 
o = f(€) as defined in this paper will apply exactly 
to anisotropic materials. 
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NOMENCLATURE 


c Constant of integration 

c¢, Constant of integration 

D_ Variable coefficient relating plastic strains and stresses 
E Modulus of elasticity 


._F End load 


k Constant value of e 

fo Internal pressure 

pi External pressure 

r Radial distance to an element in the undeformed tube 
Ro Internal radius of the undeformed tube 
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R, External radius of the undeformed tube 

u Radial displacement of an element of initial radius r 
during plastic yielding of the tube 

ug Radial displacement of the internal surface of the tube 

u; Radial displacement of the external surface of the tube 

z  Log-(e*+c/r*) 

3 Radial, tangential, and axial strains, respectively, corre- 
sponding to an element at radial distance r+u in the 

e-} deformed tube 


€1 
“f Principal strains 


v2 
€ Effective strain, «= 3 ((€1 —€2)?+ (€2—€3)?+ (€3 —€1)*)! 


p r+u 

po Rotuo 

po Ritu 

o,| Radial, tangential, and axial stresses, respectively, cor- 
o:¢ responding to the radial distance r+u in the expanded 
oz} tube 

oro —Po 

a1 

o2¢ Principal stresses 

o3 


2 
o  Effectivestress, o= aC —a2)?+ (o2—a3)*+(03—0;)*)! 


oo Constant yield stress 
@ Angle 





Gyrograms for Simple Harmonic Systems Subjected to External Forces* 


V. ROJANSKY 
Department of Physics, Union College, Schenectady, New York 
(Received November 17, 1947) 


A graphical method is described for the approximate integration of equations of the form 
md?x /dt?+-kx =f(t), which arise in several branches of physics and engineering. The graphs, 
named ‘‘gyrograms,” are an extension of the familiar ‘reference circle’’ used in representing 
simple harmonic motion. The method is applicable whenever f(t) can be adequately approxi- 


mated by a step function having a few steps. 


INTRODUCTION 


UR purpose is to describe a graphical method 
for integrating the equation 


9 


Pace wey (1) 
dt? 


when f(t) is a step function having a few steps, or 
when it can be adequately approximated by such 
a step function. For definiteness, we shall regard 
Eq. (1) as referring to the motion of the friction- 
less one-dimensional elastic system shown in 
Fig. 1 and supposed to be so far removed from the 
earth that the gravitational forces acting on it are 
negligible, and we shall take m to be the mass of 
the bob, k the combined restoring constant of the 
two (massless) springs, x the (vertical) displace- 


* The author is indebted to the Chief Scientific Adviser, 
Ministry of Works (London), for permission to publish this 
material, which was prepared in 1944, when the author was 
a member of the staff of the Research and Experiments 
Department of the Ministry of Home Security (London) as 
a representative of the Princeton University Station of 
Division 2, NDRC, under OSRD contract OE Msr-260. 
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ment of the bob from its free equilibrium position 
(that is, the equilibrium position in the absence 
of external forces, namely, forces other than 
those that are exerted by the springs), and f(t) 
the external force acting on the bob and de- 
pending explicitly only on the time ¢. But the 
method applies to any equation of the form (1), 
irrespective of its physical interpretation.! 

The present method is rather elementary and 
can hardly be new; yet no references to it seem to 
be available.” We shall describe it by examples. 


An equation of this form governs, for example, the 
bending moment along an axially compressed uniform 
beam; for a graphical method of handling the boundary 
conditions that then arise, see V. Rojansky and R. A. 
Beth, ‘““Camptograms for beams in compression,” Trans 
Am. Soc. Mech. Eng. 69, A-202 (1947). 

? For a recent method for the graphical solution of equa- 
tions of the form (1), see G. E. > a mnt “A method for 
estimating equivalent static loads in simple elastic struc- 
tures,” The David W. Taylor Model Basin, U. S. Navy, 
Report No. 507 (June 1943). Hudson’s method is based on 
a method of Kelvin? and of Meissner.‘ 

Lord Kelvin, Phil. Mag. 34, 433 (1892). This method 
was designed primarily to treat the motion of a particle 
in a plane. 

4E. Meissner, Schweizerische Bauzeitung (Ziirich) 104, 
35 (1934). 
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F 7. 


Fic. 1. A simple harmonic system and a reference circle; 
the circle is drawn for the case of free oscillation between 
the limits indicated by the dashed lines. 


FREE OSCILIATION 


We begin with the “reference circle,’ familiar 
from the elementary theory of simple harmonic 
motion. Let no external forces act on the bob in 
Fig. 1, and let the bob oscillate, about its free 
equilibrium position, between the limits indicated 
by the dashed lines. The corresponding reference 
circle is shown at the right-hand side of the figure. 
In the coordinate frame accompanying this circle, 
the ordinate x represents the displacement of the 
bob from the point x=0, and the abscissa u, 
which has the physical dimensionality of length, 
is given by the equation 


u=v/w, (2) 
where v is the velocity of the bob and 
w=(k/m)!; (3) 
that is, 
u=dx/dd, (4) 


“4 


Upr%o 








Fic. 2. Gyrogram for a free oscillation starting with the 
values x =X»9 and v=wto, where #9 =xXo>0. 
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where 

o=wt. (5) 
When expressed in terms of x and u, the energy E 
of the system, namely, $kx?+4mv*, takes a 
symmetric form: 


E=$k(x?+u’). (6) 


The natural period of the system will be denoted 
by T: 

T =2r/w. (7) 
Each of the quantities x,.v, and presently f(é), 
will be regarded as positive if it is directed 
upward. The plane of the reference circle will be 
called the xu plane. 

The basic property of the reference circle is 
that, while the bob oscillates freely, its displace- 
ment x is at all times equal to the vertical com- 
ponent of the displacement of a “representative 
point,’’ which moves in the xu plane counter- 
clockwise along the circumference of this circle 
with the constant angular velocity w. 

Let the motion begin at t=0 with the initial 
values x=Xo and v=%, so that u%p»=vo/w. Then 
the initial coordinates of the representative point 
in the xu plane are xo and uo; hence the appro- 
priate reference circle (whose center lies at QO) 
passes through this point, and therefore the 
representative point begins its motion as shown 
by the circular arc in Fig. 2, drawn for a case 
when xo is positive and u»=xo. This arc we call 
the ‘‘gyrogram” of the motion. The gyrogram of 
Fig. 2 subtends at O the angle of 120° and thus 
refers to the time interval lasting from t=0 to 
t=4T7; it can, of course, be extended along the 
reference circle indefinitely. 


x } 


ie 


U 
n 0 
{/k 











Fic. 3. Gyrogram for a case of a constant force suddenly 
applied to a system at rest. 
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CONSTANT EXTERNAL FORCE 


If the bob is acted upon by a constant (vertical) 
external force (for example, the force caused by 
gravity), so that 


f(t) =f=constant, (8) 


the motion remains simple harmonic and its 
period is still given by Eqs. (7) and (3). The 
equilibrium position, however, no longer lies at 
x=0, but at the point x=f/k, at which the 
springs balance the external force. The motion 
can thus still be described by a reference circle, 
whose center, however, now lies not at the origin 
O but at a point O’, say, the coordinates of which 
in the xu plane are 


x=f/k, u=0. (9) 


SUDDEN APPLICATION OF A CONSTANT 
EXTERNAL FORCE TO A SYSTEM 
AT REST 


While the bob is at rest in the absence of ex- 
ternal forces, the representative point stays at 
the origin O of the xu plane. At the instant when 
the constant force f is applied, the representative 
point starts moving with the constant angular 
velocity w along the appropriate reference circle. 
The center of this circle lies at the point O’, 
having the coordinates x=f/k and u=0; this 
circle passes through the point O, where the 
representative point is originally located, and 
hence the gyrogram has the form shown in 
Fig. 3, drawn for a case when f is directed 
upward. 


SUDDEN APPLICATION OF A CONSTANT 
EXTERNAL FORCE TO A SYSTEM 
IN MOTION 


Let the motion begin at the instant ¢=0 with 
the initial values x=x»9 and v=vo, and continue 
in the absence of external forces for one-third of 
the period, but at the instant t=t,;=}3T let a 
constant force f be suddenly applied to the bob. 
A gyrogram for such a case is shown in Fig. 4. 
The arc beginning at the point (#o, x9) and ending 
at the point marked f¢; is taken from Fig. 2; it 
refers to the first one-third of the period, when 
the motion is free and when the center of the 
reference circle lies at O. The rest of the gyro- 
gram has its center of curvature at O’, namely, 
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Fic. 4. Gyrogram for a case of a constant force suddenly 
applied to a system in motion. 


at the point whose coordinates are x=f/k and 
u=0, as in Fig. 3. 

The fact that the two arcs of the gyrogram in 
Fig. 4 join at the point marked ¢,; represents 
graphically the fact that x and v vary continu- 
ously at the instant ¢,, while the kink in the 
gyrogram at this point corresponds to the fact 
that the acceleration of the bob, as well as the 
higher derivatives of x with respect to ¢, change 
discontinuously at the instant 4. 


A STEP-LIKE EXTERNAL FORCE 


Let the bob be at rest until the instant ¢=0, 
when the step-like external force shown in Fig. 5 
is applied to it. The interval OT in the figure 
marks the assumed natural period of the system. 
The symbol fo: denotes the constant external 
force acting from t=0 to ¢=t,, and so on. The 
round-figure numerical data that we adopt in 
this example are given in Table 1.5 

We need not know the actual values of ¢1, fz, 
and so on, but only their ratios to the period T 
(which corresponds to 360°); accordingly, we 
give in the second column of Table I the assumed 
values of ¢;, $2, and so on (that is, wt1, we, and so 
on), expressed in degrees. 

Similarly, we need not know the actual values 
of foi, fiz, and so on, but only their ratios to the 
restoring constant k; the assumed values of these 
ratios are given in the last column of Table I. 

The graphical procedure is now as follows 
(Fig. 6): 


5 We use the c.g.s. system of units and express m in 
grams, k in dynes per centimeter, f in dynes, v in centime- 
ters per second, and x, u, and f/k in centimeters. 
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Fic. 5. A step-like external force. 


1. On a sheet of rectangular graph paper, 
draw the vertical x axis and the horizontal u axis, 
directed as shown, and choose a scale unit (the 
same unit for x and u) to represent a centimeter. 

2. Along the x axis mark the points Oo, O1, 
and so on, whose respective ordinates are equal to 
the values of f/k listed in Table I (for example, 
the ordinate of Oo; equals fo:/k, namely, 2.5 cm). 

3. Using the values of A@ listed in Table I, 
draw counterclockwise the following circular 
arcs: the arc Ot, (center Oo:, radius Oo,O0, angle 
30°) ; the arc tytz (center Oj2, radius Oj2l;, angle 
50°); the arc fet; (center O23, radius Oozfe, angle 
90°); the arc tst, (center O34, radius Ogzgfs, 
angle 60°); and, finally, the thin-line circle 
(center O, radius Ot,) which refers to the motion 
after the instant ¢4. 

The gyrogram, which enables us to estimate 
the values of x and v at any instant ¢, is now com- 
plete. We note, for example, that the first maxi- 

















1 scale unit = | cm 





Fic. 6. Gyrogram for the case of the external force de- 
scribed in Fig. 5 and Table I. 
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mum of x equals about 4.0 cm. To estimate the 
time of its occurrence, we add the angle between 
Oost and the x axis (namely, about 52°) to do, 
getting about 132°, or about 0.37 360°; hence 
the first maximum of x occurs at about the 
instant £=0.37T. The first maximum of u occurs 
at t=, and equals about 3.2 cm; hence the first 
maximum of v occurs at t=f. and equals about 
3.2~ cm/sec. And so on. 


SUMMARY OF METHOD 


If f(t) has the constant value f. during an 
interval extending from ¢=/,to0 t=, the segment 
of the gyrogram for this interval is a counter- 
clockwise circular arc; its center of curvature lies 
at the point Ow, whose coordinates are x =fa/k 
and u=0, and the angle ¢ that it subtends at 
Ow equals w(t—t,) radians. This arc begins at the 
terminus of the preceding segment of the gyro- 
gram, or, if ¢, is the initial instant of the motion, 
this arc begins at the point in the xu plane 
specified by the initial conditions. ; 

In describing the method, we have relied on an 
intuitive appeal to the familiar ‘reference 
circle.’’ To justify it more properly, two steps are 
necessary : first, to show that, if f is constant, the 
first integral of Eq. (1) can be written as 


(x—f/k)?+u? =constant =r’, (10) 


which is the equation of a circle in the xu plane; 
and, secondly, to show that, if f is constant, the 
differential of arc of the gyrogram, namely, 
[ (dx)?+ (du)? }!, equals rd¢, where r is the con- 
stant of integration appearing in (10), and ¢ is 
expressed in radians. The proofs are left to the 
reader. 


CONCLUDING REMARKS 


If f(t) is a step function, the present method is 
in principle exact and is subject only to the errors 
of draftsmanship. If f(t) is continuous, as is 
usually the case in practice, the method can still 
be used, provided that f(t) can be adequately 
approximated by a step function having com- 
fortably few stes ; but then, of course, analytical 
errors creep into the work. The decision -on the 
adequacy of the approximation is best left to the 
computer familiar with the practical aspects of 
the problem. The effect of the lengths and shapes 
of the steps of the approximating step function on 
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the graphical results can of course be studied by 
trial computations. 

If all the Ad’s involved in a computation are 
equal, the graphical work can be simplified with 
the help of a wedge of angle Ag, cut out from a 
sheet of polar graph paper so that its vertex lies 
at the printed origin of the polar coordinate 
system. After marking their centers of curvature 
on the work-sheet, the terminal points of the 
various circular arcs of the gyrogram can then be 
located by means of this wedge alone, without 
using either ruler, compasses, or protractor, and 
the arcs themselves can be drawn in later, using 
only compasses. A set of such wedges is con- 
venient to have on hand even if not all the Ad’s 
are equal, but if some of them recur sufficiently 
often. 

In certain problems the arithmetic can be 
shortened by replacing x and u by variables 
proportional to them ; for example, if equivalent 
static loads are being estimated, it is convenient 
to use kx and ku, rather than x and uy, for the 
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t f f/k 
sec o Ad dynes cm 
0 0° 

30° ihe 2.5 
ty 30° 

50° te 3.5 
te 80° 

90° te 0.0 
ts 1 70° 

60° _™ —1.5 
ts 230° 

2) io a 
vertical and the horizontal coordinate, re- 


spectively. 

Analytic expressions for various quantities can 
be obtained from gyrograms trigonometrically. 
For example, in Fig. 4 the first and the second ex- 
treme values of x arer and fk~! — (r?—2fk-'r sin15° 
+f?k-*)}, respectively; here r = (xo?+u¢7)!. 

I am indebted to Professor W. N. Thomas, 
Dr. F. Garwood, Mr. R. J. Hansen, and Pro- 
fessor H. F. Bohnenblust for helpful discussions. 
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Large Plastic Flow and the Collapse of Hollow Cylinders 


P. W. BripGMAN 
The Lyman Laboratory of Physics, Harvard University, Cambridge, Massachusetts 


(Received November 20, 1947) 


The collapse of heavy steel cylinders under external pres- 
sures running up to 412,000 p.s.i. was studied during the 
war under contract with the Watertown Arsenal. The dis- 
tortions are so large that the conventional small strain 
analysis is not adequate. The presence of large strains is 
recognized in the analysis by writing the stress equations 
of equilibrium in terms of the final displaced positions of 
the particles instead of in terms of the initial positions. The 
conventional assumption of negligible volume change is 
retained. Attention is called to a purely geometrical condi- 
tion on the strains. In the experiments the cylinders showed 
no change of length. The generality of the conditions under 


ROFESSOR MacGregor was so kind as to 

let me see the manuscript of the preceding 
paper at the same time that he was submitting 
it for publication in this journal. At that time 
he made the suggestion that, in view of the close 
connection between our respective investigations 
and the comparative inaccessibility of my own 
publication, it would be of interest if I should 
outline here the general character of my results. 
The editor, Dr. Hutchisson, has fallen in with 
this suggestion, and I am glad to avail myself of 
the opportunity thus offered. My original work 
was done during the war under contract with 
the Watertown Arsenal. It was unclassified and 
permission has been given for independent and 
duplicate publication. 

With regard to the theoretical discussions, 
Professor McGregor has already intimated that 
they are essentially the same. The details, and 
in particular the order of the analysis, were not 
always the same however, and it seems to me 
worth while to restate certain of my theoretical 
results in view of the contemplated application 
to experiment, and in order that the degree of 
generality of the assumptions on which the 
results depend may be more apparent. 

Jn the first place, there is a purely geometrical 
condition on the strains (large) independent of 
any connection between stress and strain. This 
is, in MacGregor’s notation: 


u=([(r+u)?(1+e,’)—c }i—(r+u). 
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which this may be expected is examined. A linear relation 
was found experimentally between collapsing pressure and 
log (rext/Tiat). Furthermore, P/log(rext/rint) is not constant 
but rises continually with increasing collapsing pressure. 
These results indicate that the strain hardening increases 
continually with increasing distortion, but at a rate two- or 
threefold less than would be given by experiments in simple 
tension. The collapse was pushed experimentally to such 
extremes that the geometrical symmetry was lost because 
of the grain structure. This limits the region of applicability 
of any simple theory. 


Here c is a constant of integration determined by 
the boundary conditions. In the deduction of 
this equation the only assumptions are that the 
strains satisfy the condition of constancy of 
volume 


€-te-te,=0 


and that the ordinary (not natural or logarithmic) 
axial strain e,’ is constant across the section and 
along the length, a constancy demanded by 
considerations of symmetry. Incidentally, the 
formulas become somewhat simpler if expressed 
in terms of the coordinates of the finally displaced 
particles instead of in terms of the initial posi- 
tions, as above. 

It is an experimental result, both in Professor 
MacGregor’s and my own work, that there is no 
change of length, or e,/=0=e,. In Professor 
MacGregor’s analysis this result was deduced 
under the assumption of the distortion energy 
law of yielding. The same result holds, however, 
under somewhat less restrictive conditions, 
namely, if the relation between stress and strain 
is: 

€-=alo,—}(or+¢,) ], 
e:=alo.—3(o.+0,) ], 
€,=alo,—3(¢,+0:) |. 


This is the integrated form of the relation which 
holds in the case of isotropy. These equations 
and the corresponding assumptions embrace both 
the strain hardening case, in which a is a function 
of stress or strain or both, as in MacGregor’s 
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analysis, and the non-strain hardening case, 
which is a degenerate case in which a may be a 
function of geometry and not of physical proper- 
ties. 

My experiments were experiments on the 
collapse by inward radial flow of hollow steel 
cylinders closed at both ends and exposed to 
hydrostatic pressure over the entire external 
surface. In all, 27 experiments were made on 27 
specimens of five different grades of steel. Each 
experiment consisted of exposing a virgin speci- 
men to a monotonically increasing pressure to a 
predetermined maximum, maintenance at the 
maximum for an interval, and monotonic release 
of pressure to zero. The cylinders were initially 
all of the same dimensions: 1.813’ long, 0.3125” 
o.d., and 0.0998” id. They were exposed to 
different maximum pressures varying from 71,000 
to 412,000 p.s.i. After exposure to pressure the 
permanent change of dimensions was determined. 
In general, pressure was increased to the maxi- 
mum in about 10 minutes, maintained for 10 
minutes, and released to zero in 2 or 3 minutes. 
In these limits there were no measureable time 
effects, as proved by control experiments in 
which pressure was immediately released after 
reaching the maximum, or maintained at the 


maximum for 2 hours. Of the 27 specimens, 5 
were of a 1020 steel left dead soft (designatéd as 
B in Figs. 3 and 4). The other 22 were all of the 
same 1045 steel in four different heat treatments. 
These were all heated to 1575°F. From this 
temperature a first set was quenched into salt at 
1100°F (designated as A5), a second set quenched 
into salt at 800°F (A6), a third set quenched into 
water and drawn at 800°F (A7), and a fourth 
set quenched into water and drawn at 1100°F 
(A8). The heat treatments were made at the 
Arsenal. 

The assumption of no permanent change of 
volume was checked from the measured change 
of dimensions. Actually a permanent decrease of 
density averaging around 0.2 percent was found. 
This is too small to invalidate the theoretical 
analysis; it is in agreement with other well-known 
results on the effect of cold work. The analysis 
neglects the elastic component of volume change. 
In the extreme case of any of the experiments 
the maximum component of radial displacement 
due to elastic distortion occurs at the external 
surface and is 7 percent of the displacement due 
to plastic flow. In all cases the elastic part 
vanishes in comparison with the plastic part at 
the inner surface. 





Fic. 1. The immediate 
neighborhood of the in- 
ternal cavity, under an 
original magnification of 
of 200, of the cylinder of 
1020 steel after exposure 
toa maximum collapsing 
pressure of 355,000 p.s.i. 
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The permanent changes of length varied be- 
tween extremes of +0.11 and —0.19 percent. 
This change of length is probably due to effects 
at the ends, where friction of the closing caps 
introduces some unavoidable local warping. In 
any event it would seem to be sufficiently close 
to the theoretical value zero to justify the 
analysis. 

The radial flow in 
samples was extreme, 


the case of the softer 
resulting in what to 
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Fic. 2. The immediate 
neighborhood of the in- 
ternal cavity, under an 
original magnification of 
200, of the cylinder of 
1020 steel after exposure 
toa maximum collapsing 
pressure of 412,000 p.s.i. 





superficial examination appeared as complete 
closing of the internal bore. The extremes to 
which plastic flow may be forced are shown in 
Figs. 1 and 2. These are microphotographs, made 
at the Arsenal, which show under a magnification 
of 200-fold the general region of the internal hole 
of the two specimens of 1020 steel exposed to 
maximum pressures of 355,000 and 412,000 p.s.i., 
respectively. These photographs show the enor- 
mous distortion in the original grains and make 


Fic. 3. The natural 
logarithm of the ratio of 
final external to final 
internal radius for the 
various steels described 
in the text as a function 
of “maximum collapsing 
pressure. 
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clear that even in the absence of other limiting 
factors the mathematical analysis cannot be 
pushed to extreme distortions because of failure 
of physical and geometrical isotropy and homo- 
geneity due to the grain structure of the material. 
In Fig. 3 are plotted the various values of the 
natural logarithm of the ratio of the final external 
to internal radius as a function of collapsing 
pressure. The initial value for the uncollapsed 
virgin specimen is 1.139, indicated in Fig. 3 
by a dotted line. It is seen that within experi- 
mental error, after plastic flow has once started, 
there is a linear relation between pressure and 
log(rext/rint). This holds into regions of deforma- 
tion so large that geometrical uniformity is lost 
and the granular structure begins to dominate. 
If the linear relation held indefinitely, the hole 
would become vanishingly small only for infinite 
pressure. This is contrary to some early experi- 
ments of mine! in which copper tubes were 
pushed to apparent complete collapse, and extra- 
polation indicated complete collapse of mild 
steel in the neighborhood of 320,000 p.s.i. It is 
now apparent that these results, which were 
unexpected because they indicated a strain 
softening instead of a strain hardening for ex- 
treme distortions, are to be explained in terms 
of a failure of geometrical regularity due to 
increasing dominance by the individual grains. 
The theoretical equations give the means for 
deducing the form of the curves of Fig. 3 in 
either of two cases: (a) There is strain hardening 





1P. W. Bridgman, Phys. Rev. 34, 1 (1912). 
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and the general relation between stress and 
strain is determined by some such criterion as a 
universal functional relation between ‘“‘signifi- 
cant” stress and “significant’”’ strain; (b) there 
is no strain hardening and flow is determined 
either by the maximum shearing stress criterion 
or the maximum shearing energy criterion, which 
under these particular conditions become identi- 
cal. In the latter case the result is particularly 
simple, and is P/log(rext/r int) =const., where the 
constant is the critical shearing stress. In Fig. 4 
P/log(rext/rint) is plotted as a function of P; it 
is not constant, but in all cases rises continuously 
with increasing pressure. Hence, under these 
conditions the maximum shearing stress criterion 
is not valid, but there is strain hardening, which 
continues increasingly over the entire range, 
with no lapse into a strain softening. Theoreti- 
cally, it should be possible to work backwards 
from the curves of Fig. 4 and deduce the strain 
hardening curve, in conjunction with some such 
assumption as a functionalrelation between 
significant stress and strain. However, such a 
procedure turns out not to be sensitive, and 
only rough conclusions can be drawn about the 
strain hardening curve. It is definitely indicated 
that the hardening which would account for the 
curves of Figs. 3 and 4 increases with strain at 
a rate considerably less than the rate given by 
ordinary tension tests, probably less by a factor 
of 2 or 3. This is in line with other evidence 
which indicates that under large plastic defor- 
mations strain hardening in compression is less 
rapid than in tension. 
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Letters to the Editor 








“Physical Basis of Bird Navigation” 


JOSEPH SLEPIAN 
Westinghouse Research Laboratories, Pittsburgh 30, Pennsylvania 
January 6, 1948 


N his remarkable paper on the above subject in the 

Journal of Applied Physics,* Dr. Yeagley describes the 
astonishing experimental support to his hypothesis that 
navigating birds are ‘sensitive’ to the magnetic field and 
Coriolis forces associated with the earth. Astonishing as 
this hypothesis may seem, I do not know of any funda- 
mental law of physics which is contradicted by it, or which 
would rule it out, as thus stated. 

However, Dr. Yeagley, in his paper, states his hypothesis 
in a more limited manner, and to this more limited hy- 
pothesis, it seems to me, serious objection can be made. 
Dr. Yeagley states that the birds are “sensitive” to the 
effect of motion through the earth’s magnetic field (p. 1036, 
second column), or to the effect of flying through a mag- 
netic field (p. 1037, second column). In both cases, I 
presume Dr. Yeagley means a uniform motion through the 
magnetic field. 

The restricted principle of relativity applied to this case 
tells us that in a system, such as the flying bird, in uniform 
motion relative to the earth, electrical and magnetic 
phenomena will proceed in exactly the same way as if the 
system were at rest, that is, the bird on the ground, and the 
electric and magnetic field in the system changed by the 
following scheme: 


E,'=E,, H,'=H,, 
Vz Vz 
Rh, Hy+—Es 
c c 
E, = Sf. 4? Hi, = 7 y\ ’ 
c c (1) 
Vz Vz 
E.+—H, H.—“E, 
C,! =—— H,=———_— 


Here E., E,, Ez, and H,, Hy, Hz, are the components, 
respectively, of the earth’s electric and magnetic fields, vz, 
the velocity of the bird moving in the x direction, and c, the 
velocity of light. E.’, E,’, E.’, Hz’, H,’, H.’ are the com- 
ponents of the electric and magnetic fields to which the 
bird at rest on the earth must be subjected if it is to be 
“sensitive’’ to exactly the same effects as the bird moving 
with velocity vz through the fields E,, Ey, E., Hz, Hy, H:z. 

“If we neglect v,*/c? relative to (1), and also neglect the 
effect of motion through the earth’s electric field, that one of 
the above equations which is relevant to my argument 


becomes 


E,'-E, = —"H.. (2) 
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Hence any “sensitivity” which the bird may have to its 
uniform motion, vz, through the vertical component of the 
earth’s magnetic field, H,, will be indistinguishable from its 
“sensitivity’’ to a change in the horizontal component of 
the earth’s electric field given by (2). 

Unfortunately, the intensity of the earth’s electric field is 
normally thousands of times larger than the quantity 
vz/cH, in Eq. (2). If the earth’s surface were perfectly 
smooth and horizontal, and if there were no clouds or other 
disturbing charge distributions in the air, we might argue 
that the horizontal component of the earth’s electric field, 
E,, remains always zero. But with actual terrain, local 
winds, and clouds, E, will always be large, and variable. 
Hence, we may conclude that whatever “sensitivity” the 
bird may have to its motion, v;, in the field, H,, it will be 
quite completely overshadowed and obliterated by the 
indistinguishable ‘‘sensitivity’’ to changes in the earth's 
horizontal electric field, Ey. 

I have frequently had occasion in the past to use the 
above argument to convince inventors who aspire to 
measure the velocity of an airplane by observing the effects 
of its motion through the earth's magnetic field upon 
suitable electrical apparatus on the plane. No matter how 
ingeniously they devise their apparatus, the effect of 
changes in the earth’s horizontal field is indistinguishable 
from, and overshadows, the effect they are trying to 
observe. 

It seems to me, then, that Dr. Yeagley will need to 
postulate a “‘sensitivity’’ of the navigating bird to the 
magnetic field directly, rather than a “‘sensitivity’’ to the 
effect of the uniform motion of the bird through the field. 
He will also need to assign to the necessary observation, 
which the bird must have of the ground over which the 
bird is flying, a role more primary than that of merely co- 
operating with the “‘sensitivity” to its motion through the 
magnetic field (p. 1037, second column). 


* Henry L. Yeagley, J. App. Phys. 18, 1035 (1947). 





Remarks on: “‘A Preliminary Study of a Physical 
Basis of Bird Navigation’’* 


RUSSELL H. VARIAN 


Microwave Laboratory, Department of Physics, Stanford University, 
Stanford, California 


January 22, 1948 


O support his very interesting theory of the mecha- 
nism of bird navigation, the author presents very 
strong experimental evidence for the conclusion that 
pigeons, and presumably other birds, find their way home 
over territory which is unfamiliar, by the aid of two senses. 
One is a sense of latitude, the other a sense of the strength 
of the vertical component of the earth’s magnetic field. 
Since the magnetic poles of the earth are displaced from the 
axis of rotation, the lines of equal vertical component 
intersect the latitude lines, forming a grid which may be 
used to define a location. 
The author appears to jump to unwarranted conclusions 
concerning the operation of these senses. 
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First, he postulates that the indication of the vertical 
component of the earth’s magnetic field is provided by 
some organ which can measure the e.m.f. produced in a 
conductor cutting the lines of force of the earth’s magnetic 
field. 

The primary reason for writing this letter is that the 
present writer, in his student days, conceived the idea of 
utilizing this effect to construct a ground speed meter for 
airplanes, and stirred up quite a bit of interest before he 
discovered that the method was fundamentally fallacious. 
Later, Dr. D. L. Webster produced a more general refuta- 
tion than that of the writer. This argument, which is based 
on relativity, is here outlined. 

In relativity, an identity exists between an electrostatic 
field and a magnetic field moving relative to a frame of 
reference. It then follows that any attempt to shield a pick- 
up system from spurious electrostatic forces must shield the 
pick-up from the influence of the moving magnetic field. 
This is briefly treated by W. R. Smythe, Static and Dynamic 
Electricity (McGraw-Hill Book Company, Inc., New York, 
1939), pp. 499-501. If we dispense with shielding, we are 
then confronted with the conductor. This reduces itself 
to measuring a potential difference between the conductor 
ends and the atmosphere, but the atmosphere is moving 
with respect to the magnetic field and therefore has an 
electrostatic field caused by its motion. Thus, if the 
formidable obstacle of extraneous electrostatic forces is 
successfully overcome, the end result is a measurement 
of air speed. Since pigeons cannot measure the vertical 
component of the earth’s field by measuring the rate of 
cutting magnetic lines of force, the alternative appears to 
be that they possess an organ capable of acting as a 
magnetometer. 

The author’s second postulate is that the detection of 
latitude is made by measurement of the Coriolis force. This 
he computes for his case to be 0.0061 ft. per sec. per sec., or 
about 1/6000 g. In order for pigeons to locate themselves 
within 25 miles or so, a measurement of less than 1 percent 
of this would have to be made, or very roughly a measure- 
ment of 10~* g. The only way the acceleration of gravity 
may be removed from the measurement is by a knowledge 
of its direction of action to about 10~* radian without 
reference to the force of g for determining that direction. 
This appears to be practically impossible, even though it 
may be theoretically possible. It would therefore appear 
that as far as the evidence presented in the article is con- 
cerned, the determination of latitude by the sun is more 
likely than determination by Coriolis force. 

The hypothesis that a bird measures latitude by the 
elevation of the sun should receive attention. The sun 
appears as a fairly large object, yet determination of its 
elevation to one solar diameter would give latitude within 
30 miles; thus, it appears that if solar elevation is used, no 
great precision is required to determine latitude within the 
pigeon’s apparent accuracy. The cited fact that pigeons are 
unable to navigate in thick haze or fog or complete darkness 
is concordant with this view. 

It should be further pointed out that either mutation or 
chance development of an organ for another function must 
produce an organ capable of measuring Coriolis force suffi- 
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ciently well to be of survival value before evolution can 
take hold to produce further refinement. This seems to be 
almost an impossible occurrence. ° 

The following experiments to throw further light on the 
nature of the organs involved suggest themselves: 

1. Attach magnets to the birds where they will be 
stationary with respect to the body of the bird, but will 
produce a change in vertical component of magnetic field. 
2. Determine definitely whether pigeons can navigate 
under: conditions of good visibility but heavy overcast. 
This last experiment would be conclusive if it showed 
pigeons could not navigate under heavy overcast, but if the 
result were the reverse it would not rule out an organ 
sensitive to extreme infra-red which could operate through 
considerable overcast. 


* Henry L. Yeagley, J. App. Phys. 18, 1035 (1947). 





Remarks on: ‘‘The Physical Basis of 
Bird Navigation” 


LEVERETT Davis, JR. 
California Institute of Technology, Pasadena, California 
January 19, 1948 


N a recent article! H. L. Yeagley describes experiments 
with homing pigeons which he interprets as indicating 
that homing ability is based on the detection of induced 
electromotive forces due to motion with respect to the 
earth’s magnetic field and on the detection of the Coriolis 
force due to motion with respect to the rotating earth. 
However, even if it is assumed that a pigeon can make any 
required measurements with infinite accuracy, it seems 
conceptually impossible for it to detect either the electro- 
motive force or the Coriolis force. A discussion of the 
impossibility of determining the electromotive force is 
given by Smythe,? who bases his discussion on the well 
established principal that the electric fields observed by 
virtue of motion through a uniform magnetic field are 
indistinguishable from uniform static fields. Since the static 
fields due to atmospheric electricity’ are of the order of one 
volt/cm and are quite variable, it would appear impossible 
to detect a superposed field of 10-5 volt/cm due to motion. 
It is also impossible that these fields should produce 
physiological effects, as suggested by Yeagley. The situa- 
tion with respect to Coriolis force is even worse. It is 
difficult to conceive of a device that responds to Coriolis 
force that does not respond equally to gravitational force, 
certainly no mechanical device will. Hence it is impossible 
to separate the Coriolis force from the force of gravity 
except by comparison of measurements of two observers 
having different velocities, and hence different Coriolis 
forces. It does not help the pigeon to change its speed unless 
it can, in effect, determine the position of the vertical when 
stationary on the ground and then recognize the position of 
this line by visual observations (or by the use of a gyro- 
scope), to about 0.1 second of arc for comparison with the 
apparent vertical when flying. When flying over water or 
fog where the solid earth is not visible, it would seem to be 
impossible for the pigeon to determine either its speed or 
acceleration with respect to the earth, both being needed 
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for a determination of latitude by means of the Coriolis 
force, and the speed being needed for a determination of the 
vertical component of the magnetic field from the induced 
electromotive force. 

Even though it is impossible for the pigeon to detect the 
induced electromotive force produced by flying through the 
earth’s field, and thereby determine the vertical component 
of the field, there is no such reason to suppose that the bird 
cannot measure the field directly, perhaps by means of the 
e.m.f.’s induced in moving blood. And there are many ways 
whereby latitude can be measured without measuring 
Coriolis force. Offhand, it would seem easier for the pigeon 
to observe the sun or the stars or to measure the accelera- 
tion of gravity, which depends on latitude, than to measure 
a force of the order of the Coriolis force. It is evident, there- 
fore, that Yeagley can meet the above objections without 
an essential change in his theory merely by adopting some 
other mechanisms by which a pigeon measures the vertical 
component of the magnetic field and the latitude. It is 
obvious that this would not affect any of the experiments he 
describes. 

Because it seems desirable to test such a remarkable 
theory in as many ways as possible in spite of the amazing 
experiments already reported, those of the following experi- 
ments that have not yet been tried might well be considered. 
To distinguish between Yeagley’s assumption that pigeons 
are sensitive to electric fields, the assumption that they are 
sensitive to magnetic fields, and the assumption that they 
are sensitive to neither, the most direct experiment would 
be a maze experiment or a feeding apparatus experiment in 
which the essential clues that the bird would have would be 
appropriate fields. Similar experiments in slowly ac- 
celerating elevators or rotating houses might show whether 
pigeons were sensitive to small accelerations and hence that 
such a sensitivity could be used for the determination of 
latitude. It would be interesting to know the effect on 
homing ability of a magnet fastened to the bird’s body 
rather than to its wings and of the effect of magnetic 
storms. Experiments, such as the Special Release No. 1 of 
Yeagley’s First Nebraska Experiment, in which birds flew 
toward a conjugate point to which they had never been 
would seem to be somewhat freer of unexpected alternative 
explanations than experiments in which the pigeons knew 
their loft was at the conjugate point because they had seen 
it there. 

1H. L. Yeagley, J. App. Phys. 18, 1035 (1947). 

2W. R. Smythe Static and Dynamic Electricity (McGraw-Hill Book 


Company, Inc., New York, 1939), §14.14. 
3W. F. G. Swann, Int. Crit. Tab. 6, 442. 





A New Method of Measuring Diffusion 
Coefficient in Solids with 
Radioactive Tracers 


G. C. KuczyNskI 
Metallurgical Research, Sylvania Center, Bayside, New York 
January 16, 1948 


WO experimental methods have been used up to the 
present time in the determination of the diffusion 
coefficient with the help of radioactive tracers. The older 
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one, as described by W. A. Johnson,' requires the slicing of 
a specimen previously heated to a certain temperature for a 
long period of time. From the dissolved slices the average 
concentration of the radioactive material is determined and 
the concentration versus depth curve is constructed. From 
this curve, coefficient of diffusion D can readily be de- 
termined. The other method (2) dispenses of cutting, and 
from one reading of radioactive intensity and knowledge of 
the absorption coefficient for a given material the diffusion 
coefficient can be determined. In both methods the appro- 
priate correction has to be made for radioactive decomposi- 
tion during heating time. 

The present method does not require cutting of the 
specimens, nor the knowledge of the absorption coefficient. 

One face of a rectangular specimen a few centimeters in 
length is covered with a thin layér of radioactive material 
whose diffusion is to be measured and is then heated for a 
long period of time. Afterwards, the flat surface perpen- 
dicular to the activated face is etched electrolytically. The 
specimen is then transferred to a lead box (Fig. 1) with its 
etched surface up. The sliding cover is moved so that edge 
A is at an arbitrary distance x; from the radioactive layer 
and a reading of the radioactivity J; is taken. Then the 
sliding cover is moved to the left so that edge B will take 
the position x2 such that the radioactivity reading J: from 
this portion of the specimen is equal J;. From these two 
readings only, the diffusion coefficient can be calculated as 
follows: 


o 


ez)/2(Dt) 
I, =bLeart I, e~* dz, 


where a is a factor expressing absorption, b the width of the 
specimen, and ig the original concentration of the radio- 
active material on the activated face. 
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From J, = J, denoting 
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Q 
one obtains 
o(¥1) +o(y2) = 3, 
where 
yr=x1/2(Dt)t. and = ye=x2/2(Dt). 
Since 
2—, o(n)+¢ *y,) = }. 
MyM M1 Xi 
This equation can easily be solved by graphical methods. 
When 4; is known, D is calculated from 


D=(1/t)(x1/2y1)*. 


It can be seen that the determination of D does not depend 
upon knowledge of the absorption coefficient and original 
concentration of the radioactive material on the activated 
face. 

This method is expected to produce results under the 
following conditions: 1. The radioactive layer should be 
thin. 2. The surface from which the readings are taken 
should be deeply and uniformly etched to remove all 
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radioactive material which has spread too far because of 
the possibly accelerated surface diffusion. 3. The time of 
heating should be long in order to obtain a deep penetration 
measurable in terms of x; and x2. 4. x; and x2 should be 
measured with great accuracy (0.001 cm or better). 5. The 
length of the specimen should be great in comparison with 
x, and X2. 

Arrangements are being made to test this material ex- 
perimentally in Sylvania Bayside laboratories. 


W. A. Johnson, Trans. A.I.M.M.E. 143, 107 (1941). 
J. Steigman, W. Shockley, and F. C. Nix, Phys. Rev. 56, 13 (1939) 





New Booklets 








The Esterline-Angus Company, Inc., P. O. Box 596, 
Indianapolis 6, Indiana, recently published technical bulle- 
tin No. 247, entitled Operation Recorders—Their Selection 
and Use. A lengthy list of applications is included. 16 
pages, free on request. 


Burrell Technical Supply Company, 1942 Fifth Avenue, 
Pittsburgh 19, Pennsylvania, has issued its bulletin No. 
207, describing the Burrell ‘‘Wrist-Action’’ Shaker, for use 
in all types of laboratories. This device simulates hand 
shaking. 4 pages, free on request. 


The Apparatus Department of General Electric Com- 
pany, Schenectady 5, New York, has published a 48-page 
brochure entitled Specialized Testing and Measuring Equip- 
ment Catalog. 


Tube Department of Radio Corporation of America, 
Harrison, New Jersey, has issued the following Application 
Notes: Suppression of Arc-Over, Corona, and High-Voltage 
Leakage in the 5TP4 Kinescope (4 pages) and Blower Re- 
quirements for RCA Forced Air-Cooled Tubes (12 pages). 


Brooks and Perkins, Inc., magnesium fabricators of 
Detroit 1, Michigan, issued in December the first issue of 
their selling magazine, Magnesium. The issue contained 12 
pages and included several articles on magnesium and its 
place in the manufacturing picture. Illustrated with 
photographs. 


VOLUME 19, MARCH, 1948 


Armour Research Foundation of Illinois Institute of 
Technology, Chicago 16, Illinois, has published a 16-page 
booklet entitled Precision Gage Laboratory—A nother A.R.F, 
Service to Industry. Free on request. 


Fischer and Porter Company, Department 6F-H, 
Hatboro, Pennsylvania, announces its new bulletin, No. 88, 
containing more detailed information on the company’s 
Viscorator for viscosity measurement. 4 pages, free on 
request. 


North American Philips Company, Inc., 100 East 42nd 
Street, New York 17, New York, is distributing Volume 9, 
Number 6 of Philips Technical Review, published in Holland. 
Following are the titles of articles contained in the issue: 
“A 48-channel carrier telephone system,”’ “‘An improved 
method for the air-cooling of transmitting valves,” “Im- 
provements in the construction of cathode-ray tubes,’’ and 
“On the crystalline structure of ferrites and analogous 
metal oxides.” 





Here and There 


Personnel 








Two appointments were recently announced by the 
National Bureau of Standards. Dr. Wallace R. Brode, 
chemist, has been named Associate Director of the Bureau, 
succeeding Dr. Hugh L. Dryden, who was recently ap- 
pointed Director of Aeronautical Research of the National 
Advisory Committee for Aeronautics. Dr. Brode was 
formerly at Ohio State University. John Todd, British 
mathematician, has also been appointed to the staff. On 
leave of absence from Kings College, University of London, 
he will aid in the development and application of high speed 
automatic digital computing machines. 


Dr. Samuel E. Sheppard, widely known for his research 
on the sensitivity of photographic materials, retired 
January 1 after 35 years with Kodak Research Laboratories. 


Conference on Surface Reactions 


The Pittsburgh International Conference on Surface 
Reactions will be held at the Mellon Institute for Industrial 
Research in Pittsburgh, June 7-11, 1948. The Conference 
Committee is led by Dr. Earl Gulbransen of Westinghouse 
Research Laboratories. Preliminary plans call for technical 
sessions mornings and evenings and, in the afternoons, 
visits to Pittsburgh research laboratories which are working 
on surface reactions. Scientists, engineers, and educators 
from many parts of the world have been invited to par- 
ticipate, and several will present papers. 


General Electric Fellowships 


The General Electric Company is again offering (for the 
summer of 1948) fellowships for a six-weeks course at 
Union College and at Case Institute of Technology to high 
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